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Introduction 


The  present  study  is  investigating  the  basic  neuromotor  coordination 
mechanisms  involved  in  a  rapid  elbow  flexion  movement,  and  in  developing 
mathematical  models  to  explain  the  interaction  of  these  basic  neuromotor 
coordinations  with  the  biomechanical  parameters  of  movement  speed.  Speed  of 
movement  is  being  assessed  in  a  biomechanics  mode  via  displacement,  velocity, 
acceleration,  point  of  inflection  between  acceleration  and  deceleration,  and 
the  total  time  of  an  elbow  flexion  movement.  Electromyographic  analysis 
techniques  are  used  to  monitor  the  sequential  timing  of  agonist  and  antagonist 
muscle  activity.  The  mathematical  modelling  effort  incorporates  the  biomechanical 
parameters  into  an  interface  with  the  neurophysiological  parameters  involving 
the  central  and  peripheral  nervous  systems,  and  then  extends  the  interface  to 
include  viscicelastic  perperties  of  the  muscle,  activation  delays,  and  neuronal 
pools . 

The  experimental  approach  includes  consideration  of  the  neuromotor  coordination 
mechanisms  in  both  loaded  and  unloaded  elbow  flexion  movements;  changes  in 
control  mechanisms  due  to  practice  and  learning  effects;  changes  in  control 
mechanisms  due  to  local  muscular  fatigue  induced  by  isometric  exercise  in  the 
agonist  and  in  the  antagonist  muscle  groups;  feasible  training  and  practice 
regimens  involving  artificial  means  of  enhancing  beneficial  changes  in  control 
mechanisms;  and  development  of  suitable  mathematical  models  to  explain  in 
functional  terms  the  ways  in  which  the  adaptive  mechanisms  can  account  for 
changes  in  basic  coordination  skill  and  the  breakdown  of  skilled  movement 
patterns  due  to  local  muscular  fatigue.  The  planned  series  of  studies  in¬ 
corporates  research  protocols  from  at  least  three  usually  distinct  and  isolated 
areas  of  research:  neurophysiology,  biomechanics,  and  computer  science  and 
mathematical  modelling. 


Summary  of  Part  I.  Experimental  Studies 

The  first  formal  investigation  completed  December,  1980  was  designed 
to  secure  criterion  measure  estimates  essential  for  input  into  preliminary 
mathematical  modelling  attempts  for  an  unloaded  fast  elbow  flexion  movement. 

As  the  early  mathematical  models  developed,  it  was  necessary  to  revise 
experimental  methodologies  and  to  collect  additional  criterion  measures  as 
dictated  by  modelling  needs.  The  first  stage  of  mathematical  modelling  has 
been  completed  and  is  described  in  Part  II  of  this  report. 

Beginning  in  January,  1981  the  second  formal  investigation  was  begun. 

Again  guided  by  the  insights  and  requirements  generated  by  the  mathematical 
modelling,  the  second  study  assessed  basic  neuromotor  mechanisms  involved  in 
fast  elbow  flexion:  (a)  in  loaded  and  'unloaded  conditions,  (b)  under 
conditions  of  rested  and  induced  fatigue  of  the  agonist  or  the  antagonist 
muscle  groups,  (c)  a  unique  condition  in  which  agonist  and  antagonist  muscle 
groups  are,  simultaneously  fatigued,  and  (d)  in  male  and  female  subjects. 

During  the  second  year  of  the  grant,  data  were  collected  on  make  and 
female  subjects  trained  on  isotonic  rapid  elbow  flexion  movements.  Following 
sufficient  training  to  yield  reproducible  motion  characteristics,  muscle  fatigue 
was  induced  separately  for  agonist  and  antagonist  muscles  as  well  as  for  a  combined 
agonist  and  antagonist  muscle  fatigue  regimen.  Both  high  and  low  intensity  loads 
were  employed  to  produce  isotonic  fatigue.  Testing  was  done  in  the  loaded  and 
unloaded  modes.  Biomechanical  data  on  displacement,  velocity,  acceleration, 
and  point  of  inflection  between  acceleration  and  deceleration  was  obtained. 
Concurrent  EMG  activity  of  agonist  and  antagonist  muscle  groups  was  monitored. 


The  results  are  being  used  to  further  refine  mathematical  models  of  fast  arm 


3. 


L -  A. V  ■*  "  «.■*  "  J.*”  ""  ”  ■"  "  K  v "  "  ■ 

% 


movement  as  driven  by  volitional  commands  that  were  derived  during  the  first 
year  of  the  grant.  Pilot  studies  were  initiated  to  evaluate  the  feasibility 
of  using  functional  electrical  stimulation  to  increase  speed  of  limb  movement, 
aid  in  identifying  mechanisms  of  neuromotor  coordination,  and  examine  the 
likelihood  of  imparting  movement  patterns  artificially  without  actual  physical 
practice. 


<r 


Part  I— A.  Completed  Work 

Both  Principals  Investigators  devoted  full  time  during  the  first  summer  to 
the  grant  by  (a)  identifying,  reviewing,  and  evaluating  recent  research 
publications  relevant  to  theoretical  aspects  of  the  grant  project,  (b)  planning 
the  revised  research  protocols  in  detailed  methodology,  and  (c)  identifying, 
reviewing,  and  evaluating  recent  research  publications  relevant  to  data  analysis 
and  computer  programming.  Two  research  assistants  were  employed  during  the 
month  of  August  to  expedite  these  undertakings.  Preliminary  orders  for  neces¬ 
sary  expendable  supplies  were  sent  out. 

In  September,  1380  additional  graduate  research  assistants  (RA's)  were 
recruited.  In* Dr.  Kilmer's  case  the  RA  continued  on  with  computer  programming 
for  the  mathematical  modelling  aspects  of  the  project.  In  Dr.  Kroll's  case  new 
RA's  were  indoctrinated  into  research  data  collection  protocols.  The  in¬ 
doctrination  phase  included  practice  in  apparatus  usage,  maintenance,  and 
calibration  techniques  as  well  as  pilot  testing  research  measurement  schedules 
to  insure  reliability  of  data  collection.  The  indoctrination  phase  lasted  six 
weeks  resulting  in  highly  proficient  research  assistants  capable  cf  collecting 


reliable  data. 


Beginning  the  end  of  October,  1981  data  collection  for  the  first  project 
began.  The  first  study  was  designed  to  secure  essential  criterion  measure 
necessary  for  input  into  preliminary  mathematical  models.  These  criterion 
measures  include: 

...maximum  isometric  strength  of  elbow  flexion  and  extension 

...rate  of  tension  development 

...integrated  electromyographic  activity  at  maximal  and  sub-maximal 
tensions  and  joint  angles 

...agonist-antagonist  contraction  patterns  during  unloaded  speed-of- 
movement  elbow  flexion 

...effects  of  local  muscular  fatigue  upon  neuromotor  coordination 

Preliminary  mathematical  modelling  attempts  quickly  showed  that  additional 
criterion  measures  would  be  necessary.  As  a  result  of  the  modelling  requirements 
upper  limb  volumes  and  biceps  and  triceps  reflex  time  and  rate  of  muscle  tension 
development  were  collected.  Limb  volumes  were  estimated  by  anthropometric  and 
water  displacement  techniques.  The  research  assistant  was  trained  in  these 
techniques  by  Professor  Frank  Katch  of  the  Department  of  Exercise  Science  who  is 
a  recognized  researcher  in  the  body  composition  area.  A  MEDLARS  literature 
search  revealed  only  a  few  scattered  articles  dealing  with  upper  limb  reflex 
times  with  none  of  the  data  appropriate  for  modelling  requirements.  Suitable 
apparatus  was  designed  and  built  for  upper  human  limb  reflex  testing.  Limb 
volumes  and  reflex  times  have  been  collected  on  all  subjects  tested  to  date. 

In  this  first  study  of  Phase  I,  each  of  the  12  subjects  was  tested  over 
ten  experimental  sessions  with  each  session  taking  90-120  minutes.  Sight  of 
these  sessions  involved  speed  of  movement  testing  while  the  other  two  sessions 
involved  upper  limb  reflex  and  volume  assessments.  Each  of  the  first  four 
experimental  days  included  50  trials  of  rapid  elbow  flexion  speed  of  movement 


to  a  designated  target.  Previous  research  has  shown  that  this  amount  of 
practice  insures  that  all  subjects  become  well-practiced  and  exhibit  consistent, 
stabilized  performance  both  in  the  speed  of  movement  as  well  as  in  the  neuro- 
motor  coordination  mechanisms.  Following  the  establishment  of  well-practiced 
performance,  local  muscular  fatigue  was  induced  by  two  different  effort/rest 
exercise  regimens.  A  5-seccnd  maximal  isometric  contraction  with  a  5-second 
intertrial  rest  period — designated  the  5:5  condition — for  a  total  of  30  serial 
trials  constituted  the  more  intense  fatigue  series.  A  5-second  maximal  isometric 
contraction  with  a  10-second  intertrial  rest  period — designated  the  5:10 
condition — for  a  total  of  30  serial  trials  constituted  the  less  intense  exercise 
regimen. 

These  two  exercise  regimens — the  5:5  and  the  5:10  conditions — were  designed 
to  produce  different  levels  of  fatigue  in  the  involved  muscle  group.  The  5:5 
and  the  5:10  exercise  regimens  were  administered  on  separate  occasions  to  the 
agonist  (elbow  flexors)  and  to  the  antagonist  (elbow  extensors)  muscle  groups 
using  a  balanced  order  order  of  presentation  over  subjects  and  across  test 
sessions  to  minimize  contaminating  test  order  effects. 

One  of  the  purposes  of  inducing  high  and  low  intensity  local  muscular 
fatigue  in  the  agonist  or  antagonist  muscle  groups  was  to  ascertain  the  role  of 
changes  in  the  peripheral  muscle  state  and  in  peripheral  muscle  afferent 
feedback  to  higher  nervous  system  centers.  Specifically,  the  question  of  whether 
or  not  programmed  central  commands  for  a  fast  ballistic  movement  can  be  altered 
by  the  presence  of  different  levels  of  muscular  fatigue. 

1-3 .  Work  in  Progress 

In  Phase  II  of  the  grant  proposal  an  investigation  using  the  same  testing 
schedule  outlined  for  the  Phase  I  study  has  been  replicactd  for  elbow  flexion 
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speed  of  movement.  However,  to  further  elucidate  the  basic  neuromotor 
coordination  mechanisms  a  load  is  imposed  upon  the  limb  movement  task.  The 
imposition  of  a  lead,  of  course,  greatly  affects  the  muscle  activation  time 
and  sequential  firing  of  involved  muscle  groups  and  constitutes  a  more  complex 
movement  task  with  distinctly  more  complex  nervous  system  control  system 
involvement. 

In  order  to  insure  identical  loading  for  subjects  with  different  limb 
lengths,  an  inertial  loading  technique  was  used.  The  natural  moment  of 
inertia  can  be  altered  by  placement  of  light  weights  at  long  distances  and 
heavy  weights  at  short  distances  from  the  fulcrum  point.  The  protocol  employed 
utilized  a  fixed  weight  for  all  subjects  but  varied  the  distance  from  the 
fulcrum  point  to  comply  with  the  condition  of  imposed  moment  of  inertia.  One 
load  condition  used  was  equal  to  two  times  the  natural  moment  of  inertia  while 
the  second  load  condition  was  equal  to  five  times  the  natural  moment  of  inertia. 

In  addition  to  the  two  different  exercise  regimens  of  5:5  and  5:10 
administered  to  the  agonist  and  to  the  antagonist  muscle  groups  on  separate 
occasions,  another  set  of  exercise  regimens  was  also  administered.  Based  upon 
the  mathematical  modelling  considerations,  it  was  of  interest  to  ascertain  the 
effect  of  local  muscular  fatigue  in  both  the  agonist  and  antagonist  muscle 
groups  simultaneously  upon  basic  neuromotor  coordination  mechanisms.  To 
accomplish  such  a  goal  of  producing  fatigue  in  both  agonist  and  antagonist 
muscle  groups,  a  f lexion-to-extension  contraction  sequence  was  designed. 

3y  using  a  5-second  contraction  of  the  elbow  flexors  followed  immediately 
by  a  5-second  contraction  of  the  elbow  extensors  with  no  intertrial  rest  period, 
fatigue  effects  equivalent  to  5-secor.d  contraction  and  5-second  rest  periods 


were  produced  in  both  the  agonist  and  antagonist  muscle  groups.  Similarly, 
if  the  effort  sequence  is  a  5-second  contraction  of  the  elbow  flexors 
followed  by  a  5-second  contraction  of  the  elbow  extensors  followed  by  a  5-second 
rest  period,  the  equivalent  fatigue  effects  of  a  5-second  contraction,  10-second 
rest  period  exercise  regimen  can  be  produced  in  both  the  agonist  and  antagonist 
muscle  groups  simultaneously. 

These  two  new  exercise  regimens,  designated  the  5/5:0  and  the  5/5:5, 
aloow  assessment  of  fatigue  effects  in  the  agonist  and  the  antagonist  muscle 
groups  upon  basic  neuromotor  coordination  mechanisms  and  provide  a  stringent 
test  of  the  predictive  power  of  the  mathematical  model  being  developed.  In 
actuality  the  two  new  exercise  regimens  produced  slightly  more  fatigue  than 
the  5:5  and  the  5:10  regimens  upon  which  they  were  based  because  some  degree 
of  co-contraction  of  agonist  and  antagonist  muscle  groups  occurs  in  maximal 
isometric  contractions.  Such  co-contraction  effects,  however,  are  desirable 
since  they  will  heighten  the  local  muscular  fatigue  produced. 

The  Phase  II  study  described  above  involved  12  male  and  12  female  subjects 
with  each  subject  attending  12  test  sessions  of  90-120  minutes  each.  Ten  of 
these  sessions  involved  speed  of  movement  and  exercise  conditions  while  the 
other  two  sessions  were  for  upper  limb  reflex  and  volume  testing.  One  of  the 
research  assistants  (Zulma  C.  Garcia)  used  some  of  the  data  for  her  Ph.D. 
dissertation  in  exercise  science  under  the  chairmanship  of  PI  Kroll.  Data 
collection  was  completed  March,  1981  with  subsequent  reduction  and  analysis  of 
data  taking  place.  A  Master's  thesis  was  completed  as  a  pilot  study  dealing 
with  the  effects  of  vibration  upon  speed  of  movement  and  neuromotor  coordination 
mechanisms.  This  graduate  student,  Marilyn  Teves,  was  not  a  research  assistant 
but  her  work  was  supported  in  part  by  purchase  of  essential  expendable  supplies . 


Data  collected  during  Year  1  (1  July  1980-1  June  1981)  was  reduced, 
analyzed,  and  results  inputted  into  mathematical  modelling  projects.  Research 
manuscripts  suitable  for  publication  in  scientific  journals  have  been  prepared 
(see  Appendices).  Based  upon  insights  generated  by  the  mathematical  modelling, 
refinements  in  studies  planned  for  Year  2  became  necessary  as  they  were  also 
during  Year  1.  Independent  of  the  analysis  of  data  all  ready  collected  and 
being  inputted  into  mathematical  modelling  projects,  Year  2  experimental  studies 
were  begun  the  summer  of  1981. 

Essentially  the  same  testing  schedule  which  replicated  fatigue  and 
coordination  studies  during  Year  1  was  begun  the  summer  of  1981.  The 
exercise  regimens,  however,  utilized  isotonic  rather  than  isometric  muscular 
contractions  to  induce  fatigue. 

Testing  Schedule 

Unloaded  Movement  Coordination  and  Fatigue 

Subjects  will  perform  50  trials  daily  of  rapid  elbow  flexion  speed  of 
movement.  Four  days  of  practice  will  be  necessary  to  produce  consistent  and 
stabilized  performance.  Following  establishment  of  stabilized  performance  of 
the  neuromotor  control  mechanisms  and  motion  characteristics ,  local  muscular 
fatigue  will  be  induced  as  per  the  following  exercise  regimens: 

1.  Agonist  Muscle  Fatigue 

Six  bouts  of  isotonic  elbow  flexion.  During  each  bout  subjects 
will  perform  repetitions  at  a  load  equal  to  that  necessary  to 
produce  a  six  repetition  maximum  (S-RM)  or  a  20-repetition 
maximum  (20-RM)  bout.  Preliminary  testing  on  subject  during  the 
baseline  stabilization  days  will  be  used  to  establish  the 
individual  loads  to  be  used  to  produce  a  high  intensity  (S-RI!)  and 


low  intensity  (20-RM)  fatiguing  exercising  regimens. 
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2.  Antagonist  Muscle  Fatigue 
Same  protocol  as  1  above. 

One  of  the  purposes  of  high  and  low  intensity  local  muscular  fatigue  via 
the  isotonic  exercise  regimens  on  the  agonist  and  antagonist  muscle  groups 
separately  is  to  ascertain  the  role  of  changes  in  peripheral  muscle  feedback  to 
higher  nervous  system  centers  and  determine  if  programmed  central  commands  for 
a  fast  ballistic  movement  can  be  altered  by  fatigue.  Isotonic  and  isometric 
exercise  regimens  produce  different  kinds  of  fatigue  and  comparisons  of  the 
isotonic  and  isometric  exercise  regimen  effects  upon  neuromotor  coordination 
mechanisms  is  of  particular  relevance  to  understanding  coordination  of  a  fast 
movement. 

Loaded  Movement  Coordination  and  Fatigue 

The  same  testing  schedule  proposed  above  will  be  replicated  for  elbow 
flexion  speed  of  movement  against  a  load.  In  order  to  insure  identical 
loading  for  subjects  with  different  limb  lengths,  an  inertial  loading  technique 
will  be  used.  The  protocol  will  utilize  a  fixed  weight  for  all  subjects  and 
will  vary  the  distance  from  fulcrum  point  to  comply  with  the  condition  of  moment 
of  inertia.  One  load  condition  will  equal  two  times  the  nature  moment  of  inertia 
and  the  second  load  condition  will  equal  five  times  the  natural  moment  of  inertia 

Functional  Electrical  Stimulation 

Pilot  studies  were  initiated  to  evaluate  the  feasibility  of  using 
electrical  stimulation  to  increase  speed  of  fast  limb  movement  and  examine  the 
possibility  of  imparting  movement  patterns  artificially  without  actual  physical 


oractice . 


Study  of  the  literature  on  Funcational  Electrical  Stimulation  (FES) 


coupled  with  results  from  our  own  pilot  studies  showed  that  no  electrical 
stimulation  apparatus  was  commercially  available  which  possessed  all  of  the 
control  features  needed.  Extended  negotiations  with  several  manufacturers 
eventually  led  to  one  manufacturer  agreeing  to  make  necessary  modifications  to 
their  available  stimulator.  Design  problems  still  exist  with  the  present 
simulator  and  additional  modifications  are  being  made  by  our  own  electronic 
technician. 

Reverse  Look  Theory  of  Motor  Learning 

As  a  direct  result  of  Year  2  efforts,  a  new  theory  of  motor  learning  has 
been  formulated  anf  formally  presented  at  a  neurophysiology  symposium  on 
movement  disorders  held  at  Columbia  University  in  March  of  1982  (see  Appendix 
•  for  paper) .  The  new  theory  has  been  named  the  reverse  loop  theory  of  motor 

learning  since  it  is  based  upon  sensory  imparted  learning  (SIL)  produced  by 
electrical  stimulation  of  movement  synergy  patterns  prescribed  by  an  analysis 
1  of  neuromotor  coordination  mechanisms  of  human  limb  movement.  One  article  is 

being  prepared  describing  the  reverse  look  theory  oi  motor  learning  and  sensory 
imparted  learning  for  publication  in  a  suitable  journal,  possible  the  Journal 
of  Motor  Behavior  of  Experimental  Brain  Research. 

II.  Proposed  Experimental  Studies  -  Year  3 

At  the  present  time  we  are  listing  the  effectiveness  of  several  aspects 
of  electrical  stimulation. 

1.  Wave  form:  monophasic ,  biphasic,  biphasic  compensated 

2.  Pulse  duration 

3.  Frequency  of  pulses 


4.  Intensity  level 
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The  problem  is  made  more  complex  because  of  the  likelihood  that  motor  units 
of  different  muscle  fiber  type  composition  may  require  different  stimulus  train 
parameters.  Current  efforts  in  mathematical  modelling  are  directed  at  this 
problem  and  will  provide  important  guidelines  for  stimulation  train  parameters. 
Even  at  this  time,  however,  crude  parameter  values  have  been  ascertained  for 
effective  stimulation  and  continued  work  will  serve  to  refine  the  techniques 
and  no  obstacle  to  further  progress  is  anticipated. 

Early  pilot  testing  revealed  that  electrical  stimulation  to  produce  muscle 
synergy  contraction  as  prescribed  by  a  model  of  neuromotor  coordination 
mechanises  can  produce  learning  (see  appendix  for  paper  presented  at  Columbia 
University  Neurophysiology  Symposium  on  movement  disorders) .  Studies  proposed 
for  Year  3  will  examine  the  importance  of  neuromotor  coordination  mechanisms  as 
far  as  producing  improved  performance  through  the  use  of  electrical  stimulation. 

Our  work  to  date  has  shown  several  important  neuromotor  coordination 
mechanisms  are  involved  in  fast  limb  movement: 

1.  agonist  to  antagonist  latency  of  activation 

2.  duration  of  EMG  burst  (or  envelope  of  EMG  activity)  in  agonist 
and  antagonist  modes 

3.  relative  EMG  intensity  of  agonist  and  antagonist  muscle  activity 

Series  1.  A  series  of  activities  will  manipulate  all  three  basic  neuromotor 

coordination  mechanisms  simultaneously  and  assess  the  degree  to  which  well- 
practiced  speed  of  limb  movement  can  be  improved  by  elec-trical  stimulation 
designed  to  simulate  muscle  synergy  contraction  patterns.  3y  optimizing  the 
three  neuromotor  coordination  mechanisms  to  produce  faster  limb  movement  the 
effectiveness  and  magnitude  of  electrical  stimulation  regimen  can  be  assessed. 

Series  2.  In  order  to  elucidate  the  relative  imoortance  of  each  of  these 


three  major  neuromotor  coordination  mechanisms  to  limb  movement,  electrical 


stimulation  regimens  must  be  devised  to  produce  sensory  imparted  learning 
which  manipulates  each  of  these  brain  control  mechanisms.  Following  the 
establishment  of  well-practiced,  stabilized  speed-of-movement  times  and 
associated  Ei‘G  patterns,  groups  of  subjects  will  be  electrically  stimulated 
while  optimizing  only  one  of  the  three  basic  control  factors  as  a  means  of 
improving  performance.  One  group  of  subjects, for  example,  will  be  stimulated 
electrically  with  only  the  agonist  to  antagonist  latency  of  activation 
manipulated  to  optimize  limb  movement.  Another  group  of  subjects  will  be 
stimulated  with  the  duration  of  EMG  activity  manipulated  to  optimize  limb 
movement.  A  third  group  of  subjects  will  be  stimulated  with  the  relative 
EMG  activity  of  agonist  and  antagonist  muscles  being  manipulated  to  optimize 
limb  movement.  These  three  groups  will  be  compared  to  a  control  group. 

Series  3.  A  third  series  of  studies  will  assess  the  relative  importance 
of  the  three  basic  neuromotor  coordination  mechanisms  taken  two  at  a  time; 
i.e.,  electrical  stimulation  will  be  applied  to  produce  a  muscle  synergy 
contraction  pattern  which  optimizes  two  of  the  three  basic  control  mechanisms. 
Such  a  paradigm  requires  three  treatment  conditions:  1  and  2,  1  and  3,  2  and 
3  with  comparison  against  a  control  group. 

Each  of  these  three  series  of  electrical  stimulation  studies  will  be  guided 
by  actual  experimental  results  and  by  results  of  the  optimization  studies  of 
mathematical  modelling.  The  exact  research  methodology  must  remain  flexible 
within  these  constraints,  especially  for  series  two  and  three.  It  is  impossible 
to  anticipate  when  all  three  series  of  studies  can  be  completed,  but  it  is 
unlikely  that  they  can  be  completed  in  one  year  and  a  request  for  an  additional 
year  of  grant  support  is  a  distinct  likelihood.  Stabilization  of  limb  movement 


time  and  EMG  patterns  requires  extensive  practice  for  one  or  more  weeks. 
Electrical  stimulation  treatment  sessions  must  be  given  for  at  least  2-3  weeks 
before  a  re-test,  and  several  such  2-3  week  sessions  must  be  given  consecutively 
in  order  to  determine  the  magnitude  of  effects. 

We  would  also  like  to  begin  to  investigate  the  potential  of  high  frequency 
electrostimulation  purported  by  the  Russians  to  produce  better  strength  gains 
than  normal  exercise  regimen.  Our  intention  is  to  utilize  the  electromuscular 
stimulator  (EMS)  now  commercially  available  to  determine  if  the  neuromotor 
coordination  mechanisms  as  expressed  in  EMG  patterns  are  affected  by  increases 
in  muscle  strength  produced  by  EMS.  Such  a  pilot  investigation  can  be 
accomplished  most  effectively  in  connection  with  the  proposed  coordination 
studies.  Our  equipment  request  reflects  the  cost  of  teh  such  Electrostim  180 
models  manufactured  by  the  NU-Med  Surgical  Supply  Co.  (S3, 250  each) . 
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PART  II 
■Vodel  1  inq 


Summary  of  Part  II: 

We  describe  a  two-compartment  model  of  the  neurcmuscul ar  system 
involved  in  the  voluntary  fast  arm  movement  to  a  target  discussed  by  Kroll 
in  Part  I.  7ne  first  compartment  accepts  averaged  biceps  and  triceps  EMG 
signals  as  inputs,  and  models  the  arm's  musculo-skaletal  response  by 
producing  elbow  angular  position  $(t)  and  velocity  $(t)  over  the  correspond¬ 
ing  movement  time.  This  is  called  cur  S/Q/ t(t)  model ,  and  its  defining 
equation  is  of  the  form 

$(t)  *  (extensor  torque  =  Qc)  -  (flexor  torque  «  Qp). 

A  detailed  version  of  this  equation  is  being  computer-simulated  with  the 
PASCAL  programme  listed  in  Appendix  A. 

Cur  second  compartment,  called  our  com/cont/E  model,  accepts  volitional 
commands  as  inputs,  and  models  the  nervous  system's  response  by  producing 
smoothed  flexor  and  extensor  EMG  signals  to  feed  into  the  S/Q/j  model.  Fig. 
11-4  gives  a  schematic  representation  of  the  com/ccnt/E  model,  which  is 
explained  in  detail  in  the  report. 

The  purpose  of  our  two  models  is  to  separate  the  various  control 
functions  in  the  ballistic  arm-movement  system  well  enough  to  understand 
where  adaptation  occurs  as  speed,  precision,  ar.d  coordination  of  movement 
improve  with  practice. 

This  is  the  first  publication  that  has  arisen  from  w  cr'<  on  t.n ;  s 


I! -A  Completed  Work 

Since  last  June  we  have  completed  an  initial  specification  of  two 
models  that  will  help  us  to  batter  understand  and  guide  the  direction  of 
Knoll's  experiments.  The  first  model  accepts  smoothed  EMG  inputs  to  flexor 
and  extensor  muscles  and  calculates  arm  angle  dynamics  for  the  fast  arm  move¬ 
ment  described  above  by  Knoll.  This  model  is  denoted  the  E/Q/aft)  medal 
(short  for  EMG/elbcw  torque  *  Q/elbow  angular  velocity  dQ/dt  =  ;) .  The 
second  model  feeds  flexor  and  extensor  smoothed  EMG  outputs  into  the  first 
model,  and  receives  its  input  commands  frem  a  volitional  fast  arm  movement 
center  in  the  brain,  presumably  via  the  pyramidal  tract  from  the  cerebellum 
to  the  brainstem  and  spinal  cord  (cf.  Miles  and  Evarts,  1979).  This  model 
is  denoted  the  Com/Con t/E  model  (short  for  command /neural  control/EMG) .  The 
appropriate  connection  of  our  two  models  provides  a  complete  overall  modal 
for  the  fast  arm  movement  as  driven  by  volitional  commands. 

3elow  we  describe  the  E/Q/?  model  first  and  the  ccm/cont/E  mocel 
second.  Section  II-2  then  proposes  further  studies  for  next  year. 


1)  T'ne  E/q/a  model 

This  model  arose  out  of  a  desire  to  formulate  a  quantitative  account 
of  the  E/Q/a  part  of  the  arm  movement  to  target  described  by  Kroll.  After 
reviewing  the  relevant  literature,  we  rejected  possible  model  formulations 
at  the  molecular  biophysical  level  as  too  complex  (cf.  Hatze,  1973,  1530, 
and  Oijkatra  et  al . ,  1973a,  1973b).  In  addition,  for  our  movement,  they 
would  be  too  difficult  to  specify  quantitatively.  We  also  rejected  model 
formulations  at  the  spring-mass  postural  level  (cf.  Sakitt,  1580)  because 
they  would  be  too  steady-state  oriented  to  allow  an  analysis  of  our  large 
transient  arm  motion.  Other  control  engineering  models  (cf.,  e.g.,  McRuer 
seemed  too  general  or  too  coarse-grained  for  our  purpose. 

Figure  II-l  gives  a  schematic  representation  of  the  framework  for  the 
E/ Q/a  model  we  chose.  Its  correspond ing  toe-level  analysis  formula  is 

*  Qc  -  Qc  ,  (i) 

w.  i 

where  eaoh  dot  ever  ?  indicates  a  time  derivative,  and  Qc  and  Qp  are  the 
respective  torques  acout  the  elbow  as  exerted  by  the  extensor  and  flexor 
groups  of  muscles  (all  the  muscles  in  each  group  operate  in  approximate 
temporal  and  tensional  unison  throughout  our  movement  (Laaasse,  1975)). 
is  the  arm's  moment  of  inertia  for  rotations  aoout  the  fixed  elbow  joint  a 
snown  in  Fig.  1 1 - 1  (over  a  a  change  from  60*  to  about  -30'). 

At  the  secorg  level  of  analytical  refinement,  we  let 

Qc  s  'J\jc  •  r  5e  (2) 

wnere  is  the  torque  due  to  viscoelastic  braking  of  the  movement 
(Lest’lenne,  197?'.  and  •»._  is  the  maximum  'over  -90*  <  C  <  90*1  iscmetr-'c 
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For  the  flexor  sice  we  let 

Qc  3  *»e  •  *, 


7F  •  'IF 


EMGg ,  recorded  over  the  last  second  (250  ms,  or  several  seconds  'would  serve 
as  well,  but  we  need  a  standard).  To  find  we  accept  the  claim  of 

Messier,  et  al . ,  1971,  that  under  static  loads,  "the  averaged  EMG  is 
directly  proportional  to  muscle  tension  and  the  constant  of  proportionality 
is  independent  of  muscle  length."  Admittedly,  the  force  properties  of 
muscles  undergoing  rapid  active  or  passive  stretch  or  compression  are  very 
different  than  when  these  muscles  are  being  actively  or  passively  subjected 
to  static  loads.  But  from  what  little  is  known,  the  best  we  can  do  to 
account  for  this  is  to  include  the  factor  appearing  in  (2).  is  a 

fraction  between  0  and  1  that  inserts  the  force-velocity  relationship 
suggested  by  Tharstensson  et  al.,  1976,  and  Perrine  and  Edgerton,  1973,  for 
the  human  knee,  and  by  Lagasse's  1975  observations  on  forearm  accelerat ions 
during  maximally  fast  movements  (quite  like  our  experimental  ones  only 
untargeted) . 

Our  Hyg  function  is 


^  *  1  "  ToCo 

m  mm  > 


where  ?  is  in  degrees/sec,  and  the  double-cusped  brackets  mean 


x  if  x  >0 


0  otherwise 


wn  are 


13  similar  to  t--  and  s,.c  =  Equating  to  j...-  for 

•  *  ( \  <  i  4 c  v  r 


movement  is  inaccurate  to  some  decree  because  the  flexor  undergoes  active 
contraction  and/or  passive  compression  wnsreas  the  extensor  urcercces 
passive  on  active  stretcn  (actually,  the  jo'nt  "ust  :e  stao’ so  fere 
are  ~o  purely  passive  muscle  states,  only  aoorcx imate  ones). 


Our  method,  then,  reduces  to  choosing  functions  in  (2)  that  best  fit 


the  averaged  data  before,  during,  and  after  enough  trials  by  each  subject 
( hundreds  1)  to  enable  him  to  achieve  maximum  speed.  The  resulting  functions 
should  tell  us  some  valuable  quantitative  things  about  the  muscle  properties 
involved. 

At  the  last  level  of  analytical  ref inement,  we  let 


t-u 


♦IF(t)  *  ApCOS*(t) 


kF‘Era?p(t)hF(t-T)  dt. 


t- 


(5) 


where  ApCOSi(t)  converts  flexor  force  to  torque  according  to  the  geometry  of 
Fig.  II-l;  tNiflp ( t )  is  the  smoothed  EMGp  envelope  at  time  t;  kc  is  a  constant 
of  proportional i ty  that  we  assume  converts  EMG-(t)  into  flexor  force  at  t 
(Messier,  at  al . ,  1371);  and  hp(t-?)  is  the  normalized  flexor  twitch 
response  at  time  t  to  a  unit  kpS-MC-.-  impulse  that  arrived  at  time  t. 

Our  hp(t-t)  is  defined  in  Fig.  1 1-2,  where  we  equate  hp(t-t)  with  its 
extensor  counterpart  (for  want  of  knowledge  to  do  otherwise.  Tnis  could  be 
incorrect  by  perhaps  10%  or  so,  especially  after  many  practice  trials).  The 
shape  of  our  h  function  is  a  s-truncaticn  (for  economy  in  computing)  of 
Qijkstra  et  al.'s  1973a)  twitch  response  function  for  biceps  orach ii  muscle. 
3ut  whereas  their  3^  =  .015  and  3^  =  .15,  we  began  with  3^  3  .05  arc  3~  *  .3 
As  we  continue  to  fit  the  incoming  data,  our  3 ‘ s  will  doubtless  change  in 
the  direction  of  Oijkstra  et  al.'s. 


Equation  (5)  assumes  1  inear' ty 
This  is  probab’y  scmewnat  inccr-ecc 
makes  the  sire  ’i-aarity  assumption 


in  t.me  accumulation  of  twitch  -esccnses. 
for  cur  mcveme-t.  veo  Sak;tt,  13S0, 
for  an  arm  mover  am:  similar  to  Ouns  a  no 
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obtains  good  results.  In  the  future,  however,  we  plan  to  try  letting  h(t-t) 
vary  as  a  function  of 

i 

t-T 

1  _  -a( t-T-a) 

I  fcMG(  <i)e  do,  (6) 

t-T-V 

l 

the  exponentially-weighted  S.MG  input  from  t-T  back  to  t-x-v  in  time.  This 
seems  like  the  most  expeditious  way  to  introduce  an  appropriate  nonlinearity 
f  for  the  accumulation  of  twitch  responses  into  equation  (5). 

Letting  Apkp  a  Cp  in  (5)  we  get 


♦jp(t)  a  CF  cos  $( t ) 


t-u 

% 

WF(T)hF(t-t)dT 
t -»-« 


(7) 


> 


Defining  ^  similarly  and  neglecting  in  (2)  (as  does  Sakitt,  IS81.  cf. 

also  Lestie.nne,  1979),  our  final  ref ineo  version  of  equation  fji  becomes 


5(t) 


3 


•llyCOU(t) 


u 


t-u 

* 

c^Mh^t-Odt  - 

Si 

t-u-5 


t-u 


I?!Gc(T)h-(t-T)aT 


(8) 


with  the  constraining  assumptions  g;/en  ao  eve.  To  a  first  accrcx  ir.af'cr. 

0)  ty  ren’JCi-,c  Jj  oy  the  isometric  tor  cue  founc 


first  in  the  flexor  end  then  in  the  extensor  directions  for  each  subject. 
There  are  various  good  ways  of  measuring  each  subject's  J  (e.g.,  Hatze, 

1980,  refers  to  one  he  has  developed),  but  we  have  not  yet  selected  one. 

We  now  discuss  how  the  S.MG  signals  in  (3)  are  obtained.  Figure  1 1-3 
shows  an  example  test  record.  The  E.MG.„  signals  are  rectified  and  filtered 
raw  EMG  signals.  EMG_W  signals  are  derived  in  an  on-line  computer  which 
"continuously"  leaks  off  a  fraction  of  the  present  integral  of  the  rectified 
EMG  signal.  Oenoting  the  latter  EMGr,  EMG  in  Figure  1 1 -3 ,  to  a  good 
approximation,  is  given  by 


EMG(t] 


w 

J  £MGr(t) 


•a(t-t), 


Thus  if  EMG  (t)  3  K,  a  constant,  then 


M(t)  «£[1  -  e'at] 


as  t  gets  large.  Also,  if  EMGr( t)  »  Ct,  C  a  constant,  then 

EMG(t)  *  -t  -  £,[1  -  e‘at]  -  -  t  (11) 

a  a  a 

a 

for  small  t.  Finally,  it  is  easily  shown  that  if  £MG  3  K  as  in  (10)  for  a 
long  time  and  then  EMGr  becomes  0  at  t, 


d  £MG(  t-t-e) 

dc 


-  ;<e"ac  -  -X 


for  small  c.  Comparing  (10)  and  (12)  gives  us  a  way  to  determine  3  from 
apcrccr ’ ate  :a:a  samples. 

Ecuaticrs  110'  -  1 12)  enable  us  to  find  easily  the  i.'-'Gi'tl  ouant't’er 


'•scores  such  is  shown  in 
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EMGC  signal  in  Fig.  1 1 -3 ,  we  form  a  straight-line-segmental  approximation  to 
it  by  connecting  a  to  b,  b  to  c,  etc.  If  in  this  approximation,  a  segment 
is  essentially  constant  at  K,  we  let  EMG  there  equal  Ka.  If  the  segment 
rises  with  slope  C  >  0,  we  let  cMG  rise  correspond ingly  with  slope  Ca.  If 
the  segment  falls  from  K  with  slope  -K,  we  let  SMG  there  equal  0.  Thus  our 
£MG( t )  equals  a  times  the  straight-line-segmental  approximation  to  EMG 
except  where  that  approximation  falls  off.  The  falloff  intervals  are  easy 
to  handle  because  the  EMG  signals  for  cur  movement  are  either  on  at  high 
intensity  or  nearly  off  over  most  of  every  trial. 

There  is,  however,  one  adjustment  to  the  tMG(t)  signal  derived  above 
that  we  need  to  make.  In  order  to  compensate  for  cMG  inaccuracies  during 
EMG  falloff  phases  and  intrinsic  EMG  distortions  of  nervous  excitation  to 
muscle  at  lew  signal  levels,  we  shall  augment  our  EMG  signal  with  a  *(t) 
postural  feedback  term.  This  is  justified  by  assumption  4  in  the  next 
section. 

To  make  the  adjustment,  we  denote  the  previous  IMS’  by  IMG*,  and  add  a 
feedback  term  denoted  EMG*  .  This  gives  us 

cMG(t)  alMG*(t)  +lMG**(t)  (13) 


)]  G(t 


(1*1) 


on  both  the  flexor  and  extensor  sides.  We  define: 

cS5p*(t)  3  PF1  ^(t-20)  +  PF2i(t-20)  +  ?F3?( t-20) 

where  the  double-cusped  brackets  mean  the  same  as  in  (4);  where  tne  quantity 
inside  these  brackets  is  j(t)  as  predicted  frem  >,  >,  and  5>  20  ms  earlier 
(20  ms  being  the  total  delay  around  the  loco  from  y  receptcr  to  a  stimulation 
of  muscle);  wnere  G(t)  *  0  until  t  first  reaches  1  degrees,  and  G(t)  *  1  there¬ 


after;  and  wnere  the  P-,  are  constants  (new  set  at  =  .25,  =  1C,  *  ’ 

r  l  ri  r 2  r ~ 

•  ##  j 

—  recall  that  j  is  in  degrees/ms  anc  V  is  in  aegrees/ms") .  Similarly, 
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where  P£i  *  Ppr 

Finally,  to  keep  (14)  and  (15)  from  causing  everlasting  $  oscillations, 
we  change  (1)  to 

J*(t)  *  Q£ ( t )  -  Qp(t)  -  Cci(t)  (16) 

where  the  C  viscoelastic  damping  factor  (cf.  Lestienne,  1979)  is  set  "near" 
the  level  that  just  prevents  a  second  overshoot  of  the  0  degree  target 
position. 

We  now  consider  adaptation. 

As  our  subjects  accumulate  more  experience,  we  know  from  past  results 
(Lagasse,  1975,  e.g.)  that  their  movement  times  will  decrease.  If  the 
adaptations  enabling  this  occur  partly  in  the  muscles,  this  should  require 
us  to  modify  the  h  functions  in  (8).  Since  in  previous  pilot  experiments 
our  flexor  movement  adapted  several  times  faster  and  further  than  the 
opposite  extensor  movement,  it  could  happen  that  only  h£  in  (8)  adapts  with 
practice.  It  also  might  develop  that  all  adaptation  occurs  in  the  nervous 
system  (cf.  the  ccm/ccnt/E  model,  next  Section). 

The  crucial  feature  of  our  E/O/a  model  is  that  eventually  it  should 
locate  and  characterize  the  adaptive  nonlinear  effects  involved  when  speed 
of  movement  is  increased  with  practice.  This  should  give  us  a  mechanical 
perspective  on  how  to  train  for  soeed.  precision,  grace  and  skill  —  in 
short,  for  well  coordinated  fitness. 

We  have  included  as  Appendix  A  a  listing  of  the  PASCAL  programme 

written  by  Valerie  Congdon  to  ccmpute  *(t)  and  »(t),  as  outlined  above, 

using  a  file  containing  the  SMG„( t)  and  EMG- ( t )  data  for  one  movement.  The 

r  c 

program  is  well  documented,  hierarchically  organized  into  procecural  ana 
functional  levels  icentical  to  the  analysis  levels  in  this  rapcrc,  arc 
specified  so  as  to  be  easy  to  orange.  Ms.  Congdon,  who  will  *eave  L.Mass 
in  January,  1951,  has  cone  an  exemplary  job. 
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2.  The  Com/Cont/£  Model 

This  model  concerns  the  functional  organization  of  the  neural  control 
system  for  the  E/Q/?  model  described  above.  It  is  an  attempt  to  make 
precise  only  those  conclusions  that  have  been  formulated  by  several  experi¬ 
menters  as  necessary  consequences  of  their  tests  on  fast  arm  movements 
involving  the  elbow  and  sometimes  one  other  joint.  The  com/cont/E  model  is 
amenable  to  few  significant  neuroanatcmical  interpretations.  Eventually  we 
hope  to  obtain  a  good  structural  model  that  does  not  suffer  from  this 
defect,  but  that  may  be  a  couple  of  years  away. 

Our  specification  of  the  com/cont/E  model  arose  out  of  the  following 
Assumptions 

1.  For  about  the  first  100  ms  at  least,  our  fast  arm  movement  is  entirely 
feedforward  (is  independent  of  any  feedback),  and  thus  is  unaffected  by 
stretch  reflex  or  tendon  affarents  (cf.  W adman,  et  al . ,  1979;  Desmedt 
and  Godaux,  1973) . 

2.  Selecting  the  levels  and  l,  of  flexor  and  extensor  EMG  activity 

re 

require  the  first  volitional  decisions.  (As  noted  by  Krcll  above,  all 
of  the  flexor  muscles  act  aoproximataly  in  temporal  and  tensional 
unison  with  the  biceps  during  our  arm  movement.  Likewise  for  the 
triceps  and  the  extensor  muscles.)  These  levels  are  selected  to  sat 
the  speed  and  amplitude,  M,  of  the  movement  with  allowances  for  expected 
inertial  loading  of  the  limb  and  estimated  fatigue  of  the  muscles  being 
activated  (cf.  Lestienne,  at  al.,  1979;  Lestienne,  1979;  Wadman,  et 
al . ,  1979;  Wadman,  et  al . ,  198G). 

3.  The  next  movement  decision  is  to  set  the  relative  timing  cf  the  first 


2‘ 


tion  of  magnitude  and  speed,  again  allowing  for  the  expected  loading 
and  muscle  fatigue  (cf.  references  in  2  above).  Force-velocity  and 
force-length  effects  are  probably  also  included  in  this  reckoning 
(Thorstensson,  et  al.,  1976;  Perrine  and  Edgerton,  1978),  as  well  as 
viscoelastic  braking  (Lestienne,  1979;  Maughan  and  Godt,  1979). 

4.  The  flexor  and  extensor  muscle  tension  levels  at  the  end  of  the  movement 
epoch  when  the  postural  maintenance  phase  is  entered  are  variable,  but 
the  ratio  R  between  these  two  levels,  (as  reflected  by  their  leaky- 
integrated  EMG  levels)  is  a  fixed  postural  command,  and  is  not  affected 
by  the  starting  point,  direction,  amplitude,  or  speed  of  movement 
(Lestienne,  et  al.,  1979;  Sakitt,  1980).  Presumably  the  EMG  levels 
themselves  at  the  end  of  the  movement  epoch  are  quickly  adjusted  if  need 
be  (how  fast?)  to  obtain  the  desired  postural  stretch  reflex  stiffness 
(Houk,  1976).  Since  the  biceps  tendon  reflex  loop  delay  is  about  35 
ms,  the  stretch  reflex  loop  delay  is  about  50  ms  (Marsden,  et  al . ,  1976) 
and  our  actual  am  movement  duration  is  about  300  ms,  the  presently 

.  unknown  speed  with  which  the  sat  point  of  the  stretch  reflex  is  reset 
and  the  reflex  subsequently  engaged  will  be  imortant  to  discover. 

5.  In  males  at  least,  continued  practice  of  a  fast  am  movement  reduces 
execution  time  and  alters  the  movement's  EMG  signature  (Lagassa,  1975). 


In  the  light  of  the  above  assumptions,  we  now  cive  a  point-by-coint 


description  of  the  ccm/cont/E  model  (Figure  II-A); 

1.  At  the  top  of  Figure  1 1 -4,  R,  L,  M,  ILF,  and  FF  are  commanas 
from  the  volitional  decision  center  for  our  *ast  arm  movement 
"center"  might  be  distributed  over  a  large  part  of  the  brain) 
ratio  described  in  assumption  -1  above;  L  times  gc,  L.c,-,  and 


controller  of  volitional  fast  arm  movement.  See  text 
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levels  Lp  and  Lp  described  in  2  above;  M  is  defined  in  2;  ILF,  the 
inertial  loading  factor,  and  FF,  the  fatigue  factor,  referred  to  in  2 
and  3  are  appropriately  scaled  for  multiplication  into  the  subject's 
moment  of  inertia,  J,  for  the  arm  movement  in  question.  Thus  ILF  _>  1, 
with  IFL  *  1  when  no  external  load  is  added;  and  0  <  FF  £  1,  with  FF  *  0 
interpreted  as  complete  paralysis  and  FF  =  1  as  zero  fatigue  (defined  in 
terms  of  force  development  capability). 

2.  Np  and  Np  are  loosely  interpretabl e  respectively  as  flexor  and  extensor 
motorneuronal  pools  that  drive  either  the  flexor  and  extensor  arm 
musculature  for  our  movement  or  just  the  biceps  and  triceps. 

3.  At  t  =  0  the  flexor  EMG  is  activated,  t*  is  the  time  at  the  end  of  the 
first  extensor  EMG  burst.  The  time  courses  of  R  and  L  with  respect  to 
t*  are  hypothesized  as  shown  at  the  top  of  the  figure. 

4.  At  t  *  0,  L  turns  NR  on,  which  then  turns  Np  on  (see  5  below),  all  with 
delay  less  than  de.  Nf's  output  firing  rate  rises  and  falls  somewhat 
sporadically,  but  peaks  at  the  L.gp  level.  The  stretch  reflex  pathway, 
SPp,  into  Np  is  intially  inactive  as  noted  in  assumption  I  above,  "p's 
output  equals  the  algebraic  sum  of  its  inputs.  Np  keeps  Np  shut  off 
over  the  reciprocal  inhibitory  pathway. 

5.  At  time  D,  a  shot  out  of  Ng  turns  Np  off  to  end  tne  first  flexor  EMG 

burst,  c  and  k  in  the  formula  for  D  are  scaling  constants.  Since  L 

turns  Nn  on,  the  value  of  c  would  be  0  if  the  flexor  muscle  M„  were 
U  r 

perfectly  linear  in  L  and  if  no  co-contraction  of  M_  and  Mr  were 

”  w 

necessary  to  stabilize  the  elbow  joint  during  movement.  S' nee  neither 
of  these  conditions  holds,  c  >  0.  After  N_  shuts  off,  N-  is  no  lencer 

r  - 

innibitac  arc  turns  on.  Tnen  Nr  keeps  N_  s.nut  of*  by  rec’procal 
inhibition. 
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6.  Backtracking  a  bit,  only  one  of  N^'s  outputs  is  on  at  a  time.  Initially 
l  turns  N^'s  ip  output  on,  which  is  necessary  to  keep  N-  on.  Later 
Ng's  output  shot  switches  N^'s  output  over  to  i which  is  necessary  to 
keep  on. 

7.  N^'s  input  Tp  provides  NR  with  an  estimate  of  Mp's  tension,  as  computed 

(or  simulated)  by  Nh-.  The  function  h-  is  Nh-'s  estimate  of  M-'s  twitch 

r  r  r  r 

response  function  hp  (see  Part  II-A-1)  on  tne  E/Q/*  model).  In  contrast 
with  the  nervous  system,  we  would  model  Nhp’s  output  at  time  t  with 
t 

A 

N_(o)hp(t-a)da.  3y  assumption  4  above,  the  flexor  tendon  organ 

o 

response,  TORp,  is  affectively  unavailable  to  Nhp  for  at  least  30-85  ms. 
Mutatis  mutandis,  NR's  extensor  circuits  are  the  same  as  its  flexor  ones 
3.  If  our  movement  is  entirely  feedforward,  some  such  centers  as  Nhp  and 
Nhr  are  necessary  to  control  Np  and  N-'s  development  of  an  appropriately 
stiff  arm  posture  at  the  end  of  the  movement  (assumption  4  above).  The 
possibility  that  there  are  Nh-  and  Nh-  circuits  in  the  cerebellum  seems 

r  C 

to  fit  Eccles'  concept  of  the  cerebellum  as  a  control  reference  in  fast 
voluntary  movements  (Eccles,  1977).  Lagasse,  1975 ’ s  results  showing 
that  EMG  temporal  organizations  for  fast  human  limb  movements  don't 
change  with  fatigue  support  the  idea  that  such  movements  are  entirely 
feedforward . 

A 

9.  When  N^'s  T-  input  reaches  the  desired  final  level,  N^'s  output  switches 

A 

from  i-  to  i-  in  case  T_  is  not  up  to  the  desired  final  level,  and  to 

c  i  r 

A 

in  case  T  is  up  to  its  final  level.  An  active  z  sicnals  the  end  of 
c  c 

the  movement  arc  suoseouent  'eirstatement  of  tne  normal  pcstural  mainte- 


nance  mechanisms.  The  desired  final  flexor  and  extensor  tensions  are 


* 


denoted  7p  and  7^,  and  are  assumed  to  satisfy  7-3  *  7p  =  l  (assumption  4 
above).  Thus  (Tp/T^)  2  R,  and  the  NR  component  in  Fig.  1 1-4  has  suffi¬ 
cient  inputs  to  compute  all  of  its  outputs. 

10.  If  we  assume  some  background  $?p,  SP^,  and  other  noise,  the  sequence  of 
events  in  9  will  not  always  be  as  simple  as  described.  The  first 

Np  or  Np  E.MG  burst  might  sometimes  be  cut  off  prematurely,  forcing 
another  one  or  even  two  small  bursts  later  on  (see  the  l  time  course, 
top  right  in  Figure  1 1-4) .  The  famous  third  EMG  burst  (e.g.,  Wadman, 
et  al.,  1980)  might  be  accounted  for  in  this  way  —  an  attractive 
hypothesis  because  Kroll's  experiments  to  date  with  highly  practiced 
subjects  usually  don't  show  a  third  EMG  burst.  The  ability  of  the 
control  components  in  Fig  I 1-4  to  learn  (with  lots  of  practice)  the 
amount  of  viscoelastic  braking  that  actually  occurs,  the  force-velocity 
and  force-length  relationships  of  the  muscles,  and  various  nonlineari¬ 
ties  of  the  system,  offer  other  adaptation  possibilities. 

11.  The*->symbol  in  Fig.  1 1-4  points  to  likely  adaptation  sites  as  the 
same  movement  is  practiced  over  successive  trials  and  days. 

12.  and  the  (M^,  Mhp,Nhff)  complex  could  be  ccmbined  instead  of  function¬ 
ally  partitioned  as  shewn,  but  that  would  defeat  the  purpose  of  the 
com/cont/E  model . 


im 


1 1-8  Proposed  Modelling  Studies  Next  Year 
We  wish  to  proceed  on  three  fronts: 

1.  The  E/Q /*  model  will  be  adjusted,  and  perhaps  revised  or  augmented,  to 
fit  Kroll's  experimental  results.  His  tests  with  different  loads  and 
fatigue  levels  should  be  very  helpful  to  us  during  this  phase. 

2.  The  Com/Cont/E  model  will  be  scrutinized.  In  particular,  the  and  N.n 
components  will  be  checked  for  plausibility  against  Kroll's  fatigue 

and  load  variational  test  data.  We  now  believe  that  ''efferent-copy1' 

M„  simulators  such  as  the  functions  actually  exist  (perhaps  at  the 
highly  integrative  thalamic  level  atop  the  brain  stem  reticular  forma¬ 
tion,  or  in  the  cerebellum),  but  we  shall  seek  more  evidence  to  support 
this  idea.  We  also  plan  pilot  tests  of  various  kinds  to  further 
elucidate  the  role  of  stretch  reflexes  at  different  times  during  and 
immediately  following  fast  arm  movements.  To  date,  Wadman  et  al. 

(1979)  have  shown  that  when  a  subject's  planned  fast  arm  movement  is 
unexpectedly  blocked,  the  pattern  of  EMG  activity  over  at  least  the 
first  100  ms  is  the  same  as  before  blocking.  But  we  know  of  nothing 
else  in  the  literature  on  the  dynamics  of  rapidly  changing  the  set 
points  of  stretch  reflexes. 

3.  We  shall  develop  computer  programs  in  PASCAL  for  solving  the  nonlinear 
delay-differential  equations  that  will  arise  wnen  the  E/Q/t>  model  is 
augmented  to  include  feedback  stretch  reflex  information.  Banks 
(1979,  1980)  has  recently  oevelcped  the  requisite  mathematical  methods, 
but  as  yet  PASCAL  programs  to  implement  them  do  not  exist. 

In  1  to  3  above,  we  shall  oe  especially  watchful  of  adaptation  parame¬ 
ters  and  character istics .  In  year  three  of  the  contract,  we  snail  apply 
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An  EMG-Level  Muscle  Model  for  a  Fast  Arm  Movement  to  Target* 

W.  Kilmer,  W.  Kroll,  V.  Congdon 

Departments  of  Electrical  and  Computer  Engineering, 

Exercise  Science  and  Computer  and  Information  Science 
University  of  Massachusetts 
Amherst,  Massachusetts  01003  USA 

ABSTRACT 

A  model  of  human  muscle  action  is  presented  for  a  maximally  fast,  large- 
amDlitude  forearm  movement  to  target.  The  inputs  to  the  model  are  approxi¬ 
mately  the  biceps  and  triceps  EMG  envelopes  over  a  single  movement.  The 
model’s  output  gives  the  corresponding  displacement  angle  of  the  forearm 
about  a  fixed  elbow  position  as  a  function  of  time.  The  idea  of  the  model 
is  to  conceive  of  both  EMG  input  drives  as  successions  of  millisecond  input 
pulses,  with  each  pulse  resulting  in  a  muscle  tension  twitch.  Every  twitch 
is  amplitude-scaled,  parametrically-shaped,  and  duration-limited  as  a  func¬ 
tion  of  the  muscle's  contractile  history  thus  far  in  the  movement.  The 
muscle  tension  at  any  time  t  is  the  sum  of  the  residual  tension  levels  of 
all  twitches  begun  before  t. 

The  model  was  develooed  and  tested  with  special  reference  to  two  subjects 
one,  according  to  the  model  dynamics,  was  a  comparatively  slow-twitch  type, 
and  the  other  modelled  as  a  fast-twitch  type.  Good  agreement  was  found 
between  model  output  and  subject  response  data  whenever  the  subject's  EMG's 
were  "synchronous."  The  model  can  be  used  to  characterize  each  subject's 
responses  by  a  suite  of  twitch  characteristics.  This  will  enable  us  to 
check  the  accepted  but  now  suspect  correlation  between  muscle  biopsy-  and 
performance-determined  muscle  twitch  tyne. 

*This  study  was  supported  by  contract  DAMD1 7-3C-C-01 01  from  the  U.S.  Army 
Medical  Research  and  Development  Command.  The  views,  opinions  and/or  find¬ 
ings  contained  in  this  report  are  those  of  the  authors  and  should  not  be 
construed  as  an  official  Department  of  the  Army  position,  policy,  or 
decision,  unless  so  designated  by  other  documentation. 


1 .  INTRODUCTION 


We  present  a  mathematical  model  of  the  elbow  flexor  (biceps  group) 
and  extensor  (triceps)  muscle  actions  for  the  human  arm  movement  outlined 
in  Fig.  1.  The  inputs  to  the  model  are  the  EMG  signals  referred  to  in 
the  caption  of  Fig.  2.  They  approximate  the  envelopes  of  raw  EMG  signals, 
are  are  derived  from  EMG  signals  such  as  shown  in  Fig.  2a)  and  b)  by  the 
method  described  in  Kilmer,  et.  al .  (1981).  The  output  of  the  model  cor¬ 
responding  to  one  trial  of  the  type  in  Fig.  2  is  a  calculation  of  an  arm 
displacement  function  such  as  shown  in  Fig.  2c). 

This  is  the  first  attempt  we  know  of  to  analyze  EMG(t)-to-TT(t)  trans¬ 
formations  whereT  is  the  tension  noted  in  Fig.  1)  in  muscles  controlling 
fast,  large-amplitude  limb  movements  to  a  target.  The  most  closely  related 
modelling  works  to  date  are:  Stein  and  Oguztoreli  (1981) 's  small  signal 
linear  analysis  of  the  mamnal ian-neuromuscular  and  reflex  system;  Hatze's 
(1978,  1980,  1981)  papers  at  a  more  detailed  level  involving  state  equations 
for  numbers  of  motor  units  in  different  contractile  states;  Sakitt's  (1980) 
steady-state  spring  model  for  postural  changes  similar  to  that  shown  in 
Fig.  1;  and  Dijkstra  et.  al.'s  (1973a,  1973b)  original  cross-bridge  model 
and  analog  computer  simulation  of  small  amplitude  human  forearm  movements. 

Our  model  attempts  to  better  determine  the  correlations  between:  1) 
the  SDeed  (power)  and  shape  of  maximally  fast,  large-amplitude  human  arm 
movements  to  a  target;  2)  the  proportions  of  fast-  vs.  s low-twitch  fibers, 
as  determined  by  muscle  biopsy  in  subjects  making  such  movements;  and  3) 
the  physiologically  effective  twitch  speeds  involved  at  different  stages 
of  each  such  movement  when  performed  by  highly  skilled  (i.e.,  well-practiced) 
subjects.  As  discussed  below,  some  evidence  has  accumulated  to  suggest  that 


the  fast-  vs.  slow-twitch  dogma  of  mammalian  skeletal  muscle  action  may 
be  inaccurate  to  some  degree. 

Several  investigators  have  experimented  with  human  movements  that 
differ  from  that  of  Fig.  1  primarily  in  that  they  are  performed  horizont¬ 
ally  so  as  to  eliminate  the  chanqina  effects  of  gravity.  Most  notable 
among  them  are  Wadman,  et.  al .  (197S,  1980a,  1980b),  Lestienne  et.  al . 
(1979),  and  Lestienne  (1979).  Lacasse  (1975),  and  Kroll  and  associates 
(e.g.  Kroll  (1974)),  however,  have  collected  the  most  suitable  data  for 
our  purposes  by  using  essentially  the  Fig.  1  paradigm  on  subjects  whose 
EMG  signatures  were  always  cleaned,  sharpened,  and  stabilized  through 
extensive  pre-test  practice.  Vie  chose  the  resulting  so-called  "synchronous" 
EMG's  to  analyze  with  our  model.  (We  note  below  that  the  confounding  effect 
of  aravity  in  rig.  1  is  apparently  negligible.) 

Kroll's  protocol  for  each  of  his  29  subjects  (about  half  men  and  half 
women,  all  riqht-handed)  involved  four  practice  days  with  50  trials  per  day, 
followed  by  five  exercise  (i.e.  test)  days  with  20  trials  to  establish  base¬ 
line  and  5  test  trials— all  for  each  of  three  load  levels.  We  derived  all 
our  model  inputs  from  trials  using  the  smallest  added  load,  namely  three 
times  the  arm's  moment  of  inertia  for  the  movement.  This  gave  us  signal 
plots  of  about  the  right  size  on  all  channels  for  the  modelling  accuracy 
we  sought.  All  subjects  were  told  to  "make  the  movement  as  fast  as  you  can, 
and  stOD  at  the  0°  target  position."  The  emphasis  was  on  speed  of  movement, 
not  targeting  precision.  Under  these  circumstances,  the  published  data  on 
more  or  less  fast  movements  of  other  human  body  parts  (eyes,  head,  hands, 
fingers,  legs)  are  often  at  variance  with  our  results. 


We  selected  trial  data  from  two  subjects,  Hintr  and  Chan,  both  25-year 
old  males,  to  help  explain  our  models.  In  the  next  Section,  we  present  the 
conceptual  organization  of  our  basic  model,  the  U-Model,  pictorially,  leav¬ 
ing  almost  all  of  the  mathematical  formulae  till  later.  Then,  following  a 
mathematical  specification  of  U-Model,  we  show  how  it  is  augmented  to  J-Model 
to  fit  the  data  of  Chan's  much  faster  trials.  We  conclude  with  some  possible 
applications  of  J-Model . 

2.  CONCEPTUAL  MODEL 

The  purpose  of  our  models  is  to  better  understand,  in  a  quantitative 
way,  how  nervous  drives  to  muscles  as  represented  by  their  EMG's  lead  to 
schedules  of  muscle  tension  and  thence  to  movements  of  the  kind  outlined  in 
Fig.  1.  The  U-Model,  and  its  J-Model  augmentation,  are  embodied  in  a  PASCAL 
program,  called  FASMOV,  that  accepts  piecewise-1 inear  approximations  of 
biceps  and  triceps  EMG  waves  taken  from  a  dynograph  recorder,  and  produces 
a  prediction  of  the  subject's  corresponding  arm  movement  as  a  function  of 
time.  Kilmer,  et.  al .  (1981)  specifies  how  FASMOV  computes  effective  nervous 
drives,  EMGp  and  EMGp,  from  flexor  and  extensor  dynograph  EMG  outputs,  EMGp 
and  EMGp,  'respectively. 

Fig.  3  gives  a  graphic  outline  of  the  main  concepts  employed  in  the 
U-Model.  At  the  top  of  the  Figure  we  imagine  a  constant  flexor  EMGp(t) 
drive  that  lasts  for  100  ms.  We  conceive  of  this  drive  as  a  succession  of 
100  contiguous  1  ms-pulses,  each  of  amplitude  EMGp  (thus  there  should  be 
about  ten  times  as  many  vertical  lines  as  shown  in  EMGp(t)) .  We  assume 
that  each  millisecond  EMGp  pulse  produces  an  ensuing  twitch  of  tension 
whose  time  course  is  like  that  of  a  small  wave  in  the  baseline  train  of 
Fig.  3b),  and  whose  area  is  a  scale  factor  times  the  EMGp  value  (Meissier, 
et  al.,  1971).  The  shape  of  the  leading  twitch  in  the  Fig.  3b)  train  is 
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essentially  that  used  by  Dijkstra  et  al.  (1973a).  We  denote  the  electro-, 
mechanical  delay  between  the  time  of  an  EMGp  pulse  and  when  the  correspond¬ 
ing  twitch  starts  by  yp.  A  value  of  60  ms  for  yp  has  given  us  a  best  fit  to 
the  data  in  all  our  models  so  far.  This  value  is  midway  between  the  73  ms 
and  76  ms  values  found  experimentally  for  biceps  by  Lagasse  (1975)  and  by 
Ward  (1978)  respectively,  and  the  55  ms  and  52  ms  values  found  experimentally 
for  triceps  by  Lagasse  (1975)  and  by  Morris  and  Beaudet  (1980)  respectively. 

In  Fig.  3b)  our  twitch  durations  range  from  60  ms  at  the  outset  to  90  ms 
at  the  end.  These  durations  are  generally  compatible  with  those  measured, 
e.g.,  in  human  thumb  flexor  muscles  by  Gydikov,  et  al.,  1976.  Their  fast 
fatiguable  motor  unit  twitch  times  ranged  frcm  91  to  125  ms  before  fatigue  and 
from  125-333  ms  after.  We  do  not  suppose  -that  our  90  ms  terminal  twitch  time 
is  a  "fatigued"  duration;  rather  we  refer  to  it  as  a  "saturated"  duration. 

Note  that  we  clipped  our  saturated  twitch  durations  off  prematurely  as  shown 
in  Fig.  3b);  this  was  only  to  save  computer  memory  during  simulation  runs. 

Letting  £p  denote  flexor  twitch  duration,  we  assume  that  a  millisecond 
EMGp  pulse  at  time  t  causes  a  tension  twitch  given  by 

'  (sxp ( - B^p ( t -Up ) )  _  psxp( -B2p(T-Up) ) ) 

hf(-)  -  ■  •  for  t  +  uF  s  i  it  .  »f  *  Ef 

.  0  for  alt  other  t  (2.1) 

where  ?  =  0  until  the  first  adjustment  of  twitch  shape  in  the  Fig.  3b)  train, 
and  p  =  1  thereafter. 

Except  for  the  p  factor,  (2.1)  is  Dijkstra's  (1973a)  formula.  Our  2^ 
equals  their  ^  =  0.15,  but  whereas  their  S-|  =  0.015,  we  allowed  our  £^p  to 
range  frcm  0.05  initially  to  0.01  at  saturation  in  both  our  models.  For  every 
different  3-jp  value  we  recalculated  qp  in  (2.1)  so  as  to  renormalize  the  area 
under  hp(-:)  to  nearly  1.0  (details  in  Kilmer  etal.  (1981)). 
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Fig.  3b)  depicts  a  gradual  change  in  the  h»(*)  twitch  shape,  which  was 
our  original  idea.  But  to  help  cut  FASMOV' s  run  costs  we  programmed  FASMOV 
to  change  hp ( • )  shapes  less  gradually  (i.e.,  less  often).  3ecause  qp  (among 
other  things)  has  to  be  changed  each  time  hp(*)  changes,  we  refer  to  the 
entire  set  of  hp(*)-  or  h^( *)-changing  operations  in  FASMOV  as  a  "reque". 

Fig.  4  shows  all  but  the  last  hp(*)  and  hp(-)  that  were  used  during  the  Hintrl 
run,  where  only  nine  reques  in  all  occurred  (the  details  on  when  and  how  to 
reque  appear  in  Kilmer  et  al .  (1931)). 

The  total  flexor  tension  waveTpU)  in  Fig.  3b)  is  found  by  adding  up  at 
each  time  all  of  the  hp(t)  wave  values  that  are  not  zero  at  t.  These  nonzero 
values  would  be  the  hp(-r)  values  caused  by  EMGp  pulses  occurring  between  the 
times  (t  -  ( Up  +  Cp) )  and  (t  -  up).  £p  here  is  a  twitch  duration  such  that 
the  oldest  nonzero  twitch  at  t  rose  from  zero  at  time  t  -  Cp  and  Cp  was  the 
twitch  duration  in  effect  then  (cf.  next  Section).  Our  method  of  producing 
theT^(t)  wave  assumes  linearity  in  the  accumulation  cf  twitch  effects.  This 
is  probably  somewhat  incorrect,  but  Sakitt,  1980,  e.g.,  made  the  same  assumption 
for  an  arm  movement  similar  to  ours,  and  obtained  good  results.  The  simplifica¬ 
tion  was  justified  when  the  model  outputs  agreed  very  well  with  the  subject's 
data. 

On  the  right  end  of  Fig.  3b)  we  show  the<7-p(t)  wave  jumping  up  from  its 
sum-of-h(r) ' s  value,  given  by  the  dotted  curve,  to  an  amplitude  about  three 
times  greater.  Let  S^  denote  the  multiplication  factor  in  this  jumD,  which 
occurs  when  i  changes  from  negative  to  positive,  i.e.  when  the  biceps  stops 
contracting  and  starts  stretching.  Heckathcrne  and  Childress  (1981)  show  curves 
for  slowly  lengthening  and  contracting  rates  in  human  biceps  that  are  strik¬ 
ingly  compatible  with  the  values  for  which  our  models  best  fit  the  subjects' 

data:  =  3.0  in  all  Hintr  runs,  and  =  2.5  in  all  Chan  runs.  We  use  the 
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same  SR  value  for  both  triceps  and  biceps  in  each  model  simulation  run;  but 
since  stretching  one  of  these  muscles  means  contracting  the  other,  a  jump 
of^Tpit)  to  SRTF(t)  is  always  accompanied  by  a  jump  of*Tp  (t)  to7^(t)/SR,  and 
vice  versa. 

Fig.  3c)  outlines  our  last  major  concept.  As  the  integral  of°Tp(t)  grows, 
the  muscle  becomes  less  efficient  in  converting  EMG(t)  drive  into  tension.  We 
represent  this  growing  inefficiency  by  a  multiplying  factor  forT^(t),  which 
we  denote  by  exp(-cj£)  *exp(-np) —  the  rationale  appears  below.  The  continuous 
curve  in  Fig.  3c)  represents  a  typical  exp(-a^)  during  a  movement.  The  quantity 
-ap( t)  is  just  a  scaled  leaky  integral  of<Tp('0  from  x=0  to  t=t.  The  dotted 
adjustment  to  exp(-c*)  in  Fig.  3c)  is  due  to  the  exp(-np)  factor.  Initially 
?,p  is  nonzero  (ranging  from  -0.5  to  -1.8  in  simulation  runs  so  far);  this 
quantity  represents  a  "primed"  tension-producing  state  of  the  muscle  at  tension 
onset.  At  the  first  reque  of  each  simulation  run,  °p  is  reset  to  zero  for  the 
remainder  of  the  run  (in  the  U-Model ,  but  there  is  an  extra  twist  in  the  J-Model). 
The  reset  represents  a  loss  of  prime  in  the  muscle  state. 

To  see  the  rationale  for  Fig.  3,  we  must  sketch  a  bit  of  muscle  physiology 
(cf .  e.g. ,  Stein,  1980) : 

Nervous  action  potentials  into  a  muscle  cell  cause  the  release  of 
calcium  ions  from  the  endoplasmic  reticulum.  These  ions  diffuse 
to  sites  on  actin  filaments  and  thereby  enable  myosin  heads  to 
lock  onto  the  filaments,  pulling  them  toward  the  center  of  the 
sarcomere.  Each  myosin  head  can  repeat  this  action  by  locking 
sequentially  onto  a  series  of  sites  along  an  actin  molecule.  Between 
each  successive  pair  of  such  head  attachments,  a  biochemical  sequence 
of  three  major  actomyosin  molecular  states  must  be  traversed  (Eisenberg, 
et  al.,  1980).  R.  Cooke  (in  private  correspondence  in  1979)  and  others 
have  guessed  that  this  three-state  cycle  may  last  on  the  order  of  50  to 
100  ms.  Thus,  during  a  movement  such  as  ours,  lasting  three  or  four 
hundred  milliseconds,  more  and  more  of  the  most  accessible  (and  perhaps 
most  effective)  attachment  sites  pass  into  restorative  states  that  can¬ 
not  be  influenced  by  calcium  ions  to  produce  force  increments.  As  our 
movement  progresses,  then,  smaller  and  smaller  fractions  of  newly  freed 
calcium  ions  yield  force  increments  before  gathering  again  in  the  reti¬ 
culum.  Moreover,  the  diffusion  times  cf  those  ions  that  do  generate 
forces  presumably  increase,  making  composite  muscle  twitches  caused 
by  later-arriving  nervous  action  potentials  last  longer. 
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We  represent  part  of  these  sketched  effects  in  Fig.  3b)  by  moving  to 

longer  twitch  amplitudes  for  lacer-starting  twitches  (each  twitch  area  being 

normalized  to  nearly  1.0).  The  rest  of  the  saturation  effect  is  represented 

in  Fig.  3c)  by  the  multiplier  exp(-c|).  (The  priming  factor  exp'-°p),  though 

not  yet  mechanically  interpretable,  is  necessary  for  the  models  to  fit  the 

subjects'  data.)  tf£(t)  is  determined  by 
dc*(t) 

— =  -(l-fip)af(t)  +  YFQF(t),  (2.2) 

where  (1-Sp)  (typically  about  0.015)  is  a  rate  constant  for  recovery  of 
toward  zero,  Qp(t)  is  flexor  torque  (approximately  equal  toTp(t)  after  the 
initial  stage  of  the  movement),  and  Yp  (typically  about  1.5)  is  a  rate  constant 
for  the  growth  of  a£(t)  as  a  function  of  Qp(t)  (Qp(t)  maxima  vary  from  0.0172 
in  HintrA  to  0.0342  in  Chan2).  The  c|?(t)  function  in  (2.2)  solves  to 

ft  —  ( 1  — <5p)  (t-x) 

cf(t)  =  YFQF(t)e  dx  (2.3) 

J0  r 

since  c*(0)  =  0. 

This  function  is  the  simplest  one  that  exhibits  the  desired  shape  and  that 
also  gives  good  modelling  results. 

.  Fig.  3  contains  all  flexor  terminaol ogy ,  but  everything  there  applies  also 
to  the  extensor  except  that  r,_  =  0;  the  extensor  prime  is  dissipated  before  the 
onset  of  the  big  extensor  burst  because  the  extensor  must  first  stabilize  the 
elbow  joint  for  the  flexor  to  pull  against. 

3.  THE  U-MODEL 

We  now  give  a  top-down  description  of  our  basic  model,  the  U-Model.  All 
of  its  parameter  values  apply  uniformly  (hence  the  "U")  over  the  entire  dura¬ 
tion  of  each  run.  In  the  J-Model ,  which  is  simply  an  augmented  U-Model,  seme 
parameter  values  are  adjusted  during  each  run  (cf.  the  Fig.  12  caption). 


FASMOV  runs  as  U-Model  when  the  two  parameters  and  T^p  are  set  to  values 
exceeding  the  trial  time  for  the  run  in  question.  Each  run  of  FASMOV  simulates 
one  movement,  or  trial,  of  a  subject. 

Fig.  1  gives  the  elements  of  our  model.  Though  simple,  similar  models 
have  been  employed  by  Messier,  et  al . (1971 ) ,  Sakitt  (1980),  and  others  with 
excellent  results.  Refinements  of  the  Fig.  1  scheme  lead  invariably  to  di stri 
buted  parameter  systems  that  are  very  difficult  to  analyze  in  the  terms  we  have 
chosen,  so  we  have  not  considered  them. 

We  ignored  the  effects  of  gravity  in  our  equations.  This  is  justified  by 
the  high  quality  of  our  results  as  discussed  below. 

The  top-level  analysis  formula  for  our  movement  is 

J+(t)  «  QE(t)  -  QF(t)  -  C?i(t),  (3.1) 

where  each  dot  over  <b  indicates  a  time  derivative,  Qp  and  Q£  are  the  respective 
torques  about  the  elbow  as  exerted  by  the  flexor  muscle  group  and  the  extensor 
muscle,  and  C J>  is  a  viscoelastic  damping  torque  (cf.  Lestienne,  1979).  Since 
the  flexor  group  acts  synergi stically  (cf.  Lagasse,  1975,  e.g.),  we  shall 
regard  it  as  synonymous  with  the  biceps.  J  is  the  arm's  moment  of  inertia  for 
rotations,  (over  a  <*>  change  from  50°  to  about  -30°)  about  a  fixed  elbow  joint 
as  shown  in  Fig.  1.  In  our  model  s’ we  normalized  J  to  1.0. 

At  the  second  level  of  analysis,  we  let 

Qe(t)  *  vVE($(t))  •  *IE(t)  (3.2) 

i>jE  in  (3.2)  may  be  regarded  as  the  maximum  (over  -90°  <  i  <  90°)  isometric 
torque  producible  by  the  extensor,  given  the  smoothed  extensor  EMG  envelope, 
EMGE,  over  the  recent  past.  iyE  is  roughly  interpretable  as  a  f orce-vel oci ty 
factor.  We  got  best  fit  to  our  data  with  ?yE  values  slightly  less  than  1.0 
in  all  cases.  This  is  not  incompatible  with  the  force-velocity  relations 
obtained  by  Thorstensson  et  al .  (1976)  and  Perrine  and  Edgerton  (1978)  in  their 
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studies  of  the  human  knee,  nor  with  Lagasse's  (1975)  work  on  fast  human  forearm 
movements  (quite  like  ours  only  untargeted).  Ultimately,  however,  must  be 
regarded  as  a  successful  coverup  of  the  following  considerations:  Maton  and 
Bouisset  (1977)  found  that  the  muscles  of  the  biceps  group  all  show  a  uniform 
level  of  activity  during  isometric  contraction  regardless  of  muscle  length. 
Messier,  et  al .  (1971)  report  that  under  static  loads,  "the  averaged  EMG  is 
directly  proportional  to  muscle  tension  and  the  constant  of  proportionality  is 
independent  of  muscle  length";  but  Messier's  tension  levels  were  low.  Zuniga 
and  Sirrmone  (1969)  show  that  at  higher  tension  levels  the  tensicn-EMG  relationship 
saturates  quadratically  in  tension.  Certainly  tension  production  in  muscles 
that  are  rapidly  stretching  or  contracting  is  very  different  from  that  in 
muscles  that  are  bearing  static  loads.  During  changes  of  muscle  length,  both 
passive  (visco)  elastic  and  active  molecular  force  components  arise. 

For  the  flexor  side  we  let 

Qp(t)  =  ( "t ) )  •  ^jp(t)  (3.3) 

where  ij^p  is  similar  to  ^  and  =  Vyp-  Equating  ^  to  ^p  for  cur  movement 
is  probably  slightly  inaccurate  because  flexor  and  extensor  muscle  length  change 
in  opposite  directions  and  when  one  is  most  actively  driven  the  other  is  respond¬ 
ing  most  passively  (the  joint  must  be  stabilized,  so  there  are  never  any  purely 
passive  muscle  responses,  only  approximately  passive  ones). 

At  the  third  level  of  analysis  we  let 

t-u  p 

*IF(t)  =3p  cos  ®(t)  X p  EMGp(T)?/p(t-7)  dr  (3.4) 

Jt-uF-?p 

The  integral  in  (3.4)  is  flexor  tension  at  time  (t) ,  and  the  cdp  ccs  4>(  t) 
factor  in  (3.4)  converts5?"^ t)  into  torque.  Initially  in  each  FASMOV  run, 

Q.p  s  Ap,  where  Ap  is  the  length  of  the  effective  flexor  moment  arm  as  shown 
in  Fig.  1.  Whenever  o(t)  changes  from  negative  to  positive  during  a  run. 
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Qp  switches  from  Ap  to  ApS^  as  discussed  in  the  previous  Section.  Whenever 

4> ( t)  changes  back  from  greater  than  zero  to  negative,  as  it  sometimes  does  at 

late  stages  in  a  run,  d p  switches  back  from  ApS^  to  Ap. 

Referring  to  the  Hintrl  run.'n^t-x)  in  (3.4)  is  that  hp(*)  curve  in  Fig.  4 

that  was  in  effect  when  the  muscle  twitch  corresponding  to  the  hypothetical 

EMG(x)  pulse  started  up  from  zero.  This  startup  time  was  t  +  up*  where  yp  is 

the  flexor  electromechanical  coupling  time  (i.e.,  the  time  between  a  nervous 

a 

stimulus  to  a  muscle  and  the  beginning  of  the  corresponding  tension  twitch; 
cf.  Norman  and  Komi  (1979)).  Thus  if  t-T  <  yp,*H‘p(t-t)  is  zero  because  the 


tension  twitch  for  EMGp(x)  has  not  yet  begun.  Letting  denote  the  duration 
of  a  tension-  twitch  starting  up  from  zero  at  x  +  up,  if  t-x  >  yp  +  5p,?fp(t-x) 
is  zero  because  the  tension  twitch  for  EMGp(x)  has  already  ended.  These  zero 


to  Fig.  3. 

For  the  extensor  muscle,  in  place  of  (3.4)  we  have 

ft”wE  _ 

^IE ( t)  C0S  ,  E  EMGE(x)ftE(t-x)  dx  (3.5) 

Jt-yE-CE 

The  explanation  of  (3.5)  differs  from  that  of  (3.4)  in  only  two  respects: 

"°F 

First,  n;E  =  KEe  ,  where  oE  =  a|  withc|  exactly  anal ogous  to  c£;  no  nE 
is  mentioned  because  it  would  always  be  zero.  Second,  QE  starts  equal  to  AE> 
switches  to  Ap/Sp  whenever  $(t)  changes  from  negative  to  positive,  and  switches 
back  to  AE  whenever  #(t)  changes  frcm  positive  to  negative. 

The  fourth  level  of  analysis  involves  detailed  specifications  of  how 
FASMOV  actually  calculates  c-(t)  and  c|(t),  and  how  FASMOV  calculates  or 
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selects  an  appropriate  *H( •)  function  to  use  in  each  of  its  (3.4)  and  (3.5) 
integration  operations.  These  specifications  appear  in  Kilmer,  et.  al .  (1981). 

4.  U-MOOEL  RESULTS  FOR  HINTR 

Fig.  5  shows  the  excellent  agreement  between  U-Model's  computed  response 
for  Hintrl,  and  Hintrl's  recorded  responses  shown  in  Fig.  2.  U-Model's  para¬ 
meters'  were  set  to  produce  Fig.  ‘5  as  noted  in  the  caption.  The  subject's  data 
and  U-Model's  computed  responses  for  the  HintrA  trial  are  presented  in  Figs. 

6  to  7.  The  captions  of  Figs.  5  through  7  cover  several  salient  points  regard¬ 
ing  the  production  of  U-Model's  outputs  for  Hintrl,  and  HintrA. 

Since  the  U-Model  is  basically  an  integration  model,  it  is  not  surprising 
that  Figs.  2  and  6  differ  much  more  between  them  than  Figs.  5  and  7.  If  we 
define  movement  time  as  the  time  between  the  start  of  an  arm  movement  and  the 
first  arrival  at  the  0°  target  position,  Hintr's  movement  times  are  140  ms  in 
both  Hintrl  and  HintrA  (such  consistency  is  certainly  not  the  rule!). 

Fig.  8a)  shows  that  the  second  flexor  burst  for  HintrA  is  appreciably 

smaller  than  the  first.  But  at  the  $  sign  change  at  236  ms  (cf.  Fig.  7),  the 

2 

Qp/Q^  ratio  jumps  to  (S^)  =  9-0  times  its  value  immediately  preceding  236  ms. 

The  validity  of  the  SR  idea  here,  as  suggested  by  Heckathorne  and  Childress 
(1981),  is  strikingly  borne  out  by  the  late  stages  of  U-Model's  response  in 
Fig.  7.  This  excellent  agreement  between  U-Model  and  subject  data  also  suggests 
that  the  confounding  effect  of  gravity  (cf.  Fig.  1),  which  was  not  included  in 
our  analysis,  was  insigif icant. 

Somewhat  surprisingly,  a  constant  electromechanical  delay,  y,  for  both 
flexor  and  extensor  gave  best  results  in  all  our  modelling.  This  is  contrary 
to  measurements  taken  between  isolated  experimental  points  by  Lagasse  (1975), 
for  example,  and  underscores  the  difference  between  continuous  dynamic  and  dis¬ 
crete  initial  parameter  values.  This  difference  is  especially  relevant  since 


Figs.  8  and  11  give  the  trial  record  and  J-Model  outputs  for  two  of  Chan's 
trials.  J-Model  is  a  slight  augmentation  of  U-Model,  as  explained  in  the 
caption  to  Fig.  9. 

As  noted  in  Fig.  9,  Chan's  movements  were  very  fast,  and  the  dynamic  opera¬ 
tional  characteristics  of  his  muscles  differed  markedly  from  those  of  Hintr— 
so  much  so  that  the  augmentation  of  our  U-Model  to  the  J-Model  was  necessary  to 
enable  the  model  to  fit  Chan's  data  at  all  well.  The  caption  in  Fig.  9  describes 
the  augmentation. 

From  Figs.  9  and  11  we  see  that  at  TQE  a  partial  extensor  desaturation 
occurs.  Our  interpretation  of  this  is  that  the  big  extensor  burst  somehow 
partly  "reprimes"  the  triceps. 

At  TQp  a  more  dramatic  desaturation  of  the  flexor  occurs.  This  is  accom¬ 
panied  by  a  sharp  increase  in  the  flexor  tension- to-EMG  scale  factor,  Ap,  and 
a  mild  increase  in  the  flexor  saturation  rate  constant  Yp.  The  large  increase 
in  Reque  (cf.  caption.  Fig.  9)  means  that  the  short  twitch  duration,  (i.e., 
small  value  of  associated  with  the  new  c|  at  Tgp  is  maintained  much  longer 
than  would  be  the  case  with  the  previous  Reque  of  0.3.  Thus  the  second  flexor 
burst  produces  a  predominantly  fast- twitch  effect,  but  also  an  effect  that 
saturates  relatively  quickly.  This  is  discussed  further  below.  Finally,  note 
the  small  38  value  in  Tpp's  summand;  it  implies  that  the  priming  effect  of  the 
second  flexor  burst  became  manifest  in  only  about  two  thirds  of  an  electro¬ 
mechanical  coupling  time. 
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6.  DISCUSSION 

Equations  (3.1)  to  (3.5)  give  the  basic  structure  of  our  model.  In  parti¬ 
cular,  (3.4)  and  (3.5)  show  that  our  major  idea  was  to  impose  ap-  and  c^- 
controlled  nonlinearities  on  hp(*)  and  hp(*)»  but  otherwise  to  compute  the  Qp 
and  Qp  torques  using  the  convolution  integral  of  linear  systems  theory.  Perhaps 
this  notion  has  been  entertained  by  others,  but  we  mathematized  it  and  put 
numbers  into  the  equation.  This  gave  us  empirically  inferred  values  for  some 
skeletalmuscular  functions  that  are  valid  for  at  least  our  movement.  Heretofore 
such  functions  were  determined  only  for  essentially  postural  specifications  or 
for  small  sinusoidal  movements  over  linear  ranges. 

Our  model  assumes  only  feed-forward  control  (contrary  to  the  ideas  of  Sakitt 
(1980)  and  Lestienne,  et  al .  (1979)).  This  raises  the  question:  given  that 
Hintr's  movement  time  was  always  140  ms,  how  could  the  nervous  system  have 
known,  e.g.,  that  an  change  from  0.28  on  trial  10  of  a  block  (Hintrl)  to 
0.32  on  trial  14  of  that  block  (Hintr2  of  Kilmer  et  al .  (1981))  was  needed? 

The  easiest  answer  is  that  the  nervous  system  provided  the  same  "drive  energy" 
in  both  cases,  but  different  muscle  fibers  received  it:  i.e.,  the  spinal  out¬ 
flow  and/or  muscle  inflow  distributions  changed  with  repeated  movements,  but 
the  basic  drive  signals  for  the  movement  did  not. 

This  answer  is*  not  necessarily  applicable  to  the  Ap-to-G^p  change  at  Tpp 
in  Chan's  movements.  There  we  suspect  that  the  nervous  system  fired  more  syn¬ 
chronously  at  Tpp  than  initially  (cf.  Milner-Brown ,  et  al.  (1975))  in  addition 
to  employing  somewhat  different  motor  units  (or  muscle  fibers  within  motor 
units)  at  these  times.  The  model's  parameter  changes  at  Tpp  in  Figs.  9  and  11 
are  what  led  us  to  this  suspicion.  Frequency  analysis  of  raw  EMG's  should  tell 
us  more  about  this  question  soon. 


> 


In  Chan  trials  2,  4,  7,  and  10,  the  movement  times  were  128  ms,  126  ms, 

124  ms,  and  120  ms  respectively.  This  appears  to  be  the  result  of  some  kind 
of  "warmup"  effect  related  to  Chan’s  producing  faster  twitches  throughout  each 
trial  than  Hintr  did.  Note  that  Chan's  initial  Ap  value  changed  from  .28  in 
Chan2  to  .21  in  ChanlO,  whereas  Hintr's  initial  Ap  value  changed  in  the 
opposite  direction  from  .20  in  Hintrl  to  .23  in  Hintr2.  The  foregoing  relation¬ 
ships  may  enable  us  eventually  to  sharpen  and  better  understand  the  correlations 
between  muscle  fiber  types  according  to  biopsy,  muscle  response  speeds  (or 
powers),  and  operant  muscle  twitch  characteristics,  at  least  for  our  type  of 
movement.  Kroll  (1974)  and  others  have  sparse  systematic  evidence  but  abundant 
anecdotal  evidence  that  the  fast-twitch  muscle  fibers  of  biopsy  are  not  always 
the  fast-twitch  ones  of  action.  More  work  is  in  progress  on  this  question. 

Our  model  was  quite  sensitive  to  changes  in  parameter  values,  especially 
those  that  affected  relative  timing.  The  most  disruptive  changes  of  all  were 
unbalanced  flexor  versus  extensor  ones.  The  antagonistic  muscle  systems  in 
our  bodies  must  have  some  way  of  balancing  their  changes  due  to  temperature 
shifts,  etc.;  otherwise  we  could  not  retain  the  shapes  of  our  movements  when 
cold,  hot,  sick,  or  slightly  inebriated. 

We  believe  that  our  model  performs  essentially  the  same  kind  of  integra¬ 
tion  that  our  elbcw  muscles  do,  at  least  when  the  EMG  signals  are  clean  and 
synchronous  as  in  Hintrl  and  A  and  Chanl  and  10.  We  also  believe  that  the 
best-fitting  parameter  values  arrived  at  in  our  Hintr  and  Chan  runs  are  probably 
unique  (to  two  decimals).  Allowing  for  the  coarseness  cf  our  model,  which  was 
retained  to  keep  computer  costs  within  reason,  the  fact  that  cur  model  performed 
so  well  over  many  trials  of  several  subjects  gives  it  a  credibility,  especially 
in  its  temporal  organization,  approaching  validation. 
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7.  CONCLUSIONS 


Our  0-Model  should  prove  useful  in  at  least  two  ways.  First,  it  should 
enable  us  to  check  further  on  the  presently  intrenched  but  increasingly  suspect 
ideas  on  the  role  of  biochemically  classified  fast-twitch  and  slow-twitch 
muscle  fibers  in  fast  movements  such  as  ours.  Some  evidence  to  date  suggests 
that  the  slow-  and  fast-twitch  fibers  of  biopsy  are  not  always  the  slew-  and  fast 
twitch  fibers  of  action.  Now,  by  comparing  J-Model  simulation  with  muscle  biopsy 
results,  we  can  produce  dynamically  valid  correlations  on  this  question. 

A  second  use  of  J-Model  should  be  to  aid  the  design  of  optimal  electrical 
stimulation  patterns  for  use  in  limb  rehabilitation  following  strokes  and  other 
accidents.  This  would  build  on  Boucher  and  Lagasse's  ( 1980 )  success  in  trans¬ 
ferring  a  skilled  muscle  activation  pattern  for  our  movement  type  from  "donors" 
to  unpracticed  "recipients." 
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Fig.  1.  Schematic  representation  of  the  essential  elements  involved  in  our 
arm  movement.  Initially  6  is  60°,  and  the  movement  is  to  a  target  <p  of  0°. 

The  flexor  is  the  biceps  group,  and  the  extensor  is  the  triceps.  Flexor 
torque  is  assumed  equal  to  ^pAp  cos  <?,  where  °Tp  =  flexor  tension;  and 
extensor  torque  is  taken  similarly  as  cos  $. 

Fig.  2.  Plots  from  the  dynograph  recorder  ,n  Hintr's  tenth  trial  on  exercise 
day  four:  The  top  two  traces  are  so-called  "averaged  EMG's,"  and  are  denoted 
EMG(t)  waves,  d)  biceps  EMG;  b)  triceps  EMG;  c)  9  displacement  in  degrees 
(cf.  Fig.  1);  d)  $  in  degrees/sec;  e)  $  in  degrees/sec^ .  d)  and  e)  are  time- 
shifted  to  the  right  because  of  derivational  delays.  Millisecond  time  markers 
are  written  into  the  record  at  selected  places,  starting  with  0  and  the  point 
where  the  a)  plot  rises  off  baseline.  The  small  circles  superposed  onto  the 
a)  and  b)  plots  are  for  deriving  EMG  inputs  to  the  model  from  the  dynograph 
EMG  plots.  All  Hintrl  runs  were  abbreviated  to  340  ms  to  save  computing  costs. 
EMGp(t)  and  EMGp(t)  represent  the  time  courses  of  effective  nervous  drives  on 
the  flexor  and  extensor  muscles  respectively.  The  EMG's  are  essentially  the 


same  as  the  EMG's  except  during  the  EMG's  falling  phases.  The  continuous, 
piece-wise-linear  waves  formed  by  connecting  ad-jacent  pairs  of  encircled  points 
in  a)  and  b)  are  denoted  EMGp(t)  and  EMGp(t)  respectively.  The  EMG's  are  used 

to  produce  the  EMG's  as  per  Kilmer,  et  al .  (1981).  The  idea  is  to  let  EMG  equal 

_  n  _ 

EMG  except  where  EMG  has  a  negative  slope,  and  then  to  let  EMG  be  the  level  to 

which  EMG  appears  to  be  discharging  through  the  dynooraph's  RC  pen-control  circui 


Fig.  3.  A  pictorial  outline  of  the  conceptual  basis  of  the  U-Model .  Tp  is 
the  time  at  which  i  switches  from  less  than  zero  to  greater  than  zero.  SR  is 
the  tension  switching  ratio  for  muscle  lengthening  vs.  shortening.  The  and 
r,p  values  shown  are  fairly  typical.  In  the  Hintrl  run  of  FASMOV,  our  PASCAL 
implementation  of  U-Model ,  =  3.0  and  np  =  -1.0.  The  symbol  means 

"is  set  to". 

Fig.  4.  a,  b,  and  c  show  the  successive  hp(*)  functions  used  from  the  start 
of  the  Hintrl  run  of  U-Model  until  after  a  later  reque  of  hp(*)  that  sets  it 
very  close  to  its  saturated  state,  d,  e,  and  f  show  the  same  for  hp(*)- 

Fig.  5.  The  continuous  curve  is  a  replot  of  Fig.  2c)  frcm  curvilinear  recorder 
coordinates  to  rectilinear  coordinates.  The  small  circles  identify  the  U-Model' s 
output  as  derived  from  its  EMGp(t)  and  EMGp(t)  inputs  (cf.  caption  to  Fig.  2). 

The  parameters  of  U-Model  were  first  set  to  best  fit  Hintrl' s  $(t)  response  data. 
Then  only  A^,  Ap,  and  °p  were  allowed  to  be  reset  to  best  fit  HintrA's  data. 

That  result  is  shown  in  Fig.  7. 

Fig.  6.  Similar  to  Fig.  2  except  this  Figure  is  for  Hintr's  fifteenth  trial 
of  the  second  block  of  trials  on  exercise  day  six. 

Fig.  7.  This  Figure  bears  the  same  relation  to  Fig.  6  that  Fig.  5  dees  to 
Fig.  2.  Because  HintrA's  second  flexor  burst  in  Fig.  6a)  was  relatively  large, 
we  ran  U-Model  for  420  ms  on  HintrA. 


Fig.  8.  Similar  to  Fig.  2  except  this  Figure  shows  Chan's  second  trial  on 
exercise  day  five.  Chan's  movement  times  for  trials  2,  4,  8,  and  10  of  the 
same  block  that  the  above  records  are  taken  from  were  128,  126,  124,  and  120  ms 
(Hintr's  were  140  ms  on  all  trials).  Chan  was  relatively  very  fast,  and  Hintr 
was  about  average. 


Fig.  9.  This  Figure  relates  to  Fig.  8  the  same  as  Fig.  5  does  to  Fig.  2, 
except  that  it  was  produced  by  J-Model  instead  of  U-Model .  The  two  models 
differ  as  follows"  let  the  first  numbers  in  the  Tnc  and  Tnr-  sums  denote  the 
onset  times  of  the  big  triceps  burst  and  the  second  big  flexor  burst  respec¬ 
tively.  The  second  numbers  are  empirical  parameters.  At  TQ^  we  reset  and 
at  Tgp  we  reset  five  model  parameters  as  shown  at  upper  right  in  the  Figure. 
Initial  parameter  settings  are  given  at  lower  left  in  the  Figure.  The  para¬ 
meter  values  below  the  dashed  line  in  the  initial  and  reset  parameter  lists 
are  the  same  forChan2and  for  ChanlO  in  Fig.  11;  only  the  values  above  the 
dashed  lines  differ  between  Chan2  and  Chan  10  runs.  The  value  of  Reoue  is  the 
amount  by  which  c-  (or  c^)  must  change  since  the  last  Reque  operation  (or  since 
t=0)  before  another  aj*  (or  c£)  reque  is  done.  is  discussed  in  Kilmer 

et  al .  (1981);  it  is  used  to  reque  Sp  and 

Fig.  10.  Chan's  trial  10  of  the  same  block  that  the  trial  in  Fig.  8  came  frcm. 

Fig.  11.  This  Figure  relates  to  Fig.  10  the  same  as  Fig.  9  does  to  Fig.  11. 

See  caption  for  Fig.  9  for  comments  on  the  parameter  lists  shown  in  this  Figure. 
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ON  THE  STABILITY  OF  DELAY  EQUATION  MODELS 
OF  SIMPLE  HUMAN  STRETCH  REFLEXES 

ABSTRACT 

Human  stretch  reflexes  (SR.s)  are  often  too  weak  and  ineffectual 
to  provide  adequate  postural  regulation  or  rhythmic  .movement  boosting 
(e.g.  in  ankle  pushoff  at  the  end  of  stance  phase  in  fast  running) . 
Recent  improvements  in  methods  of  artificially  enhancing  skeletomotor 
responses,  especially  in  therapeutic  regimens,  should  not  be  widely 
employed  until  the  clonus-resisting  stability  properties  of  SRs  are 
better  understood. 

We  formulate  an  idealised  linear  servo  model  of  a  segmentaiiy- 
mediated  SR  system  which  includes  the  often  ignored  electromechanical- 
coupling  delay.  For  typical  closed-loop  (delay/gain)  ratios,  the 
model  is  shown  to  be  unstable  for  all  values  of  loop  gain  when  operating 
as  a  position  servo,  but  maximally  stable  when  operating  as  a  velocity 
servo.  We  claim  that  the  velocity  servo  or  one  of  its  nonlinear 
relatives  is  a  better  model  for  some  well  studied  SRs  than,  e.g.,  Houk's 
stiff  muscle  hypothesis.  Iv'e  also  present  evidence  that  even  feeble 
and  quickly  saturating  monosynaptic  postural  servos  are  always  unstable 
if  operated  as  purs  position  regulators. 
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1.  Introduction 


Rack's  (1981)  critique  of  several  papers  on  reflexes  mediated  by 

muscle  afferents  begins  as  follows: 

"The  stretch  reflex  has  often  been  regarded  as  part  of  a 
system  for  the  servo-control  of  limb  position,  through  which 
a  discrepancy  between  the  actual  muscle  length  and  an  intended 
length  generates  activity  that  mini mites  this  error.  If  it  is 
to  work  effectively,  a  servo-control  system  of  this  type  must 
have  a  high  gain  around  the  feedback  loop,  so  that  a  small 
error  leads  to  a  powerful  correction;  all  the  available 
measurements  and  calculations  of  this  loop  gain  do,  however, 
suggest  that  it  is  rather  disappointingly  low  (Matthews,  1972; 
Vallbo,  1974)  and  this  apparent  shortcoming  has  sometimes 
puttied  physiologists .. .We  are  left  therefore  in  a  rather 
unsatisfactory  position;  there  is  ample  anatomical  and  physio¬ 
logical  evidence  that  the  stretch  reflex  pathways  exist,  but  a 
quantitative  examination  seems  to  indicate  that  they  cannot  do 
the  job  that  we  have  frequently  attributed  to  them." 

In  the  same  volume  Hagbarth  (1931),  commenting  on  gamma  loop  reflex 

functions  in  man,  says  that  it  is 

"questionable  whether  the  autogenetic  reflex  support  via  the 
(stretch  reflex)  loop  ever  exceeds  what  is  required  to  counter¬ 
balance  the  force  of  gravity.  'With  available  data  concerning 
the  position  sensitivity  of  spindles  and  their  firing  rate  during 
isometric  contractions,  it  has  been  estimated  that  the  spindle 
inflow  cannot  be  the  main  generator  of  the  motor  resrcnse  that 
enables  the  subject  to  hold  the  joint  position  constant  against 
increased  external  load." 

Similar  appraisals  of  the  mammalian  stretch  reflex  (SR)  in  its  many 
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in  several  professional  journals.  One  wonders  if  SRs  must 
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elements,  namely:  feedback,  gains,  delays,  and  muscle  mechanics,  seems 
apropos.  Though  Evarts  and  Fromm  (1981)  argue  convincingly  that  SRs  in 


intact  primates  must  always  combine  segmental  and  transcortical  con¬ 
tributions  to  the  changes  in  motomeuronal  output  following  limb  dis¬ 
placements,  we  still  have  5SR  loops  whose  elements  are  artifici.  lly 
alterable,  so  our  study  seems  justified  as  a  precursor  to  a  corresponding 
multiple-loop  study  (now  in  progress) . 

It  is  well  known  that  by  sufficiently  increasing  gain-delay  products 
around  simple  closed  loops,  one  can  always  eventually  destabilize  them. 
But  whether  the  many  weak  SSRs  that  have  been  observed  are  close  to 
instability  has  not  been  systematically  studied.  We  shall  formulate  and 
analyze  seme  simple  (but  we  believe  adequate)  53R  models  as  evidence  that 
when  typical  relative  parameter  values  obtain,  the  following  principles 
hold:  1)  If  SSRs  operate  as'pure  velocity  servos,  they  are  always 
unstable.  2)  If  SSRs  operate  as  pure  velocity  servos,  they  are  stable 
provided  the  loop  gain  does  not  exceed  a  critical  value.  Velocity  servos 
could  act  powerfully  enough  to  fit  the  experimental  data  if  their  loop 
gains  decreased  fast  enough  after  the  3SR  response  began.  3)  If  SSRs 
operate  in  equal  proportions  as  the  sum  of  position  and  velocity  servos, 
they  have  stability  properties  approximately  midway  between,  those  of  thei 
two  marts  . 


2.  Models  to  Date 


To  date  only  two  series  of  papers  provide  mathematical  frame¬ 
works  sufficient  to  accommodate  reasonably  definitive  investigations 
of  the  SSR  stability-versus -power  relationship.  Hatte's  series  is 
formulated  in  terms  of  state  equations  for  numbers  of  motor  units  in 
different  contractile  states,  and  contains  too  much  fine-grained  detai 
to  be  appropriate  for  our  purpose  (of,  e.g..  Hates  (1980)).  Stein  and 
Oguttoreli's  (1981)  paper  caps  a  series  that  analytes  loop  responses  t 
central  nervous  system  inputs;  but  they  assume  an  electromechanical 
coupling  delay  of  tero  and  do  their  stability  analysis  as  a  function 
of  delays  in  the  spinal,  cortical,  and  trans -cerebellar  pathways 
(assumed  to  be  25m s,  35ms,  and  8Sms  respecti vely) .  (Cf.  Oguttcreli 
and  Stein  (1976)  for  the  full  stability  analysis.) 

We  consider  SSRs  that  contain  an  electromechanical  coupling  delay 
of  perhaps  50ms  and  a  spinal  loop  delay  (from  muscle  receptor  to 
corresponding  5MG  input  to  the  receptor-bearing  muscle)  of,  say,  25ms. 
We  get  different  results  than  Stein  and  Oguttoreli,  mostly  because  of 
our  different  delay  structure  but  also  because  of  different  parameter 
values  and  muscle  response  shape.  Delays  such  as  we  assume  appear,  fc 
example,  in  Diets  et.  al.'s  (1979)  study  of  gastrocnemius  responses  tc 
foot  contacts  of  sprinting  humans.  The  stance  phase  in  sprinting  last 
only  about  120ms,  but  the  powerful  terminating  pushoff  is  apparently 


:he  result  or  an  active  tension  in  the  gastrocnemius  muse 


an  SSR.  For  other  examples,  we  assume  that  Hcuk's  (19“9)  and  Hcuk  et, 


in  Alluni's  (1981)  study  of  load  compensation  in  human  ankle  muscles 
under  imposed  foot  rotations;  so  this  and  similar  human  SSRs  fall 
within  the  same  general  time  frame.  Some  of  these  are  weak,  and  some 
appear  to  be  moderately  strong,  at  least  if  voluntary  background 
activation  levels  are  high  and  reflex  durations  are  brief  (Gottlieb 
and  Agarwal,  1979) . 

Of  particular  interest  to  us  because  of  our  own  experiments  on 
falling  (cf.  Kroll,  et.al.  (1981))  is  Diett,  et.al.'s  (1981)  study  of 
triceps  SRs  in  humans  landing  from  forward  falls.  They  suggest  that 
the  first  high  force  peak  just  after  contact  with  the  landing  platform 
is  produced  by  a  high  initial  muscle  stiffness  of  Houk  et.al.'s  (1981) 
type,  and  that  from  about  70ms  until  about  200ms  or  so  later,  an  SR 
fully  activates  the  triceps  to  stop  the  fall.  The  first  triceps  EMG 
burst  occurs  at  2Q-50ms,  signalling  an  S3R;  a  second  EMG  peak  occurs 
at  60-80ms  after  touchdown,  and  this  together  with  later  peaks  completes 
the  SR.  The  elbow  angle  sags  from  about  130°  initially  to  about  90° 
finally,  without  significantly  increasing  again  until  a  later  volitional 
resoonse  occurs. 


If  the  oft-invoiked  concept  of  muscle  acting  as  stiff  spring  is  used 
here  to  explain  Diet:  et.al.'s  triceps  response,  why  weren't  there 
significant  elbow  oscillations  at  the  bottom  of  the  fall”  Surely  the 
damping  was  not  great  enough  to  eliminate  them.  And  why  was  the  SR 
response  not  stronger?  Unpublished  results  by  kroll  show  that  electrical 
stimulation  to  the  triceps  just  prior  to  such  falls  can  decrease  elbow 
bending  angles  and  times  to  fractions  of  their  normal  values.  ? echoes 


times  to 


Diets  et.al.'s  SRs  were  not  stronger  because  to  have  then  so  naturally 
would  cause  the  triceps  to  often  cramp  off  unstably.  Could  Diets 
et.al.'s  tricep  SRs  be  better  regarded  as  mediated  by  elbow  angular- 
velocity  servos  of  both  SSR  and  long-loop  types  which  teroed  elbow- 
angular  velocities  as  soon  as  stably  possible  after  hand  contact?  This 
idea  is  compatible  with  the  highly  dynamic  characteristics  of  most 
human  SR  systems.  Or,  could  elbow-position  serves  that  too  feebly 
aimed  to  restore  elbow  angles  at  touchdown  better  explain  Diett,  et.al.'s 
data? 


Any  attempt  to  confirm  the  validity  of  the  above  servo  ideas,  and 
then  use  it  to  artificially  enhance  SSRs  for  whatever  purpose,  should 
proceed  cautiously.  Houk,  et.al.  (1981)  shows  that  small  stretches  give 
high  spindle  sensitivity,  hence  high  loop  gain.  Gottlieb  and  Agarwal 
C 1979)  and  Gurfinksl  and  Shik  (1972)  note  that  voluntary  background 
activity  level  increases  can  multiply  SSR  loop  gains  severalfold. 

Marsden  (1975)  found  that  increasing  the  load  on  a  finger  muscle  increased 
the  associated  SR.  Vallbo  (1981)  suggests  that  this  might  have  been 
caused  by  parallel  activation  of  skeletcmctor  and  fusimotor  systems;  but 
regardless  of  the  mechanisms,  the  net  effect  was  to  increase  SR  loop 
gain.  Proske  and  '<i  alms  ley's  (1931)  study  of  triceps  surae  stiffness  in 
cats  running  freely  on  a  treadmill  shows  that  at  treadmill  speeds  below 
5  a/ s  the  medial  gastrocnemius  (MG)  acts  somewhat  like  a  spring  in  yoke 
with  the  soleus.  3ut  at  higher  speeds,  MG  acts  much  more  powerfully, 
so  its  loop  gain  is  then  correspondingly  higher.  Data  in  Diett,  et.al. 


.979)  an-: 


.1931)  suggest  that  the  same 
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Gurfinkei  and  Shik  (19"Z)  and  A'ieser.dar.ger  (13“3)  point  out  the  complex 


tradeoffs  between  SSRs  and  long- Loop  SRs,  especially  supraspinal!;/ 
altered  SRs.  There  are  also  abundant  references  on  the  distinction 
between  ballistic  reflexes  and  slow-tracking  or  fine-postural  SRs.  Some 
believe  that  SSRs  are  meant  to  function  only  as  fine-postural  position 
servos.  Evarts  and  Fromm  (1981)  opine  that  "segmental  and  transcortical 
proprioceptive  reflexes  allow  for  servo  controls  which  are  important  in 
generating  'mild'  motor  adjustments  necessary  for  precise  active  move¬ 
ment  and  precise  postural  stability,  but  exteroceptive  (rather  than 
proprioceptive)  signals  trigger  the  intense  muscular  reactions  which 
are  often  called  for  in  response  to  stimuli  arising  from  the  environment. 

With  these  caveats  in  mind,  we  shall  formulate  an  SSR  model  that 
will  enable  us  to  investigate  some  SSR  servo  properties. 

3.  The  General  SSR  Model 

We  consider  an  SSR  system  for  extensor  and  flexor  muscle  groups 

about  a  single  joint.  Let  F  and  E  denote  flexor  and  extensor,  y  the  join 
•  • 
angle,  dy/dt  =  y  =  x  the  joint's  angular  velocity,  dx/dt  =  x  the  joint's 

angular  acceleration,  T^  the  load  torque  about  the  joint  (external  torque 

plus  damping) ,  and  Tp  and  T_  the  torques  due  to  flexor  and  extensor 

tension.  Then,  normalising  the  moment  of  inertia  for  rotation  about  the 

joint  to  1,  we  get 

y(t)  =  x(t)  (i) 

x(t)  ~p  ~  -c  ’ 

assuming  that  extensor  shortenings  decrease  y. 

Kilmer  et.al.  (1SS1)  have  shown  that  at  least  for  very  fas*,  large 
amolituds,  start-anc-stoo  arm  movements  ecc-t  the  elbow,  7r  and  T_  can 


be  quite  accurately  expressed  as  a  nonlinear  convolution  integral.  We 


use  here  the  linear  approximation  to  this  integral 


t_  =  f  EMG_cr)  h_(t-r)dr,  (2) 

E  Jt-Ap  c  E 

where  EMG_ )  is  the  envelope  of  the  extensor  E.MG  at  time  T'  ;  h_(t-T) 
is  the  extensor  tension  twitch  at  t  corresponding  to  a  unit  impulse  of 
EMG_  at  (t-  T) ;  and  A  ^  is  the  duration  of  such  a  twitch. 

A  similar  equation  obtains  for  Tp.  Focusing  only  on  the  negative 
feedback  SSR,  and  supposing  that  moment  arms,  etc.,  remain  constant 
throughout  all  rotations,  we  approximate  (2)  by 


rt 

j  <^bEy(r  -  n+cExCT  -S)^>  hcCt-r  )dT 

t-AE 


where  S  is  the  time  delay  between  an  extensor  stretch  and  the  appearance 
on  a  corresponding  extensor  E.MG,  and  ^^^is  0  if  v  is  negative  and  is 
v  otherwise. 


The  assumption  in  going  from  (2)  to  (3) , 
linear  combination  of  y(T  -  S)  and  x(T  -  S'), 
now  we  accept  it  on  the  supposition  that  y  is 
proportional  to  the  difference  between  the  pr 
neutral  reference  length,  and  x  likewise  to  t 
These  approximations  seem  tolerable  in  the  li 
is  to  characterite  the  nature  of  the  relation 
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Our  system  is  now  linear,  so  lr  we  set  n,.  =  h_,  b_  =  b  ,  c_ 

f  t  r  c  i" 


and  T  ='T„  -  T_,  we  can  rewrite  (1)  as 
c  t 

y(t)  =  x(t) 


x(t)  =-  [b/C-r  -  n+cx(T  -  n]  h(t-r) dT  . 


We  assume  in  (4)  that  T.  is  0  by  dint  of  choosing  appropriate  reference 

u 

levels  and  by  taking  damping  to  be  negligible.  We  also  assume  in  (4; 
that  the  quotient  of  muscle  tension  per  unit  EMG  input  in  muscle 
shortening  divided  by  that  in  lengthening  is  approximately  the  same  for 
flexor  as  for  extensor.  The  available  evidence  to  date  supports  this 
assumption  (cf.  Kilmer,  et.al.  (1981)).  Finally,  we  absorb  S  into  the 
h  function  in  (4)  and  change  variables  to  get 

y(t)  =  x(t) 

x(t)  =  -  f  [by  (t+  6  )  +cx(t+  e  fl  h ( 6  )  d  9  (5) 

J  -A. 

>  A 

For  analytical  convenience  we  now  rescale  t  to  t  so  that  A  =  1, 
and  approximate  the  shape  and  duration  of  h(5  )  by 


hce)  = 


(1  +  6  )h  for  -1^  9<  Th  <  0 


for  T,  <  e  <  0 
n 


where  h  is  a  cositive  constant.  In  most  of  what  follows  we  let  T.  =  1/2 
*  n 

with  the  interpretation  that  one  time  unit  on  the  't  scale  is  about  120m 
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eai  twitch  shapes,  but  Kilmer 


et.al.  (1981)  show  that  in  the  highly 


integrative  muscle  component  of  our  system,  such  shape  inaccuracies  are 
not  nearly  as  important  as  twitch  starting  time  and  area,  h  determines 
our  twitch  area,  hence  our  loop  gain.  The  final  rendition  of  our  S3R 
model  is 


y(t) 

x(t) 


(e 


+1 


)  £'°y(Zj 


e) 


►cx(t- 


■•>] 


d  5 


(7) 


=  f(by4  *  c.x*) 
z  c 

where  the  last  line  of  (7)  is  the  shorthand  we  'use  below.  Henceforth, 
we  shall  write  t  simply  as  t. 


4 .  A  Velocity  Servo 

If  b=0  and  c=l  in  (7) ,  we 
.  r  -1/2 

x(t)  =  -h  J  (6 

If  we  assume  that  the  initial 
continuous,  it  follows  from  Ha 
of  Fourier  series  that  all  sol 
are  linear  combinations  of  in 


have  the  velocity  servo 

-i)x(t*  s)de  =  f(xt) 

x(6  )  data  for  0 

is  (1377)  and  the  comp  Is 
utions  to  (3)  over  any  f 
general  infinitely  many 


(3) 


is  piecewise 
ter.ess  properties 
i.nite  t  interval 
terris  of  the  form 


exp  ('7"n  t) ,  where  ~  ^  is  the  nth  complex  eigenvalue  of  (3):  u  ^  =  T 


i/3  ,  n  =  1,2...  .  Thus  all  solution  terms  are  of  the  form 


exp(  V  t)  A  cos  /3  t  *  3  sin/3  cl 
J  n  L  n  /  n 


n  ^  n  _i 


for  acrrooriate  constants  A  and  3  . 
•  -  *  n  n 


; uprose  we  order 
X  • 


sc  rna: 


1  (h)  ^  T.(h)^  X  (h)  >  ... 


th 

Referring  to  the  solution  term  containing  c~^  as  the  n  mode,  anc  that 

containing  <7^  a.s  the  dominant  mode,  we  say  [3)  is  stable  if  Y^(h)  <■  0, 

i.e.,  if  the  dominant  mode  is  nonincreasing. 

To  find  ^(h)  we  substitute  exp(o-  t)  into  (3)  for  x(t)  and  solve 

the  resulting  eigenvalue  equation,  the  real  and  imaginary  parts  of  which 

are,  respectively : 

*  *.  * 

.  p  'S  2-n>  h  (  5  2«-  /S 2)  exp c /2)  ]  /h  = 

-Cl/2)  C  *2+/32)cos(/1/2)  -  'S  cosC/3/2)  -/3sin(/3/2)  (9a) 

-  pcos2(/3/2)  -  2/5cos(./3/2)sinC,(3 /2)-T  sin2(/3 /2)J  /ex?(y/2) 

[2 r/2  C^2w^2)exp(y/2)]  /h  = 

(l/2)sin(/3/2)  -  /5  cos(/?/2)  -  3*  sinC/3  /2)  (9b) 

♦  f/icos2(/3/2)  -  2  Tcos(/3/2)sin(y3 /2)  -/S  sin2(/3 /2)J/exp(y /2) 
The  underlined  terms  on  the  left  side  of  (9a)  and  (9b)  are  referred  to  as 
andip  respectively  below. 

Figure  1  and  Table  1  sketch  the  solutions  to  (9)  as  a  function  of  h. 
As  noted  in  the  Figure  and  Table  captions,  the  two  are  slightly  incompati¬ 
ble  because  of  calculator  roundoff  errors.  It  is  clear,  however,  that  for 
an  h  of  about  20  the  dominant  mode  of  the  solution  to  (3)  goes  unstable. 

In  order  to  translate  this  into  response  shapes  for  different  h  values  in 
a  hypothetical  idealised  falling  situation,  we  computed  the  response 
shown  in  Figure  2  on  a  CDC  Cyber  using  the  highly  accurate  algorithm 
of  Banks  and  Kappel  (19~9) .  For  comparison,  we  also  show  there  as 
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Figure  1.  Sketch  of  the  eigenvalue  loci  in  the  top  half  of  the  or  plane 

for  the  x  =  ex?(<r  t)  solutions  to  x  =  f(x^)  .  The  arrows  on  all  loci 

* 

point  in  the  direction  of  increasing  h,  as  seen  from  Table  1.  The  scale 
on  the  dashed  loci  is  distorted  to  keep  the  loci  in  the  Figure.  The 
bottom  half  of  the  ~  olane  is  the  mirror  image  of  the  top  half. 
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Figure  2.  Solutions  to  (3)  for  different  h  values  as  computed  on  a  CDC 
Cyber  using  Banks  and  Kappel's  (1979)  algorithm.  A  16-spline  approxi¬ 
mation  was  used.  The  initial  data  for  t  <  0  is  interpreted  as  the 
angular  velocity  at  an  elbow  joint  after  touchdown  (at  t  =  -.5)  of  a 
man  falling  as  shown  in  the  stick-figure  insert.  The  hypothetical 
velocity  servo  of  equation  (3)  does  not  develop  any  braking  tension  in 
the  triceps  until  t  =  0  (SCms  after  touchdown  in  the  interpretation) . 
y*(t)  is  a  sketch,  on  a  slightly  different  scale  for  readability,  of  data 
from  Diets  et.al.'s  (1981)  Figure  2.  x*(t)  is  our  crudely  derived 

approximation  to  y(t) ,  and  is  reasonably  accurate  only  for  t  >  0.  Experi 
ments  by  us  indicate  that  our  equation  (3)'s  initial  data  is  a  better 
idealisation  of  what  one  usually  finds. 
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Mote  the  comparison  to  our  x*(t)  sketch  of  the  corresponding  Dietz 
et.al.  (1981)  data.  Obviously,  by  letting  h  in  (3)  decrease  appropriate 
after  t  =  0,  we  could  fit  the  x*(t)  curve  exactly.  Any  reasonable 
interpretation  of  x*(t),  though,  would  have  to  suppose  that  the  initial 
falloff  of  x*  (from  touchdown  at  t  =  060ms  (or  t  =  -0.5)  until  the 
vicinity  of  0  was  due  to  a  passive  viscoelastic  effect,  ar.d  that  only 
the  subsequent  falloff  of  x*  could  have  been  controlled  by  our  hypo¬ 
thetical  velocity  servo.  Given  the  strong,  long-lasting  pre-innervation 
of  the  triceps  before  touchdown  (especially  in  the  blindfolded  subjects) 
a  decreasing  h  after  t  =  0  due  to  a  saturation  effect  seems  plausible 
as  a  way  to  represent  the  decreasing  (tension/EMG)  ratios  that  must 
inevitably  occur  in  such  fast,  high-powered  arm  movements  (cf.  Kilmer, 
et.al.  (1981)).  That  such  an  h  saturation  would  be  stabilizing  to  about 
the  right  degree  can  be  seen  by  comparing  the  shape  of  the  xw  response 
in  Figure  2  with  the  other  responses  there.  We  believe  that  this  is 
the  best  way  to  conceive  of  Dietz  et.al.' s  findings,  not  only  for  the 
S3R  response  component,  but  for  the  entire  SR  response.  We  warn, 
therefore,  of  the  attendant  danger  of  instability  if  artificial  stimu¬ 
lation  is  employed.  Appendix  A  proposes  an  analytically  tractable  way 


of 

letting  h  decrease 
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response 

like  Diet:  et.al.' 

3  . 
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Figure  5.  Solutions  of  x  =  f ^  (x)  for  three  dif rerent  A  s  to  show  the 
stabilizing  effect  of  reducing  the  delay  A  .  h  =  20  for  all  curves. 

The  A  -  1/2  curve  (not  shown) ,  has  an  envelope  whose  magnitude  increases 


very  slowly.  These  curves  were  obtained  using  a  hand  calculator. 
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Table  2.  Similar  to  Table  1  exceDt  now  x 


shows  some  solutions  of  (10)  for  different  h  values  to  convey  this 
more  graphically.  Thus  pure  linear  position  servos  should  never  be 
sought  in  the  body's  melange  of  SSRs.  This  should  eliminate  the  "length 
follow-up  servo"  hypothesis  of  Merton  (1953)  because  that  postulated  a 
fusimotor  equivalent  of  stretch  as  the  actuator  of  a  response  otherwise 
produced  by  a  system  essentially  the  same  as  (10) . 


6 .  A  Combined  Position  and  Velocity  Servo 


=  1 

and  c  = 

1  in  (7) ,  we 

have 

y(t) 

*  x(t) 

(12a) 

-1/2  r 

-1 

X(t) 

a-h  J 

(0  +1) 

-1  L 

/Ct+e  )  +  xCt+e  )j 

de  »  f(yt*xt)  (12b) 

The  eigenvalue  equations  for  (12)  differ  from  (9)  in  that  for  (12)  we  have 

£  =  £(i  2~3  A  2)  Cr  *D  +A2(5  2)  (13a) 

r  (Y+l)2  */32 

S,  =  /3  (3Y  2-fl2)(T  >1)  -  n>  Yf  y2-5/S  2)  (13b) 

(3  *1)"  +A 

Figure  6  and  Table  3  briefly  sketch  the  eigenspace  for  (12)  as  a 
function  of  h.  The  roundoff  errors  are  larger  than  ever,  but  the  quali¬ 
tative  features  in  Figure  5  are  correct.  Figure  7  plots  solutions  to 
(12)  for  several  different  h  values.  Cur  combined  servo  gives  us  an 
expected  blend  of  its  component  servo  response  characteristics.  Different 
proportioned  b  and  c  values  would  simply  weight  the  two  component 
characteristics  differently  to  produce  the  total  system  response. 

The  combined  servo  is  less  stable  than  the  velocity  servo,  but  in 
its  stable  range  of  gains  it  takes  the  rosit’on  error  asymutoticaliy  to 


zero,  albeit  too  slowly  and  with  too  much  oscillation  to  believe  any 
SSRs  actually  operate  that  way.  See  the  graphs  in  Appendix  3. 

Shortening  the  delay  in  (12)  from  1/2  to  &  <(1/2)  has  much  the 
same  effect  as  it  did  for  (3) .  In  fact,  if  h  is  also  greatly  decreased 
(to  slow  the  servo  down),  we  can  use  the  results  of  Driver  (1976)  to 
show  that  (12)  behaves  increasingly  like 

y(t)  *  x(t)  (14) 

x(t)  =  (-h/2) (y(t)  +  x(t) ) 

whose  eigenvalues  are  -(h/4)  ±(h/4)  (1  -4-/~~h/2)  ^  Thus  for  (h/2)"*<i  4, 


(14)  gives  stable  oscillations,  where  otherwise  its  solution's  are  of  the 
form  c^xpC-ct^t)  *■  c2exp(-a  ^t)  ,  with  c^,  C(2  >  0  and  c^,  cn  appropriate 
constants.  It  can  be  shown  in  general  that  reducing  the  product  of  gain 
and  delay  in  (12)  always  has  a  stabilizing  effect. 

Quite  possibly  SSRs  of  the  type  (12),  but  approaching  the  sluggish 
limit  (14),  operate  in  respiratory,  masticatory,  or  other  systems  where 
responses  can  be  slow  and  positions  are  not  always  corrected  by  volitional 


7.  Discussion 

Houk  et.al.  (1981b)  found  that  in  decerebrate  cats  the  discharge 
rate  r  of  spindle  receptors  to  Large  ramp  increases  in  length  obeys 
asymptotically  the  relationship 

r  -  rQ  =  K(y-y1)x0'5  (1: 

where  :<  and  v  are  constants  and  r  is  the  discharge  rate  before 
stretching  begins.  If  (15)  were  proportional  to  our  EMG  in  (I)  and  x 


and  y  were  aiwavs  oositive,  f)  would  change  to 


y.'r;/ JT .  r.jy‘t  ^  y/M  'V ’■.'  r'^vT^io;  ptt^T 


yC«)  =  x(t) 


.  -  -1/2 

x(t)  =  -h  j  K  [y-Ct* e D-xL j  j  xct+e)| 


0 . 5 


(9  +  l)d<5 


(  1 


6) 


Note  that  (16)  assumes  linear  (spindle  activicy)-to-EMG  and  EMG--fco-  m  <4  s  c  I  « 
tension  relations. 

Matthews  (1981)  has  questioned  the  generaiioability  cf  (15)  to 

other  than  ramp  stretches,  but  concedes  that  (r  -  r  )  always  increases 

o 

in  response  to  a  stretch  far  more  slowly  than  in  direct  linear  propor¬ 
tion  to  the  stretch  velocity.  This  is  contrary  to  (7),  but  Houk  et.al. 
(1981)  suggest  that  the  end  result  of  (14)  may  be  a  force-length  relation 
around  the  SSR  loop  which  is  nearly  linear  ever  a  restricted  physio¬ 
logical  range,  in  line  with  his  stiff-spring  hypothesis.  In  fact, 

Gottlieb  and  Agarwal  (1979)  found  a  linear  relationship  between  velocity 
of  imposed  stretch  and  magnitude  of  reflex  EMG  into  human  triceps  surae, 
but  the  brief  10-4Gms  duration  of  their  mvotatic  reflex  was  referred  to 
by  Houk  et.al.  (1981b)  as  a  linearly  scaled  initial  discharge  at  the 
beginning  of  the  latters '  several -second  ramp  stretch. 

We  believe  that  (7)  may  be  an  acceptable  compromise  for  at  least 
a  restricted  range  of  physiological  conditions,  but  we  also  believe  that 
(16)  merits  further  study.  Noting  that  in  Houk's  reply  to  Matthews' 
comments  above,  Houk  said  (15)  does  not  hold  during  muscle  shortening, 
we  set  K  in  (16)  to  \  when  x  ^  0  and  s  when  x<0,  with  \  =  5K 

(the  "5"  here  was  an  educated  guess) .  The  numerical  solution  of  the 
resulting  equation,  starting  from  the  initial  data  cf  Fig.  2,  is  shown 
in  Fig.  7  for  two  different  h  values.  Evidently,  (16)  behaves  essential  1/ 
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Figure  7.  Numerical  solution  to  the  equation  shown,  as  calculated 
with  a  hand  calculator  using  16  steps  per  t  time  unit  and  a 

Simpson's  integration  rule  for  the  right  side  of  the  x(t)  equation 

•  • 

Note  the  different  time  scales  in  the  h  =  10  and  h  =  1  plots. 
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factor  close  to  1).  Thus  (16)  is  probably  unstable  for  all  positive 
values  of  h. 

Since  spindle  firing  rates  reach  their  upper  limit  quite  easily, 
we  augmented  (16)  with  the  tanh  saturating  function  shown  at  lower  left 
in  Fig.  3,  and  again  solved  the  resulting  equation  numerically.  As 
Fig.  3  shows,  this  did  not  change  the  general  character  of  the  solution 
at  all.  Thus  even  mildly  sensitive  and  quickly  saturating  spindle  drives 
on  oc.  -motomeurons  are  very  unstable  in  the  position-like  servo  of  Fig.  3 
This  fact  is  in  stark  contrast  to  what  one  often  reads  in  the  physio¬ 
logical  literature. 

A  future  paper  will  treat  the  extension  of  (7)  to  include  not  only 
the  SSR  loop  but  also  the  two  long  loops  of  Cg uztoreli  (1976).  We  shall 
analyze  the  resulting  system  as  various  combinations  of  velocity  and 
position  servos  operating  with  various  combinations  of  loop  gains  and  as 
perturbed  by  stochastically  shifting  delays.  Oguztoreli  and  Stein  have 
not  considered  such  combinations  and  stochasticity  in  their  series  of 
papers . 

A  natural  test  of  (7)  would  be  to  apply  it  to  the  active  component 
of  the  gastrocnemius  push-off  tension  in  sprinting  human s .  The  challenge 
would  be  to  factor  the  stance-phase  triceps  surae  SR  into  passive  viso- 
elastic  effects,  possible  h  boosting  effects  at  the  beginning  of  each 
step's  gastrocnemius  SSR  response  (especially  likely  in  fast-stretch 
muscular  athletics),  fatigue  effects,  and  the  possible  effects  of 
riy  timed  electrical  stimulation  to  improve  perfor 


oroce 


■nance . 


Figure  3.  Similar  to  Fig.  7  only  the  :c(t)  equation  now  contains  the 
tanh  saturating  function  shown  at  lower  left.  The  h  =  10  scale  appea 
above  the  t  axis  and  the  h  =  1  scale  appears  below  it. 


3.  Conclusion 


Our  linearized  SSR  model,  equation  (7),  is  unstable  as  a  position 
servo,  and  stable  at  sufficiently  small  gains  as  a  velocity  servo.  We 
propose  that  the  velocity  serve  formulation  fits  Diett's  triceps  SR  data 
for  falling  humans  better  than  Houk's  stiff-muscle  hypothesis.  Some 
other  5Rs,  such  as  occur  in  biceps  when  an  arm  in  handshake  position  is 
suddenly  loaded,  should  also  eventually  model  best  as  velocity  servos 
because  probably  they  are  among  the  more  powerful  clonus-free  SR  systems 
of  the  type  described  by  (7) . 

A  position-like  SR  servo  suggested  by  Houk's  power  law  for  spindle 
firing  rates  was  shown  to  be  unstable  even  for  mild  spindle  sensitivities 
and  quickly  saturating  firing  rates,  contradicting  some  of  the  physio¬ 
logical  literature  (but  none  of  Houk's  papers). 

Presently  we  have  strong  analytical  and  simulation 

evidence  that  the  pure  position  servo  of  (10)  with  one  or  two  long  loops 
added  is  unstable  for  all  positive  gain  and  delay  sets,  and  that  the 
corresponding  velocity  servos  are  stable  as  long  as  the  sum  of  the  gain- 
delay  products  for  all  loops  is  less  than  a  fairly  tight  bound.  This 
suggests  the  possibility  that  all  control  of  postural  steadiness  is 

"the  Ta  I  ]  ovifla: 

mediated  by/\  steady,  centrally  produced  signals;  velocity  or  mixed 
velocity  and  position  serves  (accounting  for  the  ubiquitous  10  Hz.  tremor) 
and  volitional  "saccades". 

Gur  results  should  apply  to  the  design  of  functional  electrical 
stimulation  regimens,  anc  to  t.ne  oossibie  future  design  of  muscle-rower 
amplifiers . 
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Aopendix  A 


A  LIAPUNOV  STABILITY  ANALYSIS  FOR  FIGURE  3  ALSO  GIVES  ONE  FOR  .AN 
x* -APPROXIMATING  SYSTEM  FOR  FIGURE  2. 

The  title  of  this  Appendix  states  its  purpose. 

The  system  in  Figure  3  is 

.  *  r  ^ 

x(t)  *  -h  (e  )da  cai) 

J  -1/2- A 

We  call  such  a  system  stable  if  for  every  piecewise  continuous  initial  data 
x(0)  such  that  the  maximum  |x(©  )j  for  0  in  the  interval  from  -1/2-  A 
to  0  inclusive  is  S  ,  there  exists  an  $  (  £  )  such  that  the  maximum  |x(t)| 
over  all  finite  t  ^  0  is  ^  6 (  £  ) . 

In  Section  3.3,  Hale  (1977')  considers  the  equation 
o 


(e  fi)x(t+a  )d3 


x(t)  ■  -  f  a(- $  )g(x(t*S  ))dS 

J  -r 


where  ail  of 

a(r)  *  0,  a(t) ^  0,  a(t)  0,  a(t)  ^  0,  0  ^  t  < 
are  continuous,  and 

x 

G((x)  =  f  g(s)ds  0  as  |  x|-»  *> 

J  o 

He  shews  that  the  function 


r  0 

r  r° 

f 

< 

J  aC'e)  | 

1  J  g(x(t-s))ds 

L 

-j  .-2 

^  - 

(A3) 

! 

is  Liapunov;  that  is,  that 


6 


dV  (x(t))/dt  *  1/2  a(r 


1/2  a(r)  g(x(t*£  ))d^y 

-1/2  f  a(-S  )[  J°  g(x(t+s))dsj  dfi 


is  ^  0  for  all  finite  t >  0  (our  assumption  of  piecewise  continuous  initial 
data  is  an  obvious  generalisation  of  his  continuous  initial  data) .  Thus 
by  his  Corollary  3.1,  (A2)  is  stable. 

If  A’  0,  equation  (A  1)  is  of  the  form  (A2)-(A4);  hence  it  is  stable. 
Moreover,  if  its  initial  data  is  either  always  positive  or  always  negative, 
its  solution  is  easily  shown  to  have  always  the  same  sign  as  its  initial 
data. 

When  A  >  0,  as  it  must  be  in  our  SSR  interpretation,  Hale's 
functional  does  not  yield  (Ao) .  Applying  (A3)  to  (Al)  with  £.  =  1/2, 
r  *  1,  «.(-$)  *  (5+l)h,  a(-6)  *  h,  a(-6)  *  0  and  g(x)  *  x,  we  get 

i  V  (x(t))/<It  •  -1/:  h(-l)  '  x(u)duj  (.47) 

*  s(t)[x(t)-.x(;-i.':)] 

Routine  manipulations  show  that  d  V  (x(t))/dt<0,  implying  that  (Al)  is 
stable,  if 

x(t-l)/x(t) ^  3  h  for  all  t  >  0  (A8) 

(a  loose  sufficiency  condition) . 

Assuming  a  solution  to  (Al)  of  the  form  exp(-«  t) ,  with  (-  o*.  )  the 
dominant  negative  real  eigenvalue  of  (Al)  ,  (AS)  implies  exp(a«.  )  <  3  h, 
where  h  here  cannot  exceed  the  largest  value  remitted  bv  a  rositive  real 


in  (Al)'s  eigenvalue  equation 


2  •  r 

c*  ”  =  h  exp(oc/2)  (ext 


exp  (  c<  /2)  -!)/»■ 


HI 


This  h  is  4.346,  as  seen  from  Figure  1,  and  the  associated  satisfying 
(A9)  is  1.454.  Since  exp(3. 549)  =  3(4.346),  this  associated^  easily 
satisfies  (AS) .  A  simple  argument  also  shows  that  for  initial  data  of 
the  fora  in  Figure  3,  d  V  (x(t))/dt  in  (A")  is  <  0  for  ail  t  >  0.  Thus 
we  have  a  Liapunov  characterization  of  solutions  to  (Al)  that  go 
asymptotically  to  0  as  t  get  large.  This  characterization  will  now  be 
used  to  prove  the  stability  of  a  nonlinear  modification  of  (Al) . 

The  modification  is  given  in  Figure  Al  as  f^(x(t)).  The  idea  of 
f^  is  to  have  a  form  of  the  type  (Al)  but  such  that  h  is  now  an  original 
constant  h  times  ex?(cx(t))  for  seme  real  constant  c.  This  would  give  a 
stable  solution  to  x  =  f^(x)  by  the  above  Liapunov  argument  if  at  each 
x(t)  the  solution's  backward  continuation  over  -l^S^  0  were  exp(cx(t)) 
instead  of  its  larger  true  value.  Replacing  x(t+3)  in  (Al)  by  the 
x(t)  exp(-cx(t) )  9J  =  x(t)  ^exp(- c*t  e  )J  factor  in  f^(x(t) )  accomplishes 
precisely  this  change.  The  results  for  two  values  of  c  are  plotted  in 
Figure  Ai  for  ease  in  visualization. 


Figure  AI.  Solutions  to  x 


=  f^(x)  as  determined  by  hand  calculator  with 
A  t  step  site  equal  to  0.0S.  The  backward  curves  a  and  b  from  x(t^) 
and  x(t2)  with  c  =  -1  portray  graphically  how  the  integrand  in 
decreases  with  decreasing  x. 
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Introduction 


The  terms  muscle  reeducation,  gait  retraining,  and  neurological 
training  commonly  used  in  rehabilitation  exemplify  the  importance  of 
motor  skill  development  in  therapeutic  exercise  regimens.  In  hemiplegic 
patients,  for  example,  rehabilitation  is  considered  to  be  a  learning 
process  involving  reorgani cation  of  the  central  nervous  system  in  order 
to  relearn  lost  functions  (Swenson,  1978)  .  In  cerebral  palsy,  there  is 
a  deficiency  of  motor  control  and  the  goal  of  rehabilitation  is  to  help 
the  patient  acquire  new  motor  skills  (Halpern,  1978) .  It  is  clear,  there¬ 
fore,  that  the  success  of  therapeutic  exercise  systems  is  closely  linked 
to  motor  control  and  motor  learning  knowledge  domains.  That  is,  an  under¬ 
standing  of  the  mechanisms  of  motor  control  and  motor  learning  is  crucial 
to  the  development  and  prescription  of  therapeutic  exercise  systems. 

As  new  information  becomes  available  in  the  neurophysiology  of  movement, 
researchers  in  motor  control  and  motor  learning  attempt  to  assimilate  this 
new  information  into  existing  theories  and  models.  In  some  instances,  the 
new  information  verifies  and  clarifies  existing  theory.  In  other  instances, 
of  course,  the  new  information  necessitates  correction  of  existing  theory. 

In  still  other  instances,  the  new  information  may  suggest  the  wisdom  of 
abandoning  current  theory  and  the  development  of  new  theory. 

Recent  discovery  of  sensory  imparted  learning  has  led  to  the  postulation 
of  a  new  theory  of  motor  learning,  the  reverse  loop  theory  of  motor  learning. 
None  of  the  existing  motor  learning  theories  make  any  provision  for  sensory 
imparted  learning,  suggesting  the  potential  significance  of  the  reverse 
loop  theory  for  all  aspects  of  motor  learning.  The  reverse  loop  theory  of 


motor  learning,  furthermore,  is  not  simply  a  modification  of  existing  motor 
learning  theories,  but  an  entirely  new  type  of  motor  learning.  Coupled 
with  new  information  about  neuromotor  coordination  mechanisms  of  fast  limb 


movements  and  functional  electrical  stimulation,  it  appears  that  a  suitable 
theoretical  interface  can  be  proposed  which  has  considerable  relevance  for 
therapeutic  exercise  systems. 

It  is  the  purpose  of  this  paper  to  present  a  brief  overview  of  the 
relevant  new  developments  in  motor  control  and  learning,  neuromotor 
coordination  mechanisms,  and  functional  electrical  stimulation  as  they 
relate  to  therapeutic  exercise  systems. 

Neurophysiologically  Based  Exercise  Regimens 

A  number  of  central  nervous  system  pathologies  produce  a  variety  of 
residual  motor  deficits  which  are  aggravatingly  resistant  to  rehabilitation. 
Because  of  less  than  satisfactory  results  with  traditional  exercise  regimens 
newer  therapeutic  exercise  systems  have  been  developed  which  are  based 
largely  upon  neurophysiological  mechanisms  underlying  volitional  muscular 
activity.  Analytical  comparisons  of  the  newer  and  more  promising  therapeuti 
exercise  systems  (Harris,  1973;  Hirt,  1967;  Flanagan,  1967)  suggest  that 
these  systems  exhibit  marked  similarities  in  terms  of  actual  techniques 
employed,  and  even  greater  similarities  in  terms  of  central  nervous  system 
effects  being  sought.  As  seen  in  Table  1,  each  of  the  more  recognized 
neurophysiologically  based  exercise  systems  is  characterized  by  a  subset 
of  specific  sensory  stimulation  techniques  capable  of  facilitating  or 
inhibiting  volitional  muscular  action.  Other  systems  utilize  inherent 
movement  patterns  which  are  designed  to  simulate  particular  stages  of 
normal  motor  development. 
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Although  not  explicitly  stated,  the  basic  principle  common  to  all 
neurophysiologically  based  exercise  systems  is  the  goal  of  transmitting 
sensory  stimulus  patterns  which  are  more  intense  than  that  produced  without 
the  added  proprioceptive  facilitation.  Frenkel's  coordination  exercises 
for  tabetic  ataxia  (1902),  for  example,  emphasizes  repetition  as  a  means 
of  producing  desirable  changes  in  the  central  nervous  system: 

The  treatment  for  tabetic  ataxia  is  based  upon  the  education 
of  the  central  nervous  system  by  means  of  repeated  exercises,  where¬ 
by  it  is  enabled  to  receive  sufficiently  distant  stimuli  from  the 
limbs . . . 

Theoretically,  the  transformation  of  an  ataxic  movement  into 
a  normal  movement  takes  place  in  tabetic  subjects  according  to  the 
same  laws  as  the  acquisition  in  healthy  persons  of  a  complicated 
movement . . . 

Thus,  Frenkel's  system  emphasizes  repetitive  movement  as  a  means  of 
transmitting  appropriate  sensory  stimulus  patterns  to  the  central  nervous 
system  with  such  repetitive  movements  chosen  to  simulate  fundamental  stages 
of  normal  motor  development  in  an  appropriate  sequence.  Quite  similar 
rationales  are  espoused  by  treatment  regimen:,  prescribed  by  Fay  (1954), 

Rood  (1954),  Phelps  (1967),  and  Doman-Delacato  (Delacato,  1966:  Doman  et  al. 
1960)  . 

The  Bobath  system  (1970)  utilizes  various  postural  reflexes  and 
equilibrium  reactions  "...for  purposes  of  modifying  muscle  tone  or  eliciting 
desired  movements  (Harris,  1978)."  Muscle  stretch  and  tapping  are  used  to 
further  enhance  proprioceptive  facilitation.  Herman  Rabat  was  one  of  the 
first  proponents  of  proprioceptive  neuromuscular  facilitation  as  a  "... 
valuable  tool  in  therapeutic  exercise  to  increase  central  excitation  and 
thereby  enhance  the  voluntary  activity  of  paralyzed  or  paretic  muscles 


(Knott  and  Voss,  1962)."  Included  in  Rabat's  suggestions  were  techniques 
utilizing  reflexes,  irradiation,  and  successive  induction. 

It  seems  clear,  therefore,  that  the  production  of  enhanced  sensory 
stimulus  patterns  for  transmission  to  the  central  nervous  system  is  the 
goal  of  all  modem  therapeutic  exercise  systems  based  upon  neurophysiological 
mechanisms.  Implicit  in  such  a  treatment  rationale  is  the  importance  of 
transmitting  appropriate  sensory  stimulus  patterns  of  enhanced  intensity. 
Beneficial  effects  of  such  heightened  sensory  stimulation  transmitted  to 
the  central  nervous  system  are  presumed  to  enhance  recovery.  The  exact 
neural  mechanisms  involved  remain  poorly  defined  but  most  theoretical 
explanations  include  some  variation  of  the  concept  of  neuronal  plasticity 
(Bach-y-Rita,  1980;  Cotman,  1978)  .  The  neuronal  plasticity  argument,  of 
course,  purports  that  the  central  nervous  system  has  the  capacity  to 
create  alternate  pathways  and  thus  compensate  for  damage  to  previously 
established  central  nervous  system  pathways  for  volitional  muscular  activity. 

It  was  once  held  that  damage  to  central  nervous  system  pathways  and 
centers  was  irreversible.  Evidence  has  accumulated,  however,  which  suggests 
that  recovery  of  central  nervous  system  functional  capabilities  may  occur 
in  several  ways.  Cotman,  Nieto-Sampedro,  and  Harris  (1981),  for  example, 
have  shown  that  actual  formation  of  synapses  may  result  from  repeated 
stimulation,  and  suggested  that  such  "reactive  synaptogenesis"  may  be  a 
mechanism  involved  in  recovery.  Synaptic  plasticity  can  also  manifest 
itself  by  actual  morphologic  changes  in  synaptic  structures  as  Eccles  (1979) 
has  shown  an  increase  of  up  to  -10  percent  of  synaptic  spines  on  dendrites 
of  granule  cells  in  hippocampal  neurons.  More  extensive  dendritic  trees 


of  layer  II  pyramidal  neurons  have  been  described  in  the  elderly  (Buell 
and  Coleman,  1976)  suggesting  plasticity  even  in  the  aging  brain.  As 
early  as  1902,  Lazarus  speculated  that  recovery  of  function  in  hemiplegia 
might  result  by  compensatory  utilization  of  different  nerve  pathways  which 
take  over  for  damaged  pathways.  Both  Wall  (1980)  and  Bach-y-Rita  (1981) 
suggest  that  existing  but  unused  pathways  may  be  "unmasked"  when  the 
ordinarily  used  pathways  have  been  damaged.  Such  unmasking  is  enhanced 
by  an  increased  central  excitability  due  to  decreased  sensory  input. 
Proprioceptive  facilitation  is  also  known  to  allow  the  recruitment  of 
"subliminal  fringe  neurons"  which  cannot  be  volitionally  activated.  It 
should  be  noted  that  the  sensory  stimulus  pattern  transmitted  to  the  central 
nervous  system  plays  an  essential  role  in  all  of  the  proposed  mechanisms 
for  plasticity  described  above. 

Equally  obvious  is  the  observation  that  a  limb  (or  any  movement 
synergy)  affected  by  some  pathological  condition  can  produce  only  inefficient 
and  incorrect  volitional  movement.  Such  inefficient  and  incorrect  volitional 
movement  produces  a  sensory  stimulus  pattern  for  transmission  to  the 
central  nervous  system  which  is  itself  inefficient  and  incorrect.  Addition 
of  proprioceptive  stimulation  techniques  may  improve  the  quality  of  the 
movement  produced,  but  only  to  a  small  degree.  Thus,  even  with  the  benefit 
of  proprioceptive  neuromuscular  stimulation,  the  central  nervous  system  still 
receives  an  impoverished  sensory  stimulus  pattern.  If  recovery  of  central 
nervous  system  functioning  is  dependent  upon  the  quality  of  the  sensory 
stimulus  pattern  received,  then  recovery  is  likely  to  be  slow  and  dependent 
upon  gradual  improvement  in  the  quality  of  the  movement  synergy  which 
produces  slightly  improved  sensory  stimulus  patterns. 


What,  however,  would  be  the  effect  upon  central  nervous  system 
functioning  if  the  sensory  stimulus  pattern  transmitted  to  it  from  an 
afflicted  limb  was  not  only  correct  and  efficient,  but  an  optimized 
sensory  stimulus  pattern  representing  "perfect"  limb  movement  patterns? 

Such  a  possibility  now  seems  feasible  because  of  developments  in  functional 
electrical  stimulation,  in  the  understanding  of  neuromotor  coordination 
mechanisms  for  fast  human  limb  movements,  and  in  the  demonstration  of 
sensory  imparted  learning.  As  will  be  shown,  each  of  these  three  relatively 
distinct  areas  of  research  have  certain  elements  of  relevance  to  each  other, 
and  a  blending  of  such  elements  may  hold  promise  for  a  vastly  improved 
rehabilitation  regimen  for  recovery  of  function  in  afflicted  limbs. 

Functional  Electrical  Stimulation 
In  46  A.D.  a  Roman  physician  named  Largus  used  the  electrical  discharge 
of  torpedo  fishes  in  the  treatment  of  pain  caused  by  headache  and  gout. 

A  rival  Greek  physician,  Dioscorides,  discounted  the  value  of  torpedo 
fishes  for  headache  relief  but  recommended  its  efficacy  in  the  treatment 
of  hemorrhoids.  In  1745,  the  Leyden  jar  and  an  electrostatic  generator 
were  used  to  treat  angina  pectoris,  epilepsy,  hemiplegia,  kidney  stones, 
and  sciatica.  Benjamin  Franklin  used  an  electrical  device  to  successfully 
treat  a  young  woman  suffering  from  convulsive  fits.  Electrical  stimulation 
has  also  been  used  in  hydroelectric  baths  to  cure  chronic  inflammation  of 
the  uterus,  and,  in  the  1920s,  used  in  conjunction  with  a  spiritualist  to 
drive  out  a  possessing  spirit  from  the  body  of  a  mentally  disturbed  patient. 


Fortunately,  such  an  amusing  early  history  of  the  simple  beginnings 
of  electrical  stimulation  did  not  prevent  eventual  development  of  a  wide 


variety  of  useful  therapeutic  and  rehabilitative  techniques.  Cardiac 
pacemakers  and  transcutaneous  electrical  nerve  stimulation  (TENS)  for  the 
treatment  of  pain,  for  example,  are  two  such  valuable  developments. 

Another  electrical  stimulation  technique  receiving  major  attention  is 
that  of  functional  electrical  stimulation  (FES)  .  Work  in  FES  involves 
stimulation  of  the  nervous  system  or  the  muscle  in  order  to  restore  lost 
function.  Liberson,  for  example,  developed  both  the  technique  and  the 
device  allowing  stimulation  of  the  peroneal  nerve  to  activate  the  anterior 
tibialis  muscle  and  help  restore  a  normal  gait  pattern  in  hemiplegic 
patients  with  a  foot  drop  problem. 

At  present,  two  major  centers  of  research  and  development  are  active 
in  the  field  of  FES.  Under  the  leadership  of  Lojze  Vodovnik,  the  Faculty 
of  Electrical  Engineering  at  the  University  of  Ljublana  in  Yugoslavia  have 
concentrated  their  efforts  in  the  engineering  development  of  electrical 
stimulators  to  be  used  as  prosthetic  devices.  The  Rehabilitation  Engineerin 
Center,  Rancho  Los  .Amigos  Hospital,  University  of  Southern  California  has 
a  similar  focus  with  added  emphasis  upon  treatment  and  the  education  of 
practitioners  in  the  use  of  FES  techniques.  FES  is  currently  being  used 
in  the  treatment  of  patients  with  stroke,  arthritis,  spinal  cord  injury, 
and  post-operative  orthopedic  problems. 

FES  and  Neurological  Recovery 

Limb  paralysis  is  included  among  the  numerous  ailments  for  which 


electrical  stimulation  has  been  prescribed  ever  since  Johann  Gottlog  krueger 
speculated  about  its  use  for  restoration  of  function  in  paralyzed  limbs 
in  1743  (Licht,  1959).  As  described  by  McNeai  (197~)  ,  kratzenstein  in  174-1, 


Jallabert  in  1748,  and  Quelmalz  in  1753  all  reported  successive  restoration 
of  function  in  paralyzed  limbs  by  the  use  of  electricity.  The  modem  era 
of  the  use  of  electrical  stimulation  seems  to  have  begun  with  classical 
study  reported  by  Liberson  and  his  colleagues  in  1961 ^ .  In  1961,  Liberson 
and  associates  successfully  developed  a  technique  which  stimulated  the 
peroneal  nerve  when  the  heel  of  the  foot  of  a  hemiplegic  patient  touched 
the  floor.  Stimulation  of  the  peroneal  nerve,  activated  the  anterior 
tibialis  muscle  and  helped  restore  a  normal  gait  pattern  to  hemiplegic 
patients  with  a  foot  drop  condition.  It  was  observed  that  even  without 
benefit  of  the  electrical  stimulation  "...the  patients  acquire  the  ability 
of  dorsiflexing  the  foot  by  themseives."  Although  Liberson  and  his  associates 
noted  that  the  positive  effects  were  transitory,  the  concept  of  a  "functional 
electrotherapy"  in  which  electrical  stimulation  of  a  muscle  could  replace 
impaired  innervation  was  bom. 

Other  investigators  (Camstam,  Larsson,  5  Prevec,  1977;  Dimitrijevic, 
Gracanin,  Prevec,  5  Trontelj ,  1968)  have  described  apparent  retraining  or 
reorganization  of  the  central  nervous  system  due  to  repetitive  functional 
electrical  stimulation,  or  FES  as  it  is  now  known.  Gracanin  and  Marincek 
(1970)  reported  long  lasting  effects  from  peroneal  nerve  activation  in 
hemiplegic  foot  drop  patients  to  the  extent  that  many  patients  eventually 
walked  without  electrical  stimulation.  Waters  et  al.  (1975)  also  reported 
permanent  changes  in  walking  patterns  with  positive  changes  in  stride  length, 

^Zankel  (1960)  reported  successful  treatment  of  hemipiegic  upper  limbs  with 
stimulation  assistive  exercise  (5AE)  in  which  repetitive  sinusoidal  electrical 
stimulation  of  paralyzed  muscle  was  dene  in  conjunction  with  the  patient's  own 
volitional  exercise  effort. 


single  stance  duration,  and  walking  speed  while  Gracanin,  Vrabic,  and 
Vrabic  (1976)  showed  positive  changes  in  walking  and  posture  in  cerebral 
palsied  children  even  after  cessation  of  electrical  stimulation.  Vodovnik 
and  Rebersek  (1973)  demonstrated  greater  voluntary  torque  production  for 
dorsiflexion  when  stroke  patients  had  volitional  effort  coupled  with  con¬ 
current  electrical  stimulation.  The  latter  study  is  quite  similar  to  the 
stimulation  assistive  exercise  (SAE)  reported  by  Zankel  in  1960  which 
coupled  volitional  exercise  assisted  by  electrical  stimulation  with  the 
result  of  improvement  in  range  of  motion  of  hemiplegic  upper  limbs. 

Available  evidence  concerning  the  use  of  functional  electrical 
stimulation  as  a  technique  for  neurological  retraining  is  not  plentiful, 
and  not  without  methodological  deficiencies.  Improvements  cited  previously, 
for  example,  could  have  resulted  from  improvements  in  muscle  strength 
and  endurance  rather  chan  from  changes  in  the  neural  control  mechanisms. 

In  several  of  the  functional  electrical  stimulation  studies  the  methodology 


employed  could  be  likened  to  proprioceptive  stimulation  techniques;  i.e. 
superimposing  an  electrical  stimulation  upon  a  volitional  contraction 
(Liberson  et  al . ,  1961;  Vodovnik  and  Rebersek,  1973;  Zankel,  1960)  may  have 
acted  much  like  any  other  proprioceptive  stimulation  technique  to  enhance 
volitional  muscular  activity. 

As  it  has  been  used  to  date,  functional  electrical  stimulation  has 
been  administered  either  to  activate  a  paralyzed  muscle  to  produce  action 
where  none  can  be  produced  volitionally,  or  to  provide  additional  stimulation 
to  an  active  but  weakened  muscle.  Improvement  in  neural  control  mechanisms, 


however,  cannot  be  separated  from  imorovements  in  muscle  strength  and  enduran 


because  no  measures  of  neural  control  mechanisms  have  ever  been  collected. 
Since  neurological  retraining  is  of  more  importance  than  changes  in  the 
peripheral  muscle  apparatus,  there  is  an  obvious  need  to  develop  a  research 
methodology  which  can  differentiate  true  changes  in  neural  coordination 
mechanisms  from  changes  in  the  periphery.  Through  no  accident,  planned 
research  strategies  have  been  developed  which  provide  both  the  technique 
and  the  basic  knowledge  necessary  to  distinguish  true  changes  in  neural 
coordination  mechanisms  from  changes  in  the  peripheral  muscle  site. 


Neuromotor  Coordination  Mechanisms 


Muscular  coordination  involves  the  participation  of  many  muscles  that 
essentially  act  in  cooperation  with  each  other  to  apply  force  and  produce 
movement,  or  in  opposition  to  each  other  to  inhibit  force  application  and 
impede  movement.  In  effect,  the  properly  timed  contraction  of  muscles  at 
specified  force  levels  involving  the  agonist,  antagonist,  and  synergist 
muscles  must  be  done  in  a  sequential  order  for  movement  to  even  occur  and 
for  efficiency  of  movement  to  be  possible.  The  initiation  of  a  movement, 
control  during  the  movement,  and  the  ending  of  the  movement  all  depend  upon 
basic  neuromotor  control  mechanisms.  It  is  oovious,  therefore,  that  the 
timing  of  force  applications  by  involved  muscle  groups  is  an  inherent 
feature  of  skilled,  coordinated  human  performance. 

One  of  the  more  persistent  lines  of  investigation  dealing  'with 


coordinated  movement  and  application  of  muscle  force  has  been  that  dealing 
with  maximum  speed  of  limb  movement.  Since  the  force  required  tc  cause 
movement  of  a  limb  is  proportional  to  the  distance  travelled  over  time 
squared,  F  =  2md/t“  or  F  =  ma,  it  was  believed  that  the  greater  the  static 


force  (F)  of  a  muscle,  the  greater  the  speed  of  limb  movement  which  could 
be  produced.  'When  studied  in  humans,  however,  the  maximum  speed  of  limb 
movement  bore  no  relationship  to  the  maximum  isoemtric  strength  of  the 
agonist  muscle  group. 

Although  the  study  of  maximum  strength  of  an  agonist  muscle  group 
and  maximum  speed  of  movement  has  gone  on  for  several  decades,  the  present 
status  of  our  understanding  of  this  basic  coordinated  movement  pattern 
remains  painfully  inadequate.  An  enormous  amount  of  literature  exists 
in  support  of  the  untenable  notion  that  human  movement  does  not  obey  the 
basic  principle  of  physics  that  F  =  ma  since  maximum  strength  of  an  agonist 
muscle  group  simply  does  not  correlate  with  maximum  speed  of  limb  movement. 
Henry  and  Whitley  (1960),  Henry  (1960,  1962),  Clarke  (1960),  Smith  (1961a, 
1961b,  1969),  Rasch  (1963),  and  Nelson  and  Fahrney  (1963),  for  example, 
all  studied  the  speed  of  various  limb  movements  and  reported  non-significan 
correlations  with  maximal  strength  of  the  agonist  muscles  involved  in 
producing  the  movement. 

Work  in  our  laboratory  has  approached  the  problem  differently  and  we 
are  now  able  to  show  that  agonist  muscle  strength  does  not  correlate  with 
maximum  speed  of  movement  because  the  neuromotor  coordination  mechanism 
has  been  overlooked  by  past  researchers.  Any  attempt  to  understand  the 
mechanism  involved  in  fast,  ballistic,  coordinated  movement  must  consider 
the  sequential  timing  of  activation  of  agonist  and  antagonist  muscle  groups 
which  control  the  application  of  force  and  produce  any  resultant  movement. 

In  any  movement  of  a  limb,  the  force  exerted  by  the  prime  movers  must 
be  large  enough  to  overcome  the  inertia  of  the  limb  and  the  force  act  lied 


by  the  antagonists.  With  the  prime  movers  in  control,  movement  of  the 


limb  occurs  as  characterized  by  displacement,  acceleration,  and  velocity 
parameters.  Limb  movement  decelerates  when  the  force  applied  by  the 
antagonist  is  large  enough  to  overcome  the  force  of  the  prime  movers. 

Most  biomechanical  analyses  of  limb  speed  of  movement  have  carefully 
measured  acceleration,  velocity,  point  of  deceleration,  and  total  movement 
time.  Unfortunately,  the  biomechanical  approach  has  seldom — if  ever — 
sought  an  explanation  of  the  motion  being  analyzed  in  terms  of  mechanisms. 
Traditional  speedrStrength  investigations  have  limited  themselves  to  the 
study  of  the  relationship  between  maximum  strength  of  agonist  muscles  and 
limb  speed  of  movement  parameters.  Even  more  surprising,  the  sequential 
order  of  muscle  activation--which  in  effect  produces  and  controls  movement — 
has  been  similarly  overlooked. 

As  long  ago  as  1677  Descartes  postulated  that  the  interaction  between 
muscles  was  regulated  by  a  series  of  valves  which  adjusted  the  flow  of 
animal  spirits  into  flexor  and  extensor  muscles.  The  question  of  actual 
interaction  between  agonist  and  antagonist  muscles  in  humans  has  remained 
a  key  issue  in  the  neurophysiology  of  movement  ever  since.  Boulogne  (1367) 
and  Winslow  (1732)  believed  that  agonist  and  antagonist  muscles  contracted 
simultaneously  while  Bell  (1323)  and  Pettigrew  (1923)  contended  that  the 
contraction  of  one  muscle  was  accompanied  by  relaxation  of  its  antagonist. 

Later  work  by  Beaunis  (1389)  and  Demeny  (1390)  utilized  kymographic 
records  and  supported  simultaneous  contraction  of  agonist  and  antagonist 
muscles.  Sherrington  (1906),  however,  promulgated  his  famous  principle  of 
reciprocal  innervation  and  showed  that  contraction  of  an  agonist  was 
accompanied  by  inhibition  of  the  antagonist  in  laboratory  preparations. 


Tilney  and  Pike  (1925)  studied  many  different  actions  in  the  upper  limbs 
and  were  able  to  verify  the  principle  of  reciprocal  innervation  in  voluntary 
human  movement  as  they  showed  both  agonist  and  antagonist  active  at  the 
same  site. 

The  first  electromyographic  investigation  of  the  problem  by  Golla  and 
Hettwer  (1924)  and  later  electromyographic  studies  by  Stetson  and  Bouman 
(1935),  Hudgins  (1939),  and  Bierman  and  Ralston  (1963)  suggested  that  the 
type  of  movement  being  studied  was  important.  Antagonist  muscles  were 
usually  found  to  co-contract  during  slow  controlled  movements,  and  found 
not  to  co-contract  during  fast  ballistic  movements.  Additional  research 
on  the  interaction  between  agonist  and  antagonist  muscles  in  voluntary 
movements  supports,  in  general,  the  phenomeon  of  co-contraction  and  that 
a  relationship  exists  between  the  speed  at  which  a  limb  is  moved  and  the 
sequential  order  of  muscle  firing  patterns. 

A  better  understanding  of  the  neuromotor  control  mechanisms  involved 

in  limb  speed  of  movement  came  about  when  the  research  methodology  of 

biomechanical  analysis  of  motion  was  combined  with  the  research  methodology 

of  the  neurophysiology  of  movement.  The  first  comprehensive  study  assessing 

the  neurophysiology  and  the  biomechanics  of  limb  movement  was  done  in  our 

✓ 

laboratory  as  part  of  a  Ph.D.  dissertation  (Lagasse,  1973)  .  By  monitoring 
the  time  of  arrival  of  nervous  impulses  to  the  biceps  and  triceps  muscles, 
the  initiation  of  a  rapid  forearm  flexion  movement,  velocity,  acceleration, 
and  total  speed  of  movement,  Lagosse*  was  able  to  show  that  the  sequential 
order  of  agonist-antagonist  muscle  activation  was  a  significant  predictor 
of  speed  of  movement;  the  multiple  R  of  .63  being  the  highest  correlation 


ever  shown  between  limb  speed  of  movement  and  a  criterion  measure  related 
to  a  probable  mechanism.  The  addition  of  acceleration  time,  which  alone 
correlated  with  speed  of  movement  r  =  .79,  resulted  in  a  multiple  R  of  .90. 
Lagasse  also  showed  that  practice  influenced  the  timing  and  coordination 
of  the  neuromuscular  control  mechanism,  largely  by  a  change  in  the  agonist- 
antagonist  sequential  innervation  pattern. 

Individuals  capable  of  the  fastest  speed  of  limb  movement  do  not 
contract  agonist  and  antagonist  muscles  simultaneously  to  the  degree  that 
slower  subjects  do.  Marked  co-contraction,  therefore,  characterizes  poor 
neuromotor  control  while  definite  separation  of  agonist  and  antagonist 
muscle  contraction  patterns  with  a  minimum  of  overlap  characterizes  skilled 
neuromotor  control.  Changes  in  the  agonist-antagonist  contraction  pattern, 
furthermore,  actually  constitutes  a  pure  measure  of  motor  learning  reflecting 
an  inherent  nervous  system  capability  to  learn  which  might  very  well  be 
related  to  a  number  of  other  motor  performance  qualities  including  motor 
educability,  the  development  of  power,  and  the  capacity  to  learn  motoric 
skill  tasks. 

As  part  of  a  grant  from  the  U.  S.  Army  Medical  Research  and  Development 
Command  ( DAMD  17-30-0101) ,  we  have  been  investigating  the  neuromotor  coordination 
mechanisms  involved  in  rapid  limb  movements.  Mathematical  models  have  been 
developed  using  a  two  compartment  representation  of  the  neuromuscular  system 
involved  in  voluntary  fast  arm  movements  to  a  target.  The  first  compartment 
accepts  averaged  biceps  and  triceps  EMG  signals  as  inputs,  and  models  the 
arm's  musculo-skeletal  response  by  predicting  angular  position Q  (t)  and 
velocity  Jl  (t)  over  the  corresponding  movement  time.  This  is  called  the 


E/Q/<j)  (t)  model,  short  for  EMG/elbow  torque  =-Q/elbow  angular  velocity 
dQ/'dt  =  fj>,  and  its  defining  equation  is  of  the  form: 

|)  (t)  =  (Extensor  torque  Q£)  -  (Flexor  torque  Qp)  . 

The  second  compartment  of  the  model,  called  the  command/neural 
control/EMG  or  Com/Cont/E  model,  accepts  the  volitional  commands  as  inputs, 
and  models  the  nervous  system  response  by  producing  a  smoothed  flexor 
and  extensor  EMG  signals  to  feed  into  the  E/Q/|)  model.  Input  commands 
are  presumably  via  the  pyramidal  tract  from  the  cerebellum  to  the  brain¬ 
stem  and  spinal  cord  (cf.  .Miles  and  Evarts,  1979). 

The  purpose  of  the  two  compartment  model  is  to  separate  various 
control  functions  in  the  ballistic  arm  movement  system  well  enough  to 
understand  where  adaptation  occurs  as  speed,  precision,  and  coordination 
changes  occur  due  to  practice.  The  appropriate  connection  of  the  two 
compartment  model  provides  a  complete  model  for  fast  limb  movement  as 
driven  by  volitional  commands  (Kilmer,  Kroll,  and  Congdon,  1981). 

Analysis  of  the  empirical  data  on  limb  movement  and  parameters  in 
the  mathematical  model  describing  neuromotor  coordination  mechanisms  shows 
quite  clearly  that  two  major  mechanisms  are  operative.  The  first  mechanism 
involves  the  temporal  sequence  of  agonist  and  antagonist  firing  patterns. 

If  the  agonist  begins  its  contraction  without  co-contraction  of  the 
agonist,  faster  limb  movement  is  produced.  Antagonist  muscle  firing, 
however,  is  needed  to  allow  precision  of  movement  to  a  target.  When  the 
antagonist  fires  early  so  as  to  overlap  with  agonist  firing,  limb  movement 
is  slowed.  If  the  antagonist  fires  later,  thus  allowing  agonist  contraction 
to  manifest  itself  without  competition,  acceleration  continues  longer  and 


faster  limb  movement  is  produced.  Indeed,  one  of  the  marked  practice 
effects  found  is  that  antagonist  firing  occurs  later  during  the  movement. 

The  second  major  mechanism  involves  the  ratio  of  flexor  to  extensor 
EMG  intensity  and  duration.  It  is  self-evident  that  the  flexor  EMG 
intensity  and  duration  must  be  greater  than  comparable  extensor  EMG 
activity  if  a  limb  flexion  movement  is  to  be  produced.  As  we  have  shown 
(Kilmer,  Kroll,  and  Congdon,  1981;  Kroll  and  Kilmer,  1981),  this  flexor/ 
extensor  EMG  ratio  is  as  critical  as  the  temporal  firing  sequence  of 
agonist  and  antagonist  muscles  in  producing  fast  and  efficient  limb 
movement . 

Thus,  the  temporal  sequence  of  agonist  and  antagonist  firing,  and  the 
flexor  to  extensor  EMG  ratio  intensity  and  duration  are  the  most 
critical  neuromotor  coordination  mechanisms  involved  in  limb  movement 
speed.  Although  our  analyses  suggest  the  existence  of  other  factors 
in  limb  movement,  such  as  muscle  fiber  type  composition  and  synchronized 
versus  asynchronized  motor  unit  firing  patterns,  the  two  major  mechanisms 
identified  and  quantified  are  capable  of  excellent  prediction  of  limb 
movement.  The  accuracy  with  which  the  mathematical  model  can  predict 
limb  movement  can  be  appreciated  by  an  inspection  of  Figure  1.  In 
Figure  1,  the  solid  line  represents  actual  limb  flexion  displacement 
in  a  fast  movement  and  the  dotted  circles  are  the  predicted  positions. 

As  can  be  seen,  predictions  of  model  are  extremely  accurate. 

Based  upon  results  to  date,  it  seems  allowable  to  say  that  we 
understand  a  great  deal  about  the  major  neuromotor  coordination 
mechanisms  involved  in  fast  human  limb  movement.  The  mathematical 


I 


model  described  previously  can  predict  the  manner  in  which  fast  and  efficient 
limb  movement  can  be  accomplished.  As  mentioned  before,  functional  electrical 
stimulation  (FES)  techniques  can  be  employed  to  stimulate  human  limbs  in 
accordance  with  optimized  models  for  efficient  limb  movement.  The  final 
issue  is  whether  or  not  the  blending  of  mathematical  models  of  neuroraotor 
coordination  mechanisms  with  functional  electrical  stimulation  techniques 
can  produce  real  changes  in  limb  movement  patterns.  That  is,  can  functional 
electrical  stimulation  administered  as  prescribed  by  a  mathematical  model 
of  neuromotor  coordination  produce  any  changes  in  the  neuromotor  coordination 
mechanisms  controlling  limb  movement? 

Motor  Control  and  Learning  Theories 

The  work  of  Lagasse  and  his  associates  have  shown  that  functional 
electrical  stimulation  can  produce  improvements  in  limb  movement  patterns 
without  any  physical  practice.  That  is,  functional  electrical  stimulation 
produced  learning  via  electrical  stimulation  alone.  The  importance  of 
such  "sensory  imparted  learning",  or  SIL,  cannot  be  fully  appreciated  with¬ 
out  some  understanding  of  current  theories  of  motor  control  and  learning. 

Any  attempt  to  present  the  numerous  theories  in  an  adequate  manner  would 
require  a  monumental  effort,  and  is  certainly  not  achievable  in  a  concise 
manner.  Such  an  effort  is  necessary,  however,  to  provide  some  basis  of 
reference  against  which  to  evaluate  the  significance  of  sensory  imparted 
learning  to  motor  control  and  learning  theory. 

At  the  present  time  there  are  numerous  models  and  theories  for  motor 


control  and  learning  ranging  from  simple  descriptive  taxonomies  and 
hierarchical  control  models  to  highly  technical  information  processing 
and  cybernetic  control  models  (Singer,  19S0) .  In  the  context  of  motor 


skill  acquisition, however,  there  are  two  major  categories  into  which  current 
theories  vying  for  acceptance  can  be  classified.  One  category  is  the  open- 
loop  system  of  motor  control  and  learning  in  which  it  is  postulated  that  a 
motor  skill  is  executed  without  being  altered  by  any  feedback  or  error 
detection  mechanisms.  Thus,  a  movement  is  programmed,  and  once  started, 
the  planned  movement  is  unaffected  by  any  peripheral  feedback.  Ballistic 
movements  requiring  less  than  one  normal  reaction  time  for  execution  are 
presumably  under  open-loop  control.  A  second  category  is  the  closed- loop 
system  in  which  feedback  from  the  periphery  can  alter  the  motor  task  during 
its  actual  execution.  The  closed-loop  system  stipulates  mechanisms 
involving  feedback,  error  detection,  and  error  correction  during  the  actual 
execution  of  a  movement  task. 

Open-loop  proponents  hold  that  a  movement  sequence  obeys  a  response 
chaining  principle  in  which  each  serial  aspect  of  the  movement  is  conditioned 
by  the  proprioceptive  feedback  of  the  preceding  movement  segment.  Lash ley 
(1917,  1951)  was  an  early  proponent  of  an  open-loop  theory  and  held  that 
"...sensory  factors  play  a  minor  part  in  regulating  the  intensity  and 
duration  of  nervous  discharge."  The  most  important  feature  in  open-loop 
theory  is  the  central  command  mechanism  while  sensory  feedback  is  viewed 
as  relatively  unimportant.  Indeed,  some  open-loop  proponents  contend  that 
sensory  feedback  is  not  even  necessary.  Taub  and  3erman  (1960),  for  example, 
showed  that  monkeys  could  learn  to  squeete  a  bulb  in  order  to  avoid  a  shock 
stimulus  even  when  the  involved  limb  was  deafferentated.  If  feedback  is 
necessary  for  learning,  it  may  be  provided  by  an  efference  copy  or  corollary 
discharge  from  the  efferent  command  signal. 


Closed-loop  theories  (e.g.,  Adams,  1971;  Anokhin,  1969;  Bernstein, 

1967;  Konorski,  1967;  Sokolov,  1969;  5  Schmidt,  1975),  on  the  other  hand, 
place  greater  importance  on  concurrent  feedback  from  an  ongoing  task. 

Schmidt  (1975),  for  example,  noted  that  "In  order  for  the  subject  to  be 
able  to  receive  information  about  the  correctness  of  the  movement  in 
relation  to  the  desired  outcome,  he  must  be  able  to  compare  the  actual 
feedback  with  the  feedback  expected.”  Xeele  and  Summers  (1976)  likewise 
stated  that  "...the  skill  must  be  executed  for  learning  to  occur,  for 
only  in  that  way  will  kinesthetic  feedback  be  available  for  conditioning 
to  succeeding  movements.”  In  a  paper  proposing  a  model  for  sensorimotor 
control  and  learning,  Raibert  (1978)  stated  that  "...active  movement  is 
essential  to  motor  learning  and  sensorimotor  adaptation.”  Raibert  cited  the 
classical  work  of  Held  (1961)  and  Hein  and  Held  (1965)  in  support  of  his 
position. 

Although  this  treatment  of  motor  control  and  learning  theory  inadequate! 
represents  the  actual  complexity  of  the  topic,  it  does  allow  one  to  point 
out  the  critical  importance  attached  to  the  efferent  signals  from  the  central 
nervous  system  in  the  learning  of  a  motor  task.  Whether  viewed  from  open- 
loop  or  closed- loop  schemas,  active  movement  produced  by  efferent  signals 
programmed  by  central  command  mechanisms  is  considered  to  be  absolutely 
essential  to  motor  skill  acquisition.  Present  motor  learning  theories 
clearly  make  no  provision  for  motor  skill  acquisition  without  efferent 
signals  producing  active  movement. 

Since  sensory  imparted  learning  (SIL)  has  been  demonstrated  several 
times,  it  seems  reasonable  to  suggest  that  functional  electrical  stimulation 
applied  in  accordance  with  a  model  for  neuromctor  coordination  mechanisms 


for  limb  movement  produces  a  sensory  stimulus  pattern  acceptable  to  the 
central  nervous  system  responsible  for  motor  learning.  The  "language"  of 
this  central  nervous  system  center  obviously  includes  the  relevant  parameters 
of  (1)  temporal/sequential  firing  patterns  of  muscles  involved  in  the  move¬ 
ment  synergy,  and  (2)  the  EMG  ratio  representing  intensity  of  agonist  and 
antagonist  muscle  contraction.  Receipt  of  such  relevant  sensory  stimulus 
patterns  is  followed  by  development  of  efferent  motor  programs  which  "match" 
the  sensory  stimulus  patterns  and  eventually  become  available  under  volitional 
control.  Thus,  the  term  sensory  imparted  learning  describes  the  sequence  of 
such  motor  skill  acquisition.  For  convenience,  the  term  "reverse  loop  motor 
learning"  seems  to  have  merit  as  a  description  of  the  sensory  imparted 
learning  phenomenon  in  a  theoretical  context. 

Sensory  Imcarted  Learning 


Several  investigators  reporting  on  the  efficacy  of  functional  electrical 
stimulation  have  noted  relatively  long  lasting  effects  upon  performance 
capabilities  and  have  suggested  such  changes  represent  retraining  or  reorgani 
ration  of  the  central  nervous  system  (Gracanin  and  Marincek,  1970;  Gracanin, 
Vrabic,  and  Vrabic,  1976;  Waters,  McNeal,  and  Perry,  1975;  Vodovnick  and 
Rebersek,  1973).  Such  reports  of  neurological  retraining,  however,  cannot 
be  said  to  be  free  of  criticism.  For  example,  observed  improvements  in 
functional  capacities  could  have  resulted  from  improved  muscle  strength  and 
endurance  rather  than  from  actual  changes  in  neural  control  mechanisms.  In 
all  of  the  studies  reported,  improvement  in  neural  control  mechanisms  cannot 
be  separated  from  improvements  in  muscle  strength  or  endurance  simply  because 


no  assessment  of  neural  control  mechanisms  was  ever  made. 


It  should  be  obvious  that  the  parameters  described  in  the  mathematical 

model  of  neuromotor  coordination  mechanisms  responsible  for  control  of  limb 

movement  represent  a  suitable  set  of  criterion  measures  for  assessing 

reorganisation  of  the  central  nervous  system  or  neurological  retraining. 

Thus,  we  propose  that  any  changes  produced  in  these  model  parameters  brought 

about  by  a  therapeutic  regimen  can  be  interpreted  as  actual  changes  in  the 

central  nervous  system  separate  from  any  changes  produced  in  the  peripheral 

musculature.  Furthermore,  the  question  of  whether  functional  electrical 

stimulation  can  produce  such  changes  has  already  been  answered  in  several 

/ 

crucial  investigations  conducted  by  Pierre  Lagasse  at  Laval  University  m 
Quebec  City,  Canada. 

In  a  series  of  studies,  functional  electrical  stimulation  (FES)  applied 
in  accordance  with  neuromotor  coordination  principles  has  been  shown  capable 
of  producing  improvement  in  several  physical  performance  tasks  without 
benefit  of  actual  physical  practice.  Fleury  and  Lagasse  (1979),  for  example, 
demonstrated  FES  to  be  an  effective  means  for  improving  total  reaction  time 
and  its  central  nervous  system  component  of  premotor  time.  Actual  muscle 
contraction  time  (or  motor  time  in  fractionated  reaction  time  nomenclature) 
was  the  same  in  FES  and  physical  practice  groups  of  subjects.  In  a  subsequen 
study,  Lagasse,  8oucher,  Samson,  and  Jacques  (1979)  applied  electrical 
stimulation  to  three  different  paired  muscle  groups  involved  in  a  weight 
lifting  task,  the  clean  and  jerk  lift.  The  muscle  groups  were  stimulated 
in  a  temporal  sequential  order  equivalent  to  that  exhibited  by  a  highly 
skilled  record  holder  for  the  clean  and  jerk  lift.  Male  subjects  (N  =  12) 
stimulated  with  an  "optimised”  muscle  synergy  pattern  showed  improvement 
without  any  physical  practice,  exhibiting  acquisition  of  a  proper  motor 


pattern  for  execution  of  the  lifting  task.  Boucher  and  Lagasse  (1981)  later 
showed  that  maximum  speed  of  a  horizontal  arm  sweep  task  was  significantly 
improved  by  functional  electrical  stimulation  (FES),  again  based  upon  an 
optimized  pattern  of  neuromotor  coordination  control  mechanisms.  Significant 
changes  in  limb  velocity,  acceleration,  and  temporal  firing  patterns  of  the 
three  involved  paired  muscle  groups  were  shown.  Thus,  functional  electrical 
stimulation--with  no  actual  physical  practice--produced  changes  in  the 
central  nervous  system.  Motor  learning  was  achieved  strictly  through  a 
sensory  pathway. 

Sensory  imparted  learning  (SIL)  through  use  of  functional  electrical 
stimulation  has  been  reported  only  once  before,  and  even  then  some 
skepticism  about  the  claims  was  apparent.  Writing  in  Neural  Organization 
and  its  Relevance  to  Prosthetics  (1973),  Reswick  noted: 

It  was  also  at  about  the  same  time  in  1965  that  L.  A.  Alyeyev 
and  S.  G.  Bounimovich  reported  some  interesting  developments  in 
electrical  stimulation  from  the  Institute  of  Cybernetics  in  Kiev. 

Alyeyev  had  developed  a  multichannel  stimulating  machine  which  he 
called  the  Miotone.  He  claimed  it  operated  on  the  following  principle: 
Myoelectric  signals  of  muscles  used  in  various  functions  previously 
recorded  on  magnetic  tape  obtained  from  'donors’  were  used  to  modulate 
up  to  six  channels  of  an  electrical  stimulating  machine  to  produce 
stimulation  signals  which  were  controlled  by  the  original  EMG  patterns. 

He  claimed  to  be  able  to  produce  thereby  functional  movements  in 
various  types  of  paralysis  in  patients.  These  functional  movements 
produced  by  electricity  not  only  caused  muscle  hypertrophy  and  improved 
metabolic  functions  but  also  assisted  in  the  reorganization  of  the 
central  nervous  system  so  that  stroke  patients  were  able  to  regain 
control  of  various  paralyzed  functions.  To  date  no  controlled  study 
has  been  reported  to  confirm  these  assertions. 

It  seems  allowable  to  suggest  that  Alyeyev  and  Bounimovich  (1965)  applied 
functional  electrical  stimulation  to  hemiplegic  limbs  in  accordance  with  a 
crude  model  of  neuromotor  coordination  mechanisms,  using  EMG  patterns  from 
The  series  of  studies  by  Lagasse  and  his  associates  previously 


"donors" . 


cited  used  functional  electrical  stimulation  based  upon  an  improved  model 

of  neuromotor  coordination  mechanisms  and  successfully  produced  improvement 

/ 

without  physical  practice  in  several  performance  tasks.  Lagasse  and  his 
associates,  however,  were  using  only  an  improved  model  of  neuromotor 
coordination  mechanisms,  and  not  the  best  model.  Although  the  neuromotor 
coordination  mechanism  of  temporal  sequence  of  muscle  firing  was  controlled, 
the  second  major  mechanism  involving  agonist  to  antagonist  EMG  duration  and 
intensity  was  not  described  until  quite  recently  (Kilmer,  Kroll,  and 
Congdon,  1981).  Thus,  if  functional  electrical  stimulation  is  applied 
using  the  full  and  most  modem  model  of  neuromotor  coordination  mechanisms 
for  control  of  limb  movements,  even  greater  effects  can  be  expected. 

Summary 

It  has  been  proposed  that  the  implicit  goal  of  all  neurophysiologicaily 
based  therapeutic  exercise  systems  is  to  transmit  sensory  stimulus  patterns 
to  the  central  nervous  system  which  are  more  intense  than  that  produced 
without  benefit  of  proprioceptive  facilitation.  Implicit  in  such  a  treat¬ 
ment  rationale  is  the  belief  that  transmission  of  heightened  sensory  stimulus 
patterns  to  the  central  nervous  system  will  enhance  recovery  of  lost 
neurological  functions.  Various  mechanisms  for  neurological  function  recovery 
have  been  proposed:  collateral  sprouting,  unmasking,  reactive  synapto genes is , 
sublimal  fringe  neuron  activation,  reclaiming  of  ipsilateral  nerve  pathway 
control  tracts. 

Motor  command  centers  in  the  brain  which  have  been  damaged  by  some 
pathology  can  produce  only  inefficient  and  incorrect  innervation  commands 


to  the  muscles.  Limbs  afflicted  by  some  pathological  condition  can  produce 
only  inefficient  and  incorrect  volitional  movement  patterns.  Such  inefficient 
and  incorrect  volitional  movement  produces  a  sensory  stimulus  pattern  which 
is  itself  inefficient  and  incorrect,  even  if  enhanced  by  proprioceptive 
facilitation  techniques.  If  recovery  of  central  nervous  system  functioning 
is  dependent  to  any  degree  upon  the  quality  of  the  sensory  stimulus  pattern 
received,  the  recovery  is  likely  to  be  slow  and  dependent  upon  gradual 
improvement  in  the  quality  of  the  limb  movement  synergy  which  produces 
slight  sensory  stimulus  patterns. 

If,  however,  the  sensory  stimulus  pattern  transmitted  to  the  central 
nervous  system  from  an  afflicted  limb  was  not  an  impoverished  pattern, 
but  a  sensory  stimulus  pattern  representing  optimally  efficient  limb  movement, 
the  likelihood  of  improved  neurological  recovery  might  be  enhanced.  Such 
a  possibility  now  seems  feasible  because  of  developments  in  the  areas  of 
functional  electrical  stimulation,  neuromotor  coordination  mechanisms 
controlling  fast  limb  movements,  and  sensory  imparted  learning.  These 
relatively  distinct  areas  of  research  have  an  obvious  relevance  to  each 
other,  and  a  blending  of  their  complementary  components  holds  promise  for 
a  vastly  improved  rehabilitation  regimen  for  neurological  recovery.  Based 
upon  a  mathematical  model  of  neuromotor  coordination  mechanisms  for  optimally 
efficient  limb  movement,  functional  electrical  stimulation  can  be  applied 
to  an  afflicted  limb  which  can  produce  an  optimal  sensory  stimulus  pattern 
for  transmission  to  the  central  nervous  system.  If  neurological  recovery 
is  dependent  to  any  degree  upon  the  quality  of  the  sensory  stimulus  pattern, 
receipt  of  such  optimized  sensory  stimulus  patterns  should  be  more  effective 
than  any  present  therapeutic  exercise  regimen.  Since  the  functional  electrica 
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ABSTRACT 


Twenty  male  and  twenty  female  subjects  (ages  17  to  28  years)  were 
measured  for  maximal  isometric  elbow  flexion  and  extension  strength, 
upper  limb  volume,  lengths  and  girths.  Female  flexion  strength  was 
54. 3“;  and  extension  strength  S8.2%  of  males.  Covariance  adjustment  for 
body  weight  and  upper  limb  volume  increased  female  flexion  strength  to 
”8%  and  extension  strength  to  73%  of  males,  but  males  were  still  signi¬ 
ficantly  stronger.  Arm  girths  and  limb  volume  exhibited  larger  sex 
differences  than  strength  measures.  Limb  lengths  were  not  as  important 
as  girths;  and  body  weight  and  limb  volume  importance  in  strength  pre¬ 
diction  differed  in  males  and  females.  Multiple  R's  of  .84  to  .94  for 
strength  resulted  from  a  full  set  of  anthropometric  predictors.  Simple 
measures  of  limb  girths  and  lengths,  with  and  without  body  weight,  were 
as  good  or  better  than  segmental  limb  volumes  in  strength  prediction. 

It  was  clear  that:  (1)  the  importance  of  anthropometric  measures  as 
strength  predictors  differs  not  only  between  the  sexes  but  between 
strength  measures  as  well;  (2)  male  and  female  differences  extend  beyond 
absolute  differences  in  anthropometric  measures  and  includes  differences 
in  the  relationship  of  upper  limb  anthropometric  measures;  (5)  body  size 
factors  do  not  explain  adequately  the  strength  superiority  of  males. 


Composite  strength  indices  based  upon  representative  batteries  of  strength 


measures  have  generally  shown  moderate  but  significant  correlations  with 
various  physique  measures  (Jones,  1947;  Laubach  §  McConville,  1966,  1969; 

Sills  f,  Everett,  1953).  The  most  common  anthropometric  measure  assessed  in 
conjunction  with  strength  is  undoubtedly  body  weight  followed  in  frequency 
of  usage  by  limb  girths  and  lengths.  Because  of  the  accepted  relationship 
between  muscle  cross-sectional  area  and  strength,  the  use  of  anthropometric 
measures  in  strength  studies  seems  logical.  Unfortunately,  much  of  the 
research  dealing  with  the  relationship  of  strength  and  anthropometric  measures 
has  been  equivocal  with  many  investigators  questioning  thr  usefulness  of 
anthropometric  measures  in  strength  prediction  (Cureton,  1947;  Rasch  § 

Pierson,  1963;  Smith  §  Rovce,  1963). 

Only  a  few  studies  have  included  both  males  and  females  making  it  difficult 
to  compare  the  importance  of  anthropometric  predictors  of  strength  between  the 
sexes.  Even  fewer  studies  have  included  segmental  limb  volumes  specific  to 
the  muscle  group  being  studied.  Due  to  provisions  of  the  Equal  Employment 
Opportunities  Act,  the  increased  presence  of  females  in  the  military,  and 
the  growth  of  organized  athletics  for  females,  renewed  attention  has  been 
given  to  differences  and  similarities  of  strength  and  related  parameters  in 
females  and  males.  In  the  present  study,  equal  numbers  of  males  (N=20)  and 
females  (\'=20)  were  measured  for  maximal  isometric  strength  of  elbow  flexion 
and  extension.  Strength  measures  were  then  correlated  by  sex  with  body  weight, 
limb  lengths  and  girths,  and  segmental  upper  limb  volumes.  As  will  be  shown, 
the  importance  and  effectiveness  of  anthropometric  measures  as  strength 
predictors  differs  not  only  between  the  sexes  but  between  strength  measures 


as  well. 


MATERIALS  AND  METHODS 


Twenty  male  and  twenty  female  subjects  (ages  17  to  28  years)  were  given 
medical  clearance  and  signed  informed  consent  forms  as  prescribed  by  university 
guidelines  for  the  protection  of  human  subjects.  Each  subject  completed  three 
data  collection  sessions  within  a  one-week  period.  Half  of  the  subjects  con¬ 
tinued  on  in  a  longer  study  of  arm  movement  speed  to  be  reported  elsewhere. 

Anthropometric  measures.  In  addition  to  body  weight,  segmental  limb 
volumes  for  the  right  upper  arm,  forearm,  and  hand  were  determined  using  a 
water  displacement  technique  (Dempster,  1955;  Katch  §  Weltman,  1975).  Seg¬ 
mental  volumes  were  also  estimated  by  an  anthropometric  technique  (Katch  5 
Weltman,  1975;  Sady  et  al . ,  1978).  Measures  of  limb  lengths  and  girths 
required  for  volume  estimation,  and  employed  as  separate  predictors  of  strength, 
wore  as  follows. 

Lengths 

LI  acromion  process  to  deltoid  tubercle 
L2  deltoid  tubercle  to  olecranon  process 
L3  olecranon  process  to  styloid  process  of  ulna 

L4  styloid  process  of  ulna  to  tip  of  third  finger 

Girths 

AC  circumference  at  acromion  process 

DEL  circumference  at  deltoid  tubercle 

ELB  circumference  at  olecranon  process 

WJS  circumference  at  distal  space  to  styloid  process  of  ulna 

(wrist  joint  space) 

HND  thickness  of  base  of  hand,  cross-section  height  thenar  eminence 
and  hypothenar  eminence 

I sometric  strength.  Four  5-sec  maximal  isometric  contractions  were 
■ncured  for  right  elbow  flexion  and  extension  strength.  Each  trial  was 
separated  by  a  1-minute  rest  interval  and  the  flexion-extension  series  by 
a  7-10  minute  rest  interval.  The  order  of  flexion  and  extension  testing  was 


balanced  over  subjects  with  no  differences  in  mean  strength  due  to  order  of 
testing  being  detected.  Two  trials  were  maximal  voluntary  contractions  (MVC) 
measures  where  subjects  were  instructed  to  exert  a  maximal  contraction  for 
5-seconds.  The  other  two  trials  were  fast  MVC  (FMVC)  measures  where  subjects 
were  instructed  to  reach  maximum  tension  as  quickly  as  possible  and  maintain 
maximum  until  commanded  to  stop.  Peak  tension  output  of  each  trial  was  taken 
as  maximum  strength. 

Subjects  were  tested  in  a  sitting  position,  the  chest  braced  against 
a  padded  table  support,  with  a  wide  nylon  belt  around  the  upper  body  to 
prevent  movement  away  from  the  test  apparatus.  The  upper  arm  was  parallel  to, 
and  supported  by,  the  table  top  at  a  90  degree  angle  to  the  shoulder.  The 
forearm  was  positioned  at  angles  of  75  degrees  for  flexion  and  at  90  degrees 
for  extension.  In  accordance  with  Basmajian  and  Latif  (1957),  the  elbow 
angles  are  given  as  the  compliment  of  the  angle  between  the  forearm  and  the 
upper  arm.  A  wrist  cuff  connected  via  an  adjustable  steel  cable  to  a  Statham 
strain  gauge  which  monitored  tension  output  and  was  recorded  on  a  Beckman  Type 
R  Dynagrapii.  The  angle  of  pull  for  all  trials  was  at  90  degrees  to  the  strain 
gauge  with  the  hand  positioned  midway  between  supination  and  pronation.  Fore¬ 
arm  position  is  known  to  affect  the  expression  of  elbow  strength  but  the  use 
of  a  wrist  cuff  effectively  eliminates  such  effects  (Provins  5  Salter,  1955). 
Alteration  in  forearm  position  from  the  midline  during  maximal  efforts  would 
thus  not  affect  the  expression  of  maximum  strength. 


TABLE  1  TO  APPEAR  ABOUT  HERE 
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RESULTS 

As  seen  in  Table  1,  maximum  isometric  elbow  flexion  strength  was  higher 
than  extension  strength  by  20%  in  males  and  12.1%  in  females.  The  superiority 
of  elbow  flexion  strength  is  a  common  finding  in  males  (Knapik  8  Ramos,  1980; 
Singh  8  Karpovich,  1968),  in  females  (Kroll,  1972),  and  in  elite,  trained 
weight  lifters  and  middle  distance  runners  (Tornvall,  1963).  Higher  elbow 
flexion  strength  prevails  even  though  the  brachial  triceps  typically  possesses 
a  higher  percentage  of  fast-twitch  muscle  fibers  (Buchthal  5  Schmalbruch,  1970; 
■Johnson  et  a!.,  1973),  and  has  a  larger  physiological  cross-sectional  area 
An  et  al.,  1981).  One  factor  enabling  the  elbow  flexors  to  produce  higher 
strength  is  a  decided  mechanical  advantage.  The  major  flexors  of  the  forearm 
biceps  brachii,  brachialis,  and  brachioradialis)  have  force  arms  two  to  four 
tines  longer  than  the  brachial  triceps  (Ikai  8  Fukunaga,  1968;  Morris,  1948; 
Wilkie,  1950). 

Males  were  significantly  stronger  than  females  on  both  elbow  flexion  and 
elbow  extension  isometric  strength  as  other  investigators  have  reported 
, Petrofsky  5  Phillips,  1980;  Singh  8  Karpovich,  1968).  Female  flexion  strength 
was  54.3%  and  extension  strength  58.2%  when  expressed  as  a  percentage  of  male 
strength.  Hettinger  (1961)  reported  female  elbow  flexion  and  extension 
strength  as  being  55%  of  males,  which  is  comparable  to  the  present  results, 
while  Singh  and  Karpovich  (1968)  reported  smaller  values  of  48%  (extension) 
and  44%  i flexion) .  Nordgren  (1972)  found  female  elbow  flexion  strength  to  be 
54%  and  extension  strength  52%  of  male  strength. 

There  were  no  significant  differences  between  MVC  and  FMVC  strength  in 
either  mates  or  females.  Observed  differences  were  less  than  one  Kg  with 


males  exhibiting  higher  MVC  and  females  higher  FMVC  values.  Caldwell  et  U  . 
i 1974)  and  Rupp  and  Baranowski  (1934)  also  reported  no  difference  between 
slew  and  fast  expression  of  maximum  grip  strength.  In  elderly  males,  however, 
a  fast  maximal  voluntary  contraction  produces  a  si gnificant 1 y  higher  isometric 
strength  by  some  19%  compared  with  a  gradual  expression  of  isometric  knee 
extension  strength  (Clarkson  et  al.,  1981).  MVC  and  FMVC  strength  correlated 
highly  in  males  and  females  with  r's  ranging  from  .88  to  .98.  Since  there 
sere  no  differences  between  MVC  and  FMVC  absolute  values,  and  since  the 
measures  were  so  highly  correlated,  all  further  data  analysis  reports  only  the 
MVC  results. 


Anthropometric  measures.  As  seen  in  Table  2,  males  and  females  differed 
significantly  on  all  anthropometric  measures  except  LI  (acromion  process  to 
deltoid  tubercle)  and  L3  (olecranon  process  to  styloid  process  of  ulna).  Males 
were  14.1  kg  or  23.5%  heavier  than  females  in  body  weight  while  total  upper 
limb  volume  was  1,191.8  ml  or  43.9%  larger  in  males.  Differences  in  arm  girths 
ranged  from  13.6%  to  18.2%  and  in  limb  lengths  from  5.25%  to  10.9%  in  favor  of 
the  males.  Thus,  arm  girths  and  limb  volume  exhibited  larger  sex  differences 
than  limb  strength  measures. 

In  a  study  involving  41  adult  males  Roberts  et  al.  (1959)  reported  that 
upper  arm  girth  correlated  negatively  with  four  longitudinal  measures  (height, 
upper  arm  length,  forearm  length,  and  hand  length).  The  same  general  relation¬ 
ship  is  seen  for  males  in  the  present  study  (see  Table  4)  with  AC  girth 
(circumference  at  acromion  process)  and  DEL  girth  (circumference  at  deltoid 


Table  2 


>  Means  for  Body  Weight  (kg).  Girths  (cm).  Length  (cm),  and  Limb  Volume  (ml) 


Male 

Female 

(N=20) 

(N=20) 

Mean 

SD 

Mean 

Body  Weight 

74.0 

6.1 

59.9 

7.8 

Girths 

AC' 

41.0 

2.8 

34.7 

2.5 

DEL* 

29.6 

2.4 

25.9 

2.0 

ELB* 

27.8 

1.3 

23.7 

1.2 

W'JS  * 

17.6 

0.6 

15.4 

0.6 

HMD* 

5.0 

0.3 

4.4 

0.5 

Length 

Li 

14.4 

2.1 

13.5 

1.5 

i.2* 

15.2 

1.0 

13.7 

1.1 

Li 

24.3 

2.3 

23.1 

1.6 

L4* 

21.0 

0.9 

19.3 

0.9 

Limb  Volumes 

SEG1* 

1443. 3 

292.1 

1012.6 

224.3 

.SEG2' 

1003.5 

160.4 

682.4 

121.1 

SEG3* 

1034.5 

106.9 

718.1 

106.5 

SEG4  * 

428.4 

61.5 

304.8 

41.0 

Lota  i  * 

3909 .  - 

479.1 

2717.9 

435.9 

"Significant  difference  between  males  and  females,  p  <  .01 


tubercle)  each  failing  to  correlate  significantly  with  any  of  the  four  limb 
length  measures.  In  females,  however,  five  out  of  the  eight  possible  correla¬ 
tions  between  these  two  girth  measures  and  the  four  limb  length  measures  were 
significant.  Differences  between  males  and  females  thus  extend  beyond  absolute 
differences  in  anthropometric  measures  and  includes  differences  in  the  relation¬ 
ship  of  anthropometric  measures  in  the  upper  limb  as  well. 


TABLE  3  ABOUT  HERE 


Observed  strength  differences  between  men  and  women  are  presumed  to  be 
strongly  influenced  by  over-all  body  size  rather  than  by  muscle  quality. 
Comparison  of  male  and  female  strength,  therefore,  is  more  prudent  if  body  size 
factors  can  be  controlled.  In  the  present  study,  both  body  weight  and  total 
upper  limb  volume  demonstrated  significant  relationships  with  strength  but  the 
relationship  was  different  in  the  two  sexes.  A  covariance  analysis  was  there¬ 
fore  made  to  control  for  tile  influence  of  both  body  weight  and  upper  limb  volume 
upon  strength..  The  adjusted  means  for  male  and  female  elbow  flexion  and  exten¬ 
sion  strength  are  presented  in  Table  3.  The  effect  of  covariance  adjustment  for 
body  size  is  best  observed  when  female  strength  is  expressed  as  a  percentage 
of  naie  strength.  Initiallv  female  flexion  and  extension  strength  was  54.3% 
and  58.2%  of  male  strength.  Adjusted  for  body  weight  and  upper  limb  volume, 
female  flexion  and  extension  strength  increased  to  78%  and  73%  of  male  strength. 


In  agreement  with  others  (Hoffman  et  al.,  1979;  Morrow  &  Hosier,  1981),  however, 
male  strength  was  still  superior  even  after  covariance  adjustment  for  body 
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Table  3 


Adjusted  Means  for  Flexion 

and  Extension  Strength  in  Males  and  Females 

Male 

Female 

%  Strength 

(N=20) 

(N=20) 

(Female/Male) 

Flexion 


Raw 

47.00 

25.50 

54% 

Adjusted 

40.51 

31.90 

78% 

Extension 

Raw 

39.10 

22.70 

58% 

Adjusted 


35.60 


26.30 


73% 


[II 


TABLE  4  ABOUT  HERE 


Inter-correlation  matrix.  Total  upper  limb  volume  in  females  correlated 
significantly  with  all  other  variables  (see  Table  4).  In  males,  however, 
total  upper  limb  volume  failed  to  correlate  with  body  weight  (r=.35,  ns),  hand 
girth  (r=.24,  ns),  length  3-forearm  (r=-.04,  ns),  or  with  any  of  the  strength 
measures  (r's  from  .05  to  .34,  ns).  Body  weight  correlated  significantly  with 
MVC  flexion  strength  in  males  (r=.69,  p  <  .01)  and  in  females  (r=.59,  p  <  .01) 
but  failed  to  correlate  with  MVC  extension  strength  in  either  sex.  Girth 
measures  had  more  and  higher  correlations  with  strength  measures  in  females  than 
in  males.  Of  the  eight  possible  correlation  coefficients  for  limb  length  with 
strength  measures,  two  correlations  were  significant  in  females  (both  with 
flexion  strength),  and  only  one  correlation  was  significant  in  males.  Limb 
length  3  (forearm)  and  length  4  (hand)  correlated  with  MVC  flexion  strength 
fr=.00  and  .“1)  in  females  while  in  males  these  limb  lengths  showed  no  relation¬ 
ship  with  strength.  The  general  impression  given  by  the  correlation  matrix 
suggests  that  limb  lengths  are  not  as  important  as  limb  girths,  and  the  importance 
>f  body  weight  and  limb  volume  differs  between  males  and  females  in  the  prediction 
of  elbow  flexion  and  extension  strength. 

Mu  1 1 i p 1 c  correlation  analysis .  A  better  understanding  of  the  relative 
importance  of  anthropometric  measures  in  the  prediction  of  isometric  elbow 
flexion  and  extension  strength  can  be  realized  by  multiple  correlation  analysis. 
Using  stepwise  multiple  correlation,  analysis  was  made  in  a  sequential  stage- 
si  sc  fashion  with  each  stage  including  a  different  set  or  combination  of  anthropo¬ 
metric  predictors.  In  the  first  stage,  body  weight  (BWT)  and  total  upper  limb 


Values  below  diagonal  for  males;  values  above  diagonal 


volume  (VTOT)  were  the  only  predictors.  Next,  body  weight  and  segmental  limb 

volumes  were  used  as  predictors  followed  by  segmental  limb  volumes  and  total 

upper  limb  volume.  Limb  lengths  and  girths,  with  and  without  body  weight  and 

segmental  li'iib  volumes,  were  used  as  predictor  sets  to  contrast  the  predictive 

efficiency  of  simple  length  and  girth  measures  with  limb  volume  and  body  weight 

measures.  The  final  analysis  used  the  full  set  of  anthropometric  measures  to 

predict  elbow  flexion  and  extension  MVC  strength  in  males  and  in  females. 

Results  of  the  multiple  correlation  analyses  are  presented  in  Table  5.  The 

2 

decision  rule  for  the  inclusion  of  a  predictor  was  that  it  must  increase  R  by 

2 

at  least  five  percent  or  that  deletion  must  not  decrease  R  by  more  than  five 
percent . 


TABLE  5  ABOUT  HERE 


Limb  volume  would  logically  be  expected  to  correlate  more  highly  with  limb 
strength  than  the  more  global  anthropometric  measure  of  body  weight.  Given  the 
two  predictors  of  body  weight  (BWT)  and  total  upper  limb  volume  (VTOT),  however, 
BWT  is  a  better  predictor  of  male  flexion  strength  (r=.69,  p  <.01)  and  extension 
strength  (r=.24,  ns).  In  females,  however,  VTOT  is  a  better  predictor  than  BWT 
for  flexion  strength  (r=.72,  p  < ,01)  and  extension  strength  (r=.63,  p  <  .01). 
A!i  subsequent  multiple  correlation  analyses  produced  higher  R's  than  the  above 
which  used  only  BWT  and  VTOT. 

Multiple  correlations  using  BWT  plus  segmental  limb  volumes  were  as  good 
or  better  than  prediction  equations  employing  segmental  limb  volumes  plus  VTOT. 
in  both  sets  of  equations  elbow  flexion  strength  was  better  predicted  than 


vxtension  strength.  Only  limb  volume  segment  2  (deltoid  tubercle  to  elbow) 
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Table  5 


Multiple  Correlation  Analysis* 


Body  Weight  (BWT)  *  Volume  Total  (VTOT) 


Flexion  BWT  (47.6%) 
Extension  BWT  (5.S%) 


Bernal e 


Flexion  VTOT  (51.5%) 

Extension  VTOT  (36.1%)  +  BWT  (36%) 


Blv;  +  Limb  Volumes 


Flexion  BWT  (47.4%)  +  SEG2  (10.1%)  +  SEG3  (7.7%) 
Extension  SEG2  (13.7%)  +  SEG3  (20.2%) 


Female 


Flexion  SEG3  (60.1%) 

Extension  SEG3  (45.1%)  +  BWT  (11.0%) 


Limb  Volumes  *  VTOT 


Flexion  SEG2  (44.1%)  +  SEG3  (10.3%) 
Extension  SEG2  (13.7%)  +  SEG3  (20.2%) 


Femah 


Flexion  SEG3  (60.1%) 
Extension  SEG3  (45.1%) 


r.'-"  ■■  1  <«■  «•  l»  *  ■  r»  L»  ?■  J  *  l>  ,' !■.' I*1 ".!■■«.> 

» 


% 


Limb  Lengths  +  Limb  Girths  R 

Males 

Flexion  DEL-G  (46.3)  +  LI  (6.8)  +  L4  (8.8)  .79 

Extension  DEL-G  (31.4)  +  ELB-G  (11.8)  +  HND-G  (13.0)  + 

AC-G  (13.1)  .83 

Females 

Flexion  ELB-G  (53.1)  +  L4  (16.4)  .83 

Extension  ELB-G  (64.6)  +  WJS-G  (4.8)  +  HND-G  (5.5)  .87 


Limb  Lengths  +  Limb  Girths  +  BWT 
Males 

Flexion  BWT  (47.4)  +  DEL-G  (13.8)  +  WJS-G  (4.6)  + 


Li  (6.0)  +  AC-G  (5.3)  +  HND-G  (4.1)  .90 

Extension  DEL-G  (31.4)  +  ELB-G  (11.8)  +  HND-G  (13.0)  + 

AC-G  (13.1)  .83 

Females 

Flexion  ELB-G  (53.1)  +  L4  (16.4)  .83 

Extension  ELB-G  (64.6)  +  BWT  (5.5)  +  HND-G  (4.8)  .87 

Limb  Lengths  +  Limb  Girths  +  Segmental  Volumes  +  VTOT 

Ma  1  es 


Flexion  DEL-G  (46.3)  +  LI  (6.8)  +  L4  (8.8)  +  SEG3  (3.1)  + 


1.2  (7.4)  .85 

Extension  DEL-G  (31.4)  +  SEG3  (18.2)  +  AC-G  (13.0)  + 

SEG4  (7.7)  .84 

Females 

Flexion  SEG3  (60.1)  +  ELB-G  (3.1)  +  WJS-G  (4.0)  + 

LI  (2.2)  *  AC-G  (5.2)  +  HND-G  (7.8)  .94 

Extension  ELB-G  (64.6)  +  WJS-G  (4.8)  +  HND-G  (5.5)  .S7 


A.V. 


■‘J-* 


*  *A  M 


All  Predictors 


Males 


Flexion  BWT  (47.4)  +  DEL-G  (13.3)  +  WJS-G  (4.6)  + 
L4  (5.9)  +  SEG1  (5.5) 

Extension  DEL-G  (31.4)  +  SEG3  (18.2)  +  AC-G  (13.0) 
SEG4  (7.7) 

Females 

Flexion  SEG3  (60.1)  +  ELB-G  (3.1)  +  WJS-G  (4.0)  + 

LI  (2.2)  +  AC-G  (5.2)  +  HND-G  (7.7) 

Extension  ELB-G  (64.5%)  +  BWT  (5.5%)  +  HND-G  (4.8%) 


'Values  in  parentheses 


indicate  contribution  to  R 


2 
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and  limb  volume  segment  3  (forearm)  were  selected  from  the  volume  predictors. 
BWT  was  added  to  SEG  2  and/or  SEG  3  predictor  sets  raising  the  multiple  R 
from  .74  to  .81  for  male  flexion  strength  and  from  .67  to  .75  for  female 
extension  strength. 

Limb  length  and  girth  predictor  sets  produced  generally  higher  multiple 
R's  than  predictor  sets  employing  segmental  limb  volumes,  VTOT,  and  BWT.  All 
extension  strength  predictors  in  males  and  females  were  girth  measures  while 
prediction  equations  for  flexion  strength  had  fewer  predictors  and  contained 
both  limb  length  and  limb  girth  measures.  Unlike  prediction  equations  using 
segmental  limb  volumes,  prediction  equations  using  limb  lengths  and  girths 
produced  higher  multiple  R's  for  extension  than  for  flexion  strength.  The 
addition  of  BWT  to  limb  length  and  girth  predictor  sets  improved  the  multiple 
R  only  for  male  flexion  strength,  r  increases  from  .79  to  .90. 

When  all  predictors  are  included  in  the  multiple  R  analysis  (BWT,  limb 
lengths  and  girths,  segmental  limb  volumes,  and  total  upper  limb  volume),  the 
multiple  R  for  male  flexion  strength  containing  five  predictors  of  .88  is 
actually  less  than  the  mutiple  R  of  .90  based  upon  six  predictors  without 
segmental  limb  volumes.  For  male  extension  strength,  a  multiple  R  of  .84 
with  four  predictors  including  segment  limb  volumes  is  quite  similar  to  the 
multiple  R  of  .83  based  upon  four  limb  girth  predictors.  A  similar  situation 
exists  for  the  females.  Regardless  of  whether  segment  limb  volumes  are 
included  in  the  regression  analysis,  the  same  three  limb  girth  measures  are 
chosen  for  extension  strength  prediction.  A  six  measure  prediction  equation 
including  segmental  limb  volumes  produces  a  multiple  R  of  .94  for  female 
flexion  strength  while  a  six  measure  prediction  equation  including  only  limb 
lengths,  girths,  and  BWT  produces  a  multiple  R  of  .90. 
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DISCUSSION 


The  present  results  agree  with  many  previous  reports  (Borchardt,  19b8; 
Pctrofsky  5  Lind,  1975;  Rasch  &  Pierson,  1963;  Royce,  1958;  Smith  6  Royce, 

1963)  which  have  failed  to  show  a  consistent  relationship  between  body  weight 
and  isometric  strength.  The  present  study  showed  a  significant  correlation 
between  BWT  and  MVC  flexion  strength  in  both  males  (r=.69,  p<.01)  and  females 
(r=.S9,  p  <.01)  but  BWT  failed  to  correlate  with  MVC  elbow  extension  strength 
in  either  sex.  Body  weight  correlated  significantly  with  total  upper  limb 
volume  (VTOT)  in  females  (r=.77,  p<.01)  but  not  in  males  (r=.3S,  ns).  VTOT 
correlated  significantly  with  female  elbow  flexion  strength  (r=.72,  p<.01) 
and  female  extension  strength  (r=.62,  p<.01),  but  not  with  male  strength 
measures.  Thus,  the  relative  importance  of  body  weight  and  total  upper  limb 
volume  as  strength  predictors  differed  between  males  and  females  and  between 
elbow  flexion  and  extension  strength. 

In  agreement  with  others  (Laubach  §  McConville,  1966;  Tomvall,  1963), 
limb  girth  measures  were  better  than  limb  lengths  as  strength  predictors. 

•h  wever,  this  generalization  requires  some  qualification.  Although  it  is 
true  that  no  limb  length  measures  appeared  in  any  prediction  equations  for 
elbow  extension  strength,  several  limb  length  measures  did  appear  in  elbow 
flexion  strength  prediction  equations.  The  appearance  of  L4  (hand  length) 
as  a  significant  predictor  of  flexion  strength  in  several  equations  is  of 
interest.  Roberts  et  al .  (1959)  also  found  that  hand  length  correlated 
significantly  with  male  flexion  strength  (r=.59)  and  extension  strength  (r=.42) 
Hand  length  (L4)  appeared  in  six  of  the  16  multiple  correlation  analyses  using 
limb  girths  and  lengths  as  part  of  the  predictor  set.  In  every  instance 


hand  lergth  appeared  in  prediction  of  flexion  strength  (four  times  for  males, 


twice  for  females),  but  never  in  prediction  equations  for  extension  strength. 
Since  all  strength  measures  were  secured  using  a  wrist  cuff,  the  hand  could 
not  participate  in,  or  contribute  to,  the  actual  expression  of  strength. 

The  importance  of  hand  length  in  elbow  flexion  strength  prediction 
equations  could  reflect  some  inherent  feature  of  skeletal  size,  lever  length, 
or  somatotype  related  to  isometric  strength.  The  fact  that  hand  length  did 
not  appear  in  any  extension  strength  prediction  equations,  however,  weakens 
support  for  such  a  speculation.  In  addition,  body  weight  correlated  signi¬ 
ficantly  only  with  flexion  strength  in  both  males  and  females.  Coupled  with 
the  finding  of  flexion  over  extension  strength  superiority  in  both  sexes, 
the  importance  of  hand  length  may  be  related  to  a  recent  observation  by 
Goslin  and  Charteris  (1982)  concerning  the  arboreal  heritage  of  humans. 
According  to  these  authors,  flexion  strength  superiority  would  be  logical  in 
arboreal  locomotion  where  the  elbow  is  chiefly  used  for  suspension.  On  similar 
phylogenetic  grounds,  flexion  strength  superiority  would  be  advantageously 
served  by  a  longer  hand  for  climbing  and  arboreal  locomotion. 

The  absence  of  limb  length  measures  in  elbow  extension  strength  prediction 
equations  was  also  reported  by  Roberts  et  al.  (1959)  and  may  be  due  to  ana¬ 
tomical  features  of  the  triceps  brachii  muscle.  It  is  known  that  the  three 
triceps  brachii  components  have  a  larger  volume  but  a  shorter  fiber  length 
as  compared  to  the  biceps  brachii  muscle  (An  et  al.  1981).  Being  a  bipennate 
muscle  with  shorter  and  more  obliquely  running  fibers,  the  triceps  girth 
measure  would  logically  be  more  important  than  length  measures.  The  biceps 
brachii  muscle,  on  the  other  hand,  is  a  fusiform  muscle  with  its  fibers 
x t e:ui i  ng  tac  entire  muscle  length,  and  one  might  expect  limb  length  measures 
t  '  be  of  more  importance  in  the  expression  of  elbow  flexion  strength. 
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In  a  study  involving  41  adult  males,  Roberts  et  al.  (1959)  found  that  the 
multiple  R's  were  slightly  but  consistently  lower  for  prediction  of  elbow 
flexion  than  for  elbow  extension  strength.  No  such  consistent  pattern  in  the 
present  study  was  found  although  prediction  equations  for  elbow  flexion  strength 
generally  contained  a  greater  number  of  predictor  items  than  extension  strength. 
Although  Roberts  etal.  (1959)  measured  height,  weight,  three  limb  lengths, 
and  two  limb  girths,  their  multiple  correlation  analyses  reported  only  predictor 
sets  of  three  measures;  two  of  the  three  measures  were  always  height  and  weight. 
The  highest  multiple  R's  they  reported  for  prediction  of  elbow  flexion  ( R= .73) 
and  extension  strength  (R=.76)  were  lower  than  all  multiple  R's  obtained  in 
the  present  study  which  included  body  weight  and  limb  lengths  and  girths  as 
part  of  the  predictor  set.  The  multiple  R's  in  the  present  study,  ranging 
from  .84  to  .94,  for  the  full  predictor  set  analyses,  are  the  highest  multiple 
R's  reported  to  date  for  elbow  flexion  and  extension  strength  using  anthropo¬ 
metric  measures  as  predictors. 

I 

Based  upon  all  analyses  conducted,  the  present  results  strongly  support 
the  contention  that,  compared  to  segmental  limb  volumes,  limb  lengths  and 
girths  are  better  predictors  of  elbow  flexion  and  extension  strength  in  both 
male  and  female  subjects.  The  difference  in  multiple  R's  obtained  with  and 
without  segmental  limb  volumes  are  small  and  not  significantly  different. 
Anthropometric  assessment  techniques  for  estimating  segmental  limb  volumes, 
moreover,  require  all  of  the  anthropometric  limb  girth  and  length  measures 
reported  here  to  be  obtained.  As  far  as  the  results  of  this  study  are  con¬ 
cerned,  simple  measures  of  limb  girths  and  lengths,  with  and  without  body 
weight,  are  as  good  or  better  than  segmental  limb  volume  measures  in  the 
prediction  of  elbow  strength  in  females  and  males.  It  is  also  rather  clear 
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